SPECTRAL PROPERTIES OF POLYHARMONIC OPERATORS WITH 
LIMIT-PERIODIC POTENTIAL IN DIMENSION TWO. 



YULIA KARPESHINA 1 AND YOUNG-RAN LEE 1 

o i 

In memory of our colleague and friend 
Robert M. Kauffman. 



(N 



5 , 

^ : 

Abstract. We consider a polyharmonic operator H = (— A) +V(x) in dimension two with / > 6 
and a limit-periodic potential V(x). We prove that the spectrum of H contains a semiaxis and 
qq ■ there is a family of generalized eigenfunctions at every point of this semiaxis with the following 

properties. First, the eigenfunctions are close to plane waves e^ k ' x ^ at the high energy region. 
Second, the isoenergetic curves in the space of momenta k corresponding to these eigenfunctions 
g ; have a form of slightly distorted circles with holes (Cantor type structure). 

vo 

O ■ 1. Introduction 

r- 1 

We study an operator 

JU\ H = (-A) l + V(x) (1.1) 

in two dimensions, where I > 6, V(x) is a limit-periodic potential: 

a' 
> 



x 

5-H 



V{x)=Y,V r {x), (1.2) 

r=l 

{T^-}^]^ being a family of periodic potentials with doubling periods and decreasing Loo-norms; 
namely, V T has orthogonal periods 2 r b%, 2 r £>2 and ||K-||oo < Cexp(—2 vr ) for some rj > 2 + 
64/ (2Z — 11). Without loss of generality we assume that (7 = 1 and Jq V r (x)dx = 0, Q r being 
the elementary cell of periods corresponding to V r (x). 

The one-dimensional analog of (jl.ljl . ()1.2() with I = 1 is already thoroughly investigated. It 
is proven in P~[Zj that the spectrum of the operator H\u = —u" + Vu is a Cantor type set. 
It has a positive Lebesgue measure The spectrum is absolutely continuous 0-0- 

Generalized eigenfunctions can be represented in the form of e lkx u(x), u(x) being limit-periodic 
The case of a complex- valued potential is studied in ^01- Integrated density of states is 
investigated in |llj-|14j. Properties of eigenfunctions of discrete multidimensional limit-periodic 
Schrodinger operators are studied in J^]. As to the continuum multidimensional case, it is proven 
|14j . that the integrated density of states for (jl.lj) is the limit of densities of states for periodic 
operators. Here we study properties of the spectrum and eigenfunctions of (jl.lj) . (jl.2j) in the high 
energy region. We prove the following results for the case d = 2, I > 6. 
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(1) The spectrum of the operator (jl.lj) . p.2|) contains a semiaxis. A proof of the analogous 
result by different means is to appear in the forthcoming paper In more general 
case, 8/ > d + 3, d ^ l(mod4), is considered, however, under additional restriction on the 
potential: the lattices of periods of all periodic potentials V r need to contain a nonzero 
vector 7 in common, i.e., V(x) is 7-periodic. 

(2) There are generalized eigenfunctions ^^(k, x), corresponding to the semiaxis, which are 
close to plane waves: for every k in an extensive subset Soo of M 2 , there is a solution 
fyoo^k, x) of the equation -ff^oo = ^oo^oo which can be described by the formula: 

*«,(£, x) = (l + Uoo (k, xj) , (1.3) 



hooh^) = Ofl^r 71 ), 71 >0, (1.4) 

fc|^oo ^ ' 

where Uoo(k,x) is a limit-periodic function: 

00 

Upotyix) = y^ y u r (k,x), (1.5) 

r=l 

u r (k,x) being periodic with periods 2 r ~ 1 bi, 2 r ~ 1 &2- The eigenvalue \co(k) corresponding 
to ^oo(k, x) is close to \k\ 21 : 

A M (fc) ,= \k\ 21 + O (\k\^A , 72 >0. (1.6) 

fc|— >oo ^ ' 

The "non-resonant" set Soo of vectors fe, for which ()1.3j) - (|1.6() hold, is an extensive Cantor 
type set: Soo = ^%L\Sn, where {Sn}^=i is a decreasing sequence of sets in M 2 . Each Sn 
has a finite number of holes in each bounded region. More and more holes appears when 
n increases, however holes added at each step are of smaller and smaller size. The set Soo 
satisfies the estimate: 

|SoonB R | = |B R |(i + o(ir^)), 73 >o, (1.7) 

fi^oo 

where Br is the disk of radius R centered at the origin, | • | is the Lebesgue measure in 
ro 2 



(3) The set D^A), defined as a level (isoenergetic) set for Xoo(k), 

©00(A) = {feGSoo: Aoo(fc) = A}, 



is proven to be a slightly distorted circle with infinite number of holes. It can be described 
by the formula: 



D 



00 



(A) = {£ : k = Xoo(A, v)u, v G 23oo(A)| , (1. 



where Boo (A) is a subset of the unit circle Si. The set Boo (A) can be interpreted as the 
set of possible directions of propagation for almost plane waves ([1.3)1 . The set (A) has 
a Cantor type structure and an asymptotically full measure on Si as A — > 00: 

L(03 oo (A)) = 27T + O (A-W2A (1.9) 
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here and below L(-) is a length of a curve. The value Xoo(A, u) in (|1.8|) is the "radius" of 
foo(A) in a direction V. The function xr 00 (A,i7) — A 1 / 2 ' describes the deviation of D 00 (X) 
from the perfect circle of the radius A 1 / 2 '. It is proven that the deviation is small: 

Xoo(A, v) = A 1 / 2 ' + O (A~ 74 ) , 74 >0. (1.10) 

We will prove absolute continuity of the branch of the spectrum (the semiaxis) corresponding 
to ^oo(k,x) in a forthcoming paper. 

To prove the results listed above we develop a modification of the Kolmogorov-Arnold-Moser 
(KAM) method. This paper is inspired by ^HEIEII; where the method is used for periodic 
problems. In the KAM method is applied to classical Hamiltonian systems. In |18[ I19j 
the technique developed in ^Tj\ is applied for semiclassical approximation for multidimensional 
periodic Schrodinger operators at high energies. 

We consider a sequence of operators 

M„ 

H = (-A) 1 , H {n) =H + ^2v r , n > 1, M n -> oo as n -> oo. 

r=l 

Obviously, \\H - # (n) || -> as n -» oo and ffW = + W„ where W n = E^=M n _ 1+ i *r- 

We treat each operator H^ n \ n > 1, as a perturbation of the previous operator H^ 11 ^. Every 
operator i?( n ) is periodic, however the periods go to infinity asn-> oo. We will show that there 
is a A*, A* = A*(V), such that the semiaxis [A*,oo) is contained in the spectra of all operators 
For every operator there is a set of eigenfunctions (corresponding to the semiaxis) 
close to plane waves: for any k in an extensive subset Sn of R 2 , there is a solution ^ n (k, x) of the 
differential equation H^ n ^ n = X^ n ^ n , which can be represented by the formula: 

^ n (k,x) = e 1 ^ (l + u n (k,x)), \\un\\ Loom = 0(\k\~^), 71 >0, (1.11) 

where u n (k,x) has periods 2 Mn ~ 1 b\,2 Mn ~ 1 b2- 1 The corresponding eigenvalue A^(fc) is close to 
|£| 2 <: 

AW(fc) = \kf l + O (\kr< 2 ) , 72 >0. 

|fc|^oo V ' 

The non-resonant set S n is proven to be extensive in M 2 : 

|3„nB R | = |B R |(l + 0(R-^)). (1.12) 

The set D n (A) is defined as the level (isoenergetic) set for non-resonant eigenvalue \( n \k): 

V n (X) = {ke9n- A< n >(fc) = A}. 

This set is proven to be a slightly distorted circle with a finite number of holes (see Fig. ^ 
the set £>i(A) being strictly inside the circle of the radius \ 1 / 21 for sufficiently large A. The set 
D n (A) can be described by the formula: 

D n (A) = [k : k = x n (A, z7 e £ n (A)} , (1.13) 
^Obviously, u n (k,x) is simply related to functions u r (k,x) used in 11.51 : u n (k,x) = X^r^M 1+1 u r(^> 
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Figure 1. Distorted circle with 
holes, Hi (A) 



Figure 2. Distorted circle with 
holes, T) 2 (A) 



where 23 n (A) is a subset of the unit circle Si. The set B n (A) can be interpreted as the set of 
possible directions of propagation for almost plane waves It has an asymptotically full 

measure on Si as A — » oo: 

L(S„(A)) = 2ir + o(\-^/ 21 ) . (1.14) 

The set 1> n (A) has only a finite number of holes, however their number is growing with n. More 
and more holes of a smaller and smaller size are added at each step. The value x n (A, v) — A 1 / 2 ' 
gives the deviation of D n (X) from the perfect circle of the radius X 1 / 21 in the direction v. It is 
proven that the deviation is asymptotically small: 

x n (A, Z 7)=A 1 / 2; + 0(A-^), dXn ^ V) =0{\-^) 74,75 > 0, (1.15) 

tp being an angle variable, v = (cos ip, simp). 

On each step more and more points are excluded from the non-resonant sets Sn, thus {Sn}$?Li 
is a decreasing sequence of sets. The set Soo is defined as the limit set: Soo = r\™ =1 $ n . It has an 
infinite number of holes, but nevertheless satisfies the relation (|1.7|) . For every k € Soo and every 
n, there is a generalized eigenfunction of of the type (|1.11|) . It is proven that the sequence 
of ^ n (k,x) has a limit in L DO (]R 2 ) when k £ Soo- The function ^^(kjx) = \im. n ^ o:> s S n {k,x) is 
a generalized eigenfunction of H. It can be written in the form (|1.3j) - (|1.5jl . Naturally, the 
corresponding eigenvalue Xooik) is the limit of X^ n \k) as n — > oo. 

It is shown that {"B n (X)}^ =1 is a decreasing sequence of sets, on each step more and more 
directions being excluded. We consider the limit Boo(A) of H n (X): 

oo 

Boo(A) = fl B„(A). (1.16) 

n=l 

This set has a Cantor type structure on the unit circle. It is proven that B^A) has an asymptoti- 
cally full measure on the unit circle (see (|1.9|l ). We prove that the sequence x n (A, i>), n = 1,2, 
describing the isoenergetic curves D n (A), quickly converges as n — > oo. Hence, Doo(A) can be 
described as the limit of D n (X) in the sense (|1.8[) . where Xoo(A, V) = lim^^oo x n (A, V) for every 
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V E Soo(A). It is shown that the derivatives of the functions x n (\,i>) (with respect to the angle 
variable on the unit circle) have a limit as n — ► oo for every V € < $> oa (X). We denote this limit by 



Thus, the limit curve D M (A) has a tangent vector in spite of its Cantor type structure, the tangent 
vector being the limit of corresponding tangent vectors for D n (\) as n — > oo. The curve 2) 00(A) 
looks as a slightly distorted circle with infinite number of holes. 

The main technical difficulty overcome is construction of non-resonance sets lB n (A) for every 
fixed sufficiently large A, A > Ao(V), where Ao(V) is the same for all n. The set !B n (A) is 
obtained by deleting a "resonant" part from 2$ n _i(A). Definition of S n _i(A) \ 23 n (A) includes 
Bloch eigenvalues of fl*" -1 ). To describe B n -i(A) \ B n (A) one has to use not only non-resonant 
eigenvalues of the type 1)1.6)1 . but also resonant eigenvalues, for which no suitable formulae are 
known. Absence of formulae causes difficulties in estimating the size of S n _i(A) \ H n (\). To deal 
with this problem we start with introducing an angle variable <p € [0, 2ir), v = (cos (p, sin (p) G 
Si and consider sets B n (A) in terms of this variable. Next, we show that the resonant set 
'Bn-i(A) \23 n (A) can be described as the set of zeros of determinants of the type Deb\I + S n (ip)) , 
S n ((p) being a trace type operator, 



where x n ~i((p) is a vector-function describing D n _i(A): i? n _i((/?) = x n -i(\,u)v. To obtain 
'Bn-i(A) \ S n (A) we take all values of e in a small interval and vectors b in a finite set, b ^ 0. 
Further, we extend our considerations to a complex neighborhood of [0, 2tt). We show that the 
determinants are analytic functions of ip in <Pq, and, by this, reduce the problem of estimating the 
size of the resonance set to a problem in complex analysis. We use theorems for analytic functions 
to count zeros of the determinants and to investigate how far zeros move when e changes. It 
enables us to estimate the size of the zero set of the determinants, and, hence, the size of the non- 
resonance set <P n C which is defined as a non-zero set for the determinants. Proving that the 
non-resonance set is sufficiently large, we obtain estimates 1)1.12)1 for S n and ()1.14)) for 23 n , the 
set "Bn being the real part of <& n . To obtain $> n we delete from ^0 more and more holes of smaller 
and smaller radii at each step. Thus, the non-resonance set # n C <?o has a structure of Swiss 
Cheese (Fig. |S1 pages QUI US)) . Deleting resonance set from at each step of the recurrent 
procedure we call a "Swiss Cheese Method". The essential difference of our method from those 
applied in similar situations before (see e.g. |17j-|2U|) is that we construct a non-resonance set 
not only in the whole space of a parameter (k G 1R 2 here), but also on isoenergetic curves D n (X) 
in the space of the parameter, when A is sufficiently large. Estimates for the size of non-resonance 
sets on a curve require more subtle technical considerations than those sufficient for description 
of a non-resonant set in the whole space of the parameter. 

The restriction I > 6 is technical, it is needed only for the first two steps of the recurrent 
procedure. The requirement for super exponential decay of ||Vr|| as r — ► 00 is more essential, 
since it is needed to ensure convergence of the recurrent procedure. It is not essential that the 
potential V r has doubling periods, the periods of the type q r ~ 1 bi, q r ~ 1 b2,Q £ N, can be treated 
in the same way. First consideration of the problem was given in |23j . 



<9xoo(A,i7) 
dip 



. It follows from 1)1.15)1 that 




(1.17) 



I + S n (tp) = {H( n - l \x n _i{ip) + $) -A-e) (iJ (x n _iM + &)+A) 
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The plan of the paper is the following. Sections 2-6 describe steps of the recurrent procedure. 
Discussion of convergence of the procedure and proofs of the results 1-3, listed above, are in 
Section 7. 

Acknowledgement The authors are very grateful to Prof. G. Stolz and Prof. G. Gallavotti for 
useful discussions. 

2. The First Approximation 

2.1. Operator H^ l \ We introduce the first operator H^ l \ which corresponds to a partial sum 
in the series 1)1 .2(1 : 

Mi 

#(i) = (_ A )' + Wi, W^^Vr, (2-1) 

r=l 

where Mi is chosen in such a way that 2 Ml ns k Sl 2 for a jfc > 1, sx = (21 - ll)/32. Obviously, 
the periods of W\ are (a t ,0) = 2 Ml ~ 1 (&i,0) and (0,a 2 ) = 2 Ml ~ 1 (0,6 2 ), and ai « k Sl b 1 /2, a 2 « 
/c Sl 6 2 /2. Note that ||Wi||oo < E^i l|K||oo = 0(1) as A; -> oo. 

It is well-known (see e.g. that spectral analysis of a periodic operator can be reduced 
to analysis of a family of operators H^\t), t £ .Ki, where K\ is the elementary cell of the dual 
lattice, K\ = [0, 27ra 1 ~ 1 ) x [0, 27ra 2 ~ 1 ). The vector t is called quasimomentum. An operator 
H^'it), t £ A'i, acts in ^2(^1)5 Qi being the elementary cell of the periods of the potential, 
Qi = [0, a{\ x [0,02]. The operator H^\t) is described by formula (|2.1|) and the quasiperiodic 
conditions for a function itself and its derivatives: 

u(ai,x 2 ) = exp(^iai)u(0,x 2 ), u(xi, a 2 ) = exp(zt 2 a 2 )u(xi, 0), (2.2) 

it^(ai,x 2 ) = exp(itiai)n^(0,x 2 ), w^(xi,a 2 ) = ex.p(it 2 a 2 )u^ (xi, 0), 

< j < 2Z. Each operator H^(t), t E Ki, has a discrete bounded below spectrum A^(t): 

A«(t) = U~ ^(t), A«(i) -w, 00. 

The spectrum A^ 1 ) of the operator is the union of the spectra of the operators H^(t): 

AW = U t6 KiA(i) = U nEN)tEK An\t). The functions A^(t) are continuous, so A^ has a band 
structure: 

A« =U- 1 [g«,QW], ? W = minA«(i), Q« = max A«(t). (2.3) 

Eigenfunctions of H^(t) and ii^ 1 ) are simply related. Extending all the eigenfunctions of the 
operators H^'(t) quasiperiodically (see (|2.2|) ) to M 2 , we obtain a complete system of generalized 
eigenfunctions of .H"' 1 ). 

Let Hq 1 ^ be the operator Ijl.lj) corresponding to V = 0. We consider that it has periods 
a\,a 2 and that operators Hq (t), t € K\, are defined in L 2 (Q\). Eigenfunctions of an operator 
Hq (t), t G Ki, are plane waves satisfying (|2.2|) , They are naturally indexed by points of Z 2 : 
*!/!■(£, a;) = |Qi| _1//2 expi{pj(t),x), j £ Z 2 , the eigenvalue corresponding to ^(t,x) being equal to 



2 We write a(k) « b(k) when the inequalities \b{k) < a(k) < 2b(k) hold. 
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pf(t), here and below Pj(t) = 2-irj/a + t, 2-Kj/a = (2irjx/ai, 2^2/^2), j € Z 2 , |Qi| = aia2 and 

P f(t) = \m\ 21 - 

Next, we introduce an isoenergetic surface 3 Sb(A) of the free operator Hq . A point t belongs 
to S'o(A) if and only if H^ft) has an eigenvalue equal to A, i.e., there exists a j € 1? such that 
p 2i (i) = A. This surface can be obtained as follows: the circle of radius k = A21 centered at the 
origin is divided into pieces by the dual lattice {Pq(0)} q ez 2 ; an d then all pieces are translated in a 
parallel manner into the cell K\ of the dual lattice. We also can get So (A) by drawing sufficiently 
many circles of radii k centered at the dual lattice {Pq(Q)}q£j? an d by looking at the figure in 
the cell K\. As a result of any of these two procedures we obtain a circle of radius k "packed 
into the bag Ki" as it is shown in the Fig. El Note that each piece of So (A) can be described by 
an equation pf(t) = A for a fixed j. If t € So (A), then j can be uniquely defined from the last 
equation, unless t is not the point of a self-intersection of the isoenergetic surface. A point t is a 
self-intersection of So(A) if and only if 

pf(t)=pf(t)=k 21 (2.4) 

for at least on pair of indices q, j, q 7^ j. 

Note that any vector x in M. 2 can be uniquely represented in the form x = Pj(t), where j € Z 2 
and t G K\. Let %i be the parallel shift into Kf. X 1 : R 2 -» K x , %ipj(t) = t. Suppose SlcR 2 . 
In order to obtain DCifi it is necessary to partition f2 by the lattice with nodes at the points p q (0), 
g£Z 2 and to shift all parts in a parallel manner into a single cell. It is obvious that |3Ci$l| < |0| 
for any $7. If f2 is a smooth curve, then 

L(X 1 Q) < L(fi). (2.5) 



"surface" is a traditional term. In our case, it is a curve. 



s 
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Figure 4. The first non-resonance set Xi(X,S) 

For any pair of sets J7i and Q2, 3Ci(f2i U O2) = U XiQ.2- Obviously, 3C\Sk = So (A) an d 

L(S (X)) = L(S fc ) = 2nk, k = \2i,S k being the circle of radius k centered at the origin. 

The operator has the following matrix representation in the basis of plane waves: (t) mq 
Pm{t)5 mq + w m - q , 4 here and below w m - q are Fourier coefficients of W\, the coefficient wq being 
equal to zero. 

2.2. Perturbation Formulae. In this section we consider operator H^\t) as a perturbation of 

the free operator H^'(t). We show that for every sufficiently large A there is a "non-resonant" 
subset Xl(A) of So (A) such that perturbation series for an eigenvalue and a spectral projection of 
converge when t E Xi(A)- The set xi(A) is obtained by deleting small neighborhoods of 
selfintersections of So (A), see Fig. |1J The selfintersections are described by (|2.4|) and correspond 

to degenerated eigenvalues of (t). The size of the neighborhood is k~ 1 ~ 4si ~ s , k = X^i f 5 being 
a small positive number. The set Xi(A) is sufficiently large: its relative measure with respect to 
So (A) tends to 1 as A — > 00. The precise formulation of these results is given in the next Geometric 
Lemma. The lemma is proven by elementary geometric considerations in |22j . 

Lemma 2.1 (Geometric Lemma). For an arbitrarily small positive 5, 25 < 21 — 2 — 4si, and 
sufficiently large A, A > Ao(<5), there exists a nonresonance set xi(A, 8) C So(A) such that: 

(1) For any point t G Xl the following conditions hold: 

(a) There exists a unique j G such that Pj(t) = k, k = X^i . 

(b) The inequality holds: 

mm\p 2 j (t)-p 2 i (t)\>2k-^- S . (2.6) 

(2) For any t in the (2/c _1_4si_25 ) -neighborhood of the nonresonance set in C 2 , there exists a 
unique j € I? such that 

\p 2 {t) - k 2 \ < 5k- 4si ~ 2S (2.7) 
and the inequality \2. b}) holds. 



5m q is the Kronecker symbol. 
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(3) The nonresonance set xi has an asymptotically full measure on Sq(X) in the following 
sense: 

Corollary 2.2. If t belongs to the (2fc _1_4si_2<5 ) -neighborhood of the nonresonance set xi(A,<5) 
in C 2 , then, for any z lying on the circle 

d = {z:\z- k 2l \ = k 21 ' 2 -^- 5 } (2.9) 

and any i in J?, the following inequality holds: 

2\p 2l (t)-z\>k 2l - 2 - 4si - s . (2.10) 

Corollaru \2. t A Let t £ xi(A,<5). Taking into account the relation pf{t) = k 21 and the definition 
of Ci, we see that 

\pf(t) - z \ = k 2l ~ 2 ~^- 5 . (2.11) 
Thus, the estimate (j2.1Uj) is valid for i = j. From (|2.6|) . it follows that min^j \p 2l (t) — pf{t)\ > 

-2Z-2-4si-<5 



2k Sl ~ . Using the last inequality and (J27TTJ) , we obtain inequality (|2.10j) for i ^ j. It is easy 
to see that all the estimates are stable under a perturbation of t of order 2/c _1_4si_2<5 . Therefore, 
the estimate ()2.1Uj) can be extended to the (2/c _1_4si_2<5 )-neighborhood of xi- I 

Let Ej(t) be the spectral projection of the free operator, corresponding to the eigenvalue pf{t) : 
(Ej) rm = 5j r 5j m . In the (2/c _1_4si_2<5 )-neighborhood of 5), we define functions gr(k, t) and 
operator-valued functions Gj^\k, t), r = 1,2, • • • as follows: 

gW(k,t) = ^Tr t - z^WxYdz, (2.12) 

J G~\_ 

G?Hk,t) = I {{H W {t) _ z )^ Wl y(H^\t) - z)' l dz. (2.13) 

J C\ 

To find gi 1] (k,t) and G^\k,t) it is necessary to compute the residues of a rational function of 
a simple structure, whose numerator does not depend on z, while the denominator is a product 
of factors of the type ipfit) — z). For all t in the nonresonance set within C\ the integrand has 
a single pole at the point z = k 21 = pf{t). By computing the residue at this point, we obtain 

explicit expressions for gi 1 \k,t) and Gy For example, g^\k,t) = 0, 

g£\k,t)= £ \w q \ 2 (pf(t)-p 2l +q (t))- 1 

v k g | 2 (2pf(t)-^ g (t)- P 2 L 9 (t)) 

G[ \k,t) rm = 2in\ j W77\^ + 21 /4\ J \u+\^ m ii G[\k,t)jj = 0. (2.15) 

It is not difficult to show that g^\k,t) > for sufficiently large A. For technical reasons it is 
convenient to introduce parameter a in front of the potential W\. Namely, H$ = (-A) l + aWi, 
< a < 1, h{ 1] = HW. 
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Theorem 2.3. Suppose t belongs to the (2k~ 1 ~ Asi ~ 2& )-neighborhood in K\ of the nonresonance 
set xi(A, 6), < 25 < 21 — 2 — As\. Then for sufficiently large A, A > Ao(||V||, bi, 62, 5), for all 
a, — 1 < a < 1, there exists a single eigenvalue of the operator Ha\t) in the interval ei(k,S) := 

(k 21 — ^2/— 2— 4si— 5^ ^21 _|_ £,2/-2-4si-<5)_ j£ ^ g^ ven fry fo e ser i es: 

00 

Af (a, t) = pf(t) + 2 « r ^ 1} (fc, <), (2-16) 

r=2 

converging absolutely in the disk \a\ < 1, where the index j is determined according to articles 
1(a) and 2 of Geometric Lemma. The spectral projection, corresponding to Xj (a,i), is given by 
the series: 

00 

Ef\a,t) = Ej + Y,a r Gi 1] (k,t), (2.17) 
r=X 

which converges in the class Si uniformly with respect to a in the disk \a\ < 1. 
For coefficients gi 1 \k,t), Gr (k,t) the following estimates hold: 

|<7«(M)| < ^-2-4 S1 - 70 r- 5) ||G«(M)||l < k^ or , (2.18) 

where 70 = 21 — 2 — 4s± — 25. 

Corollary 2.4. For the perturbed eigenvalue and its spectral projection the following estimates 
are valid: 

\xf\a,t)-pf(t)\ < 2a 2 k 2l ~ 2 - 4s ^°~ 5 , (2.19) 

\\Ef ] (a,t) - EjWx < 2\a\k-^°. (2.20) 

Remark 2.5. The theorem states that X^\a,t) is a single eigenvalue in the interval £\(k, 5). 
This means that \xf\a,t) - k 2l \ < fc2Z-2-4 Sl -5_ Formula (j2~19l) provid es a stronger estimate on 
the location of Xj\a,t), the right-hand side of (|2.19|) being smaller than the size of e±. 

Proof. The proof of the theorem is based on expanding the resolvent (Ha\t) — in a per- 
turbation series for z, belonging to the contour C\ about the unperturbed eigenvalue p 2l (t). 
Then integrating the resolvent yields the formulae for the perturbed eigenvalue and its spectral 
projection. 

It is obvious that 

- z)- 1 = (2f«(t) - z)~\l - aA 1 y 1 , A, := -W,{H^\t) - z)' 1 . (2.21) 

Suppose z G C\. Using (|2.1UI) . we obtain: 

\\(H^(t) - z)- 1 !) < 2k- 2l+2+4si+s , < 2\\W 1 \\k-^ - 5 . (2.22) 

Thus, ||^4i|| <C 1 for sufficiently large k. Expanding (/ — aAi) -1 in a series in powers of a.A\ we 
obtain: 

00 

(Hj»(t) - z)- 1 - (H^(t) - z)- 1 = J2 oT(H^(t) - zy'Al. (2.23) 

r=l 
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Since 21 — 2 — 2s\ > it is not hard to show that (H^ (t) — z)^ 1 6 Si. Taking into account 
estimates (|2.22|) . we see that the series ()2.23(l converges in the class Si uniformly with respect 
to a in the disk \a\ < 1. Substituting the series into the formula for a spectral projection 

E^p {a,t) = — §q (Ha\t) — z)~ 1 dz and integrating termwise, we arrive at ()2.17|) . Formula 
()2.16j) is obtained in the analogous way. Estimates (|2.18|) follow from (|2.22|) . for details see |22j . 

■ 

Next, we show that the series (j2.16|) . (j2.17|) can be extended as holomorphic functions of t in 
a complex neighborhood of Xl'i they can be differentiated with respect to t any number of times 
with retaining their asymptotic character. 

Let us introduce the notations: 

Q\m\ 

T(m) := g f mig t m 2 > m=(mi,m 2 ), \m\ := m% + m 2 , ml := mi!m 2 !, T(0)f := f. 

Lemma 2.6. The coefficients g? '(A;,i), and (k,t) can be continued as holomorphic functions 
of two variables from the real (Ik -1-481-25 ) -neighborhood of the nonresonance setxi to its complex 
(2/c~ 1 ~ 4si ~ 2<5 ) -neighborhood and the following estimates hold in the complex neighborhood: 

\T(m)gW(k t t)\ < m !fc2/-2-4 Sl - < 5- T „r+|m|(l+4s 1 +2 ( 5) ? 

||r(m)(?W(fc,i)||i < m \k-^ r+ ^ (1+ ^ +25 \ (2.25) 

Proof. Since (|2.1()|) is valid in the complex (2fc _1_4si_2(5 )-neighborhood of the nonresonance set, 
it is not hard to see that the coefficients g^{k,t), and G^\k,t) can be continued from the 
real (2/c _1 ~ 4si_25 )-neighborhood of t to the complex (2A;~ 1_4si ~ 2<5 )-neighborhood as holomorphic 
functions of two variables, and inequalities (|2.18|) are hereby preserved. Estimating by means of 
the Cauchy integral formula, the value of the derivative with respect to t in terms of the value of 
the function itself on the boundary of the (2/c _1_4si_2,5 )-neighborhood of t (formulas (|2.18|) ). we 
obtain (t2~2l1) and (JT25J). ■ 

From this lemma the following theorem easily follows. 

Theorem 2.7. The series h2.1b]) . \2.17\ ) can be continued as holomorphic functions of two 
variables from the real (2k~ 1 ~ Asi ~ 2S )-neighborhood of the nonresonance set \i to its complex 
(2A; _1_4si_2<5 ) -neighborhood and the following estimates hold in the complex neighborhood: 

\T(m)(\f\a,t) -pf(t))\ < 2 m \a 2 k 2l - 2 -^- 5 - 2 ^ + ^ (1+ ^ +25 \ (2.26) 



\T{m){Ef\a,t) - E^lll < 2m!a£;- 70+H(1 + 4si + 25) . (2.27) 



Corollary 2.8. 



V\f\a,t) = 2lk 2l - 2 k + o [k 21 - 1 ^ , k = pj(t), (2.28) 

T(m)Xf\a,t) < Al 2 k 2l ~ 2 , if \m\ = 2. (2.29) 

The next lemma will be used in the second step of approximation where the operator H^'(t) 
will play a role of the initial (unperturbed) operator. 
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Lemma 2.9. For any z on the circle C\ given by \2. <J\) and t in the complex (2fc~ 1_4si_2<5 )- 
neighborhood of \i, 

|| (ff « (t) - zY 1 1| < Ak~ 2l+2+4si+s . (2.30) 
Proof. Considering the following Hilbert relation 

(2T«(t) - zT 1 = (fT«(t) - z)' 1 + (2f«(t) - z)-\-W x ){H^\t) - z)' 1 
and (t2~TuT) . i.e., IK-H^t) - 1 1| < 2£r 2/+2+4si+<5 , we obtain 

\\(HW(t) - z)- 1 ]] < — IK^ W-^r 1 !! < 2 \\(H^(t) - z)- 1 )! < 4 A- 2 '+ 2 + 4 ^. 

i - ik^-w (*) - 

■ 

2.3. Nonresonant part of the Isoenergetic Set of i/^ 1 ). Let Si (A) 5 be the isoenergetic 
surface of the operator H a 

Si (A) = {t € i^i : 3n G N s.t. A^(a,i) = A}, (2.31) 

where {An (ck, i)}^=i is the complete set of eigenvalues of H a l \t). Next, we construct a "nonres- 
onance" subset x*(A) of Si (A), which corresponds to nonresonance eigenvalues A^ (a,i) given by 
the perturbation series. 

Let us note that for every t belonging to the non-resonant set Xi(A, $) described by Lemma l2.11 
there is a single j G Z 2 such that Pj(t) = k, k = X^l. This means that the function t — ► Pj(t) maps 
Xi(A, S) into the circle S&. The image of Xi(A, S) in S& we denote by r Do(X) nonres . Obviously, 

Xi(X,S) =X 1 V (X) nonres , (2.32) 

where %\ establishes one-to-one relation between two sets. Let 23i (A) is a set of unit vectors 
corresponding to I'o(A) rionr . es : 

Bi(A) = {V G Si : kv G D (A)„onr e ,}. 

It is easy to see that 'Bi(A) is a unit circle with holes, centered at the origin. We denote by 0i(A) 
the set of angles ip in polar coordinates, corresponding to Si(A): 

®i(A) = {^€ [0, 2tt) : (cos sin 99) G Bi(A)}. (2.33) 

Let x G 1>o(k) nonres . Then there is j G Z 2 , t G Xi(A, <5) such that x = Pj(t). 6 Obviously, 
t = %\x and, by (|2.32|) . t G Xi(A, S). According to Theorem 12.31 for sufficiently large k, there 
exists an eigenvalue of the operator H a (t), t = < a < 1, given by (|2.16|) . It is convenient 
here to denote Xy(a,t) by X^\a, if), we can do this since there is one-to-one correspondence 
between i? and the pair (t,j). We rewrite ()2.16|1 in the form: 

oo 

A«(a,x) = x 2i + / 1 (a,x), x=|x|, /i(a,if) = J> r 5 «(x), (2.34) 

r=2 



^5i(A) definitely depends on aW\, however we omit this to keep the notation simple. 

^Usually the vector Pj(t) is denoted by k, the corresponding plane wave being e^*'*^. We use less common 
notation k, since we have already other fe's in the text. 
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<7r (ir) being denned by (|2.12|) with j and t such that Pj(t) = x, and fi(a,x) satisfying the 
estimates: 

|/i(a,x)| < 2a 2 k 2l - 2 - 4si -^°' 5 , (2.35) 

|T(m)/i(a,x)| < 2a 2 k 2l - 2 - Asi -^°- 5+ ^ l+Asi+25 \ (2.36) 

By Theoreml2~3l the formulae (t2~34l - (t2~35l hold in the (2Ar 1_4si ~ 25 )-neighborhood of T> (\) nonre& 
i.e., they hold for any xv such that v £ 'Bi(A), \x — k\ < k~ l ~ 4si ~ 2& . We define Di(A) as the 
level set of the function A^ 1 ) (a, x) in this neighborhood: 

Dl(A) := {it = itv : V € Si(A), \x - k\ < Ar 1 - 451 " 25 , A (1) (a, it) = A}. (2.37) 

We prove in Lemma 12.111 that T>i(X) is a distorted circle with holes, which is close to and inside 
of the circle of radius k, Fig. ^ First, we prove that the equation \^\a, x) = A is solvable with 
respect to x = \x\ for any v = ~ € Bi(A). 

Lemma 2.10. For every v G 5H(A) and every a, < a < 1 and sufficiently large \, there is a 
unique x\ = Xi(A,i7) in the interval I\ := [Ac — A;~ 1_4si ~ 2<5 , + /c _1 ~ 4si_2<5 ], A; 2 ' = A, sitc/i i/iai 

A^(a, xiu) = A. (2.38) 
Furthermore, \x\ — k\ < c ^- 1 - 4s i- 2 7o-<5 > 

Proof. 1) Let us show that there exists such an x\. Given v £ Bi(A), by Theorem 12.31 there 
exists an eigenvalue A^ (a, xv) given by (|2.34l) for all x in the interval I\ with a sufficiently large 
k. Obviously, \W(Q,kv) = pf(t ) = k 21 = A, \W(0,xv) = pf(t) = x 21 , pj(t ) = kv, pj(t) = xv, 
where j, t, and to are uniquely determined by xv and kv. Let H^\a, v) := {A = A^ 1 ) (a, xv) : 
x £ Using the definition of I\ and considering that A^^(0,xz?) = x 21 , we easily obtain 
£W(0,i?) D [k 2l -5i,k 2l +5i], 5i = Cl k 21 - 2 - 4 ^- 25 , < a + ci(fc). Since \^{a,xv) is continuous 
in x and (|23B1 is valid, we have £W(a, v) D [£ 2Z - <5i + S 2 , k 21 +S X - S 2 ], S 2 = C2 k 2l - 2 ' As ^- 2 ^- 5 . 
Since 2 7o - 5 > 0, we have £W(a, i?) D [A; 2i - 5i/2, A; 2 ' + £ x /2]. Thus for every v € 'Bi(A), there 
exists a xi such that \^ l \a, x\v) = A; 2 ', x\ E I\. 

2) Now we show that there is only one x\ satisfying (|2.38() in the interval 1%. Differentiat- 
ing (j2.34j) with respect to x and using (|2,36[) . we get: 

aA(1) ^' Xi7) = 2/x 2 '- 1 + Oik 21 - 1 - 2 ^ 5 ). (2.39) 
ox 

Therefore, 

dX^(a,xv) 



d 



x 



> 21k 21 - 1 (1 + o{l)) , (2.40) 

and this implies that \^ {a, xv) is monotone with respect to x in I\. Thus, there is a single 
x\ G Ji satisfying (|2.38|) . 

3) Now let us estimate \x\ — k\. By (|2.34j) and (j2,35|) we have 

A«(a, kv) = k 21 + /i(a, fcz?) = A; 2 ' + ( fc 2Z-2-4 Sl -2 7o -^ ( 2 41 ) 

Using (H3U, fl23EJ), l|Q0|) . and the Mean Value Theorem, we obtain 

cA .2/-2~4 Sl -2 70 ^ > ^(a,*!*?) - A« («, kv) \ > 
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mm 



xv) 



and hence we arrive at \x\ 



dx 



x x -k\> ck 2l ^\x x - k\ 



Here and below, we use the notations xi(cp) := xi(A, V), v = (cos ip, sin ip), x.\{ip) = x\{ip)v, h\{ip) 
x x (ip) - k. 



Lemma 2.11. (1) The set X>i(A) is a distorted circle with holes which is strictly inside the 
circle of the radius k: it can be described by the formula: 

T>i(A) = {x G M 2 : x = xi(<p)v, ^G'Bi(A)}, 

k + h\{ip) and h\{ip) obeys the inequalities: 

dhi 



(2.42) 



where x\{ip) 



\hi\ < ck' 



-1— 4si— 270— 8 



dip 



< ck 



-270+1 



(2.43) 



hi{ip) < when W\ 7^ 0. 

(2) The total length of ¥>i(\) satisfies the estimate: 

= 2vr(l + 0{k- 5/2 )). (2.44) 

(3) Function h\(ip) can be extended as a holomorphic function of ip to the complex (2fc _2 ~ 4si_25 ) 
neighborhood of each connected component of 0i(A) and estimates \2.43\) hold. 

(4) The curve T>i (A) has a length which is asymptotically close to that of the whole circle in 
the following sense: 



Proof. 



(1) Let x G Di(A), then x = 



\=k 

x\V. Using Lemma l2.1()| we get |/ii(<^)| 



L(Di(A)) = 2irk(l + 0(k- 5 ' 2 )), A 

A— >oo 



cA ,_l_4 Sl -2 7 o-<5 and 



dip 



dtp 

A = xf + fi(xi cos ip, xi sin ip) 
and, therefore, by implicit differentiation, 



(2.45) 
\x\ — k\ < 

From (|2.34j) and the definition of Di(A), we have 

(2.46) 



dhi 




dx\ 




dip 




dip 





dfi 

dip 



2lxf- 1 + (V/,z7) 



(2.47) 



Using (|2.36|) . we easily obtain |V/i| < ck 21 270 l , and, hence, ()2.43|) for the derivative. 
It is easy to show that g^\k,t) given by (|2. 14f) is positive when W\(x) is non-zero, and, 
moreover, fi(a, x) > 0. It easily follows that h\ < 0. 

(2) By definition, , B 1 (A) is the set of directions corresponding to r Do(X) n onresi the latter set 
being a subset of the sphere of radius k. Formula (|2.32j) establishes one-to-one correspon- 
dence between Xi(^,$) an d "DoWnonres, their length being equal. Considering (|2.8|) . we 
obtain L {V (X) nonres ) = 2vrA:(l + 0(k" 5 / 2 )). Hence, (l2~Hl) holds. 

(3) This part of the lemma easily follows from the Implicit Function Theorem applied to 
(El 



the above estimate for |V/i| and the fact that f\(a,x) is a holomorphic function 
of two variables (Theorem 12. 7|) . 
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Figure 5. The set Xi(A,<5) 



(4) Considering that L(Di(A)) = Jg 1 m v (<£>) + ^Kv 9 ) 2 ^V 9 an d taking into account (|2.44|) 
and (|2.47|) . we obtain 



L(D (A))-L(!Di(A)) = / fcdp+/ ffc - Jxf^) + x^) 2 ) dip 

-'[o,27t]\0i(a) Je~ v v 7 



< / / (/iiM + ^iM) 

J[0,27r]\ei(A) J0iW 

< ck ■ k- 5 / 2 + O(jt-l-4 S1 -270-5) + O ( fc -270+1) = f, . 0(k~ 5 / 2 ), 

here and below, Do (A), corresponding to the free operator, is the perfect circle of radius 
A 1 / 2 '. Hence, (Z3g) is valid. 

■ 

Next, we define a nonresonance subset x*(^) °f isoenergetic set Si (A) as the parallel shift of 
D X (A) into Ki (Fig. EJ): 

X*(A) := 3CiDi(A). (2.48) 

Lemma 2.12. The set x*(A) &ekm#,s to £/ie (ck^ 1 - 431 - 2 ^ " 8 ) -neighborhood o/xi(A) m ifi. // 

i G Xi(^); i^-en £/ie operator (t) has a simple eigenvalue Xn\a,t), n G N ; egita/ to A ; no 
oi/ier eigenvalues being in the interval e\(k,S), £\{k,5) := (k 21 — k 2l ~~ 2 ~ 4si ~~ s ,k 21 + fc 2i ~ 2 ~ 4si_<5 ). 
This eigenvalue is given by the perturbation series h2.1b}) . j being uniquely defined by t from the 
relation pf{t) G E\(k,5). 

Proof. By Lemma l2~TTl T>i(A) is in the (cfe _1 "~ 4si_2 ' ft_5 )-neighborhood of D (A). Consider- 
ing that XiW = OCiDo(A) and XiM = ^lfi(A), we immediately obtain that XiM is m the 
(cfc~ 1 ~ 4si ~ 270 ~ <5 )-neighborhood of XiM- The size of this neighborhood is less than /c~ 1_4si-2<5 , 
hence Theorem 12 .31 holds in it, i.e. for any t G X\W there is a single eigenvalue of (t) in the in- 
terval £\(k, 5). Since XiW C 5*i(A), this eigenvalue is equal to A. By the theorem, the eigenvalue 
is given by the series (j2.16|) . j being uniquely defined by t from the relation pf{t) G E\(k, 5). | 



Lemma 2.13. Formula \2.4$j) establishes one-to-one correspondence between x*W and T>i(\). 
Remark 2.14. From geometric point of view this means that XiW does not have self-intersections. 
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Proof. Suppose there is a pair if 1,1, i?i,2 G I*i(A), such that %\5k\^\ = Xiifi^ = t, t £ X*(A). By 
the definition (|2,37|) of Di(A), we have A^ (a, ifi,! = A, i.e., the eigenvalue A of 

Ha\t) is not simple. This contradicts to the previous lemma. ■ 

3. The Second Approximation 

(2) . (2) 

3.1. The Operator Ha . Choosing s 2 = 2s\, we define the second operator H a by the formula: 

H (2) = H {l) + aW ^ (0 < a < 1), W 2 = Yl V r> i 3 - 1 ) 

r=Mi+l 

where is defined by ()2.1j) . M 2 is chosen in such a way that 2 M2 ~ k S2 . Obviously, the periods 
of W 2 are 2 M2 ~ 1 (b 1 , 0) and 2 M2 ~ 1 (0, b 2 ). We will write them in the form: iVi(ai, 0) and iV"i(0, a 2 ), 
where ai,a 2 are the periods of W\ and JVi = 2 M2 ~ Ml , l^ 2 " 51 < TVx < ik S2 ~ Sl . Note that 

||W 2 ]|oo < J2 W V r\\oo< Yl exp(-2^) <exp(-F Sl ). (3.2) 

r=M 1 + l r=Mi+l 

3.1.1. Multiple Periods ofWi(x). The operator = Hq + Wi(x) has the periods a\,a 2 . The 
corresponding family of operators, {H^ {t)}teKx , acts in L 2 (Q\), where Q\ = [0,a±] x [0, a 2 ] 
and Ki = [0, 2ir/a\) x [0, 2ir/a 2 ). Eigenvalues of H^(t) are denoted by An (t), n G N, and its 
spectrum by A^>(t). Now let us consider the same W\(x) as a periodic function with the periods 
Niai, Nia 2 . Obviously, the definition of the operator does not depend on the way how we 
define the periods of W\. However, the family of operators {H^\t)}t£Ki does change, when we 
replace the periods a±, a 2 by N\a±, N\a 2 . The family of operators {H^ (t)} tg ^ 1 has to be replaced 
by a family of operators {H^ 1 \t)} t& k 2 acting in L 2 (Q 2 ), where Q 2 = [0, Nim] x [0, iVia2] and 
K 2 = [0,2vr/iViai) x [0, 2vr/iVia 2 ). We denote eigenvalues of H^\t) by >${j), n G N and its 
spectrum by A.W(r). The next lemma establishes a connection between spectra of operators 
H^\t) and H^'(t). It easily follows from Bloch theory (see e.g. |21j). 

Lemma 3.1. For any r G K 2 , 

A«(r) = (J A<%), (3.3) 

where 

P = {p= ( Pl ,p 2 ) G 1? : < pi < iVi - 1, < p 2 < JVi - 1} (3.4) 
and t p = (t p ,i,tp, 2 ) = (n + 2npi/Niai,T2 + 2irp 2 /Nia 2 ) G iTi, see Fi#. £. 

We defined isoenergetic set «Si(A) C K\ of i^ 1 ) by formula (|2,31|) , Obviously, this definition 
is directly associated with the family of operators H^ l \t) and, therefore, with periods ai,a 2 , 
which we assigned to W\{x). Now, assuming that the periods are equal to N\a\, N±a 2 , we give an 

analogous definition of the isoenerg etic set Si (A) in K 2 : 5i(A) := {r G K 2 : 3n G N : )&\t) = 
A}. By Lemma |3~T1 Si (A) can be expressed as follows: 

5i(A) = {r G K 2 : 3n G N, p G P : A^ (r + 2vrp/iVia) = a}, 
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Figure 6. Relation between r (x) and L 



2np/N 1 

a — ( ATi ai ' jVi 02 ) ' ^^e relation between Si(A) and Si(A) can be easily understood from 
the geometric point of view as Si (A) = X2Si(\), where %2 is the parallel shift into K2, i.e., 
%2 ■ -» #2, K 2 (r + 2vrm/iVia) = r, m G Z 2 , r £ K 2 . Thus, S X (A) is obtained from Si (A) by 
cutting Si (A) into pieces of the size K2 and shifting them together in K2. 

Definition 3.2. We say that r is a point of self-intersection of Si(A), if there is a pair n,n G N, 
n^h such that A„ 1} (r) = \%\t) = A. 

Remark 3.3. By Lemma 13. 11 r is a point of self-intersection of Si (A), if there is a pair p,p G P 
and a pair n,n G N such that \p— pj + jn — ra| 7^ and A n (T + 27rp/iVia) = X^ 1 {r + 2Ttp / N\a) = A. 

Now let us recall that the isoenergetic set Si (A) consists of two parts: Si (A) = x*(A) U 
(Si(A) \Xi(A))> where Xi(A) is the first non-resonance set given by ()2.48j) . Obviously X2Xi(ty c 
3C2Si(A) = Si (A) and can be described by the formula: 

K 2 xt(A) = {r G K 2 : 3p G P : r + 2^p/iVia G Xi(A)} . 

Let us consider only those self-intersections of Si which belong to 3<C2Xi(A), i.e., we consider the 
points of intersection of X2X1W both with itself and with Si (A) \ %2Xi(ty- 

Lemma 3.4. A self-intersection r of Si (A) belongs to 3C 2 x*(A) if and only if there are a pair 
p,p G P, p 7^ p and a pair n,h G N such that r + 2irp/Nia G X*(A) an d An ( T + 2-Kp/N\a) = 
A^(r + 27rp/iVia) = A ; the eigenvalue A^ (r + 2np/Ni a) being given by the series \2. 1 b]) with t = 
T + 2irp/Nia and j uniquely defined byt from the relation p 21 (t) G S\, E\ = (k 21 — k 2l ~ 2 ~ 4si ~ s , k 21 + 

Proof. Suppose r is a point of self-intersection of Si (A) belonging to 3C 2 Xi(A)- Since, r G 3<C2Xi(A), 
there is a p G P such that r + 2np/Nia G X*(A). By Lemma 12.121 there is an eigenvalue 
An^(r + 2-Kp/N\a) = A, which is given by the series (|2.16|) with t = r + 2irp/Nia and j uniquely 
defined by i from the relation pf{i) G e%. By the same lemma \u\t + 2irp/N\a) is a unique 
eigenvalue of (r + 2-Kp/Nia) in ej_. This means 



AW(T + 2vrp/Afia) - a1 1} (t + 2TT P /N ia ) > ^-2-4^-^ when n/n. 



(3.5) 
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Since r is a point of self-intersection of Si (A), An (r + 2irp/Nia) = aI (t + 2irp/Nia) = A for 
some h,p: \p — p\ + \n — h\ ^ 0. Considering the above inequality, we obtain that the equality 
holds for some p ^ p. The converse part of the lemma is trivial. ■ 

To obtain a new non-resonance set x 2 (A) we remove from 3C 2 x*(A), shortly speaking, a neigh- 
borhood of its intersections with the whole Si (A). More precisely, we remove from 3C 2 Xi(A) the 
following set: 

fti(A) = {T€K 2 xi(A):an,n€N, P,P G P, P + V ■ A„ 1} (r + 2n P /N ia ) = A, 

r + 2np/N ia G X J(A), |A«(r + 2^p/N ia ) - \f{r + 2ttp/N 1 o)\ < e t }, (3.6) 

where t\ = e~i kV 1 . We define X2(A) by the formula: 

X2(A)=X 2 (xt(A))\Oi(A). (3.7) 

3.2. Perturbation Formulae. Before proving the main result, we formulate the Geometric 
Lemma: 

Lemma 3.5 (Geometric Lemma). For an arbitrarily small positive 5, 75 < 21 — 11 — 16si and 
sufficiently large X, A > Xo(V,8), there exists a non-resonance set X2(A,5) C %2X\ such that: 

(1) For any r G \i, the following conditions hold: 

(a) There exists a unique p G P such that r + 2irp/Nia G Xi- 7 

(b) The following relation holds: 

xf\T + 2irp/N ia ) = k 21 , 

where X^ {t + 2irp / N\a) is given by the perturbation series \2. 1 b}) with a = 1, j being 
uniquely defined by t = r + 2irp/Nia as it is described in Part 2 of Geometric Lemma 
for the previous step. 

(c) The eigenvalue A^- (r + 2irp/Nia) is a simple eigenvalue of H^\t) and its distance 
from all other eigenvalues X^\t + 27rp/Nia), n G N of H\{t) is greater than e\ = 
e -k kVSl : 

\xf\r + 27rp/N ia ) - A^(r + 2n P /N ia )\ > e v (3.8) 

(2) For any r in the (eik~ 2l+1 ~ s ) -neighborhood in C 2 of X2, there exists a unique p G P such 
that t + 2irp/N\a is in the (e\k~ 2l+1 ~ s ) -neighborhood in C 2 of xt an d 

\xf\r + 27rp/N ia ) - k 2l \ < e^k' 5 , (3.9) 

j being uniquely defined by r + 2irp/Nia as it is described in Part 2 of Geometric Lemma 
for the previous step. 

(3) The second non-resonance set X2 has an asymptotically full measure in xi in the following 
sense: 

L{^2X*l\X2)) ,-2-2 Sl (o in\ 

L(xl) < ■ (3 - 10) 

The proof of the lemma will be presented in section 13.31 



From geometric point of view this means that Xz{X) does not have self-intersections. 
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Corollary 3.6. Ifr belongs to the complex (e\k~ 2l+1 ~ 5 ) — neighborhood of the second non-resonance 
set X2(A, 5), then for any z lying on the circle C2 = {z : \z — k 2l \ = e\/2}, the following inequalities 
hold: 

|| (j ^(D (r) _ 2; )-i|| < ±, (3.11) 

||(iJ«(r) - < ^ , co = 326x62. (3.12) 

ei 

Corollary is proven in Appendix lA.il 

Remark 3.7. Note that every point 2-Km/N\a (m G Z 2 ) of a dual lattice corresponding to the 
larger periods N\ai, N\CL2 can be uniquely represented in the form 2irm/Nia = 2nj/a + 2irp/Nia, 
where m = N\j +p and 27rj'/a is a point of a dual lattice for periods 01, 02, while p G P is 
responsible for refining the lattice. 

Let us consider a normalized eigenfunction i/j n (t,x) of H^(t) in L 2 {Q±). We extend it 
quasiperiodically to Q2, renormalize in L2(Q2) and denote the new function by tp n {T, x), r = D^i. 
The Fourier representations of ip n (t,x) in L>z(Qi) and tp n (r,x) in L2{Q-2) are simply related. 
If we denote Fourier coefficients of ip n (t,x) with respect to the basis of exponential functions 
|Q 1 |-i/2 e <(t+27ri/a,a.) > j e Z 2 ? in l 2 {Qi) by C nj , then, the Fourier coefficients C nm of i> n {r,x) with 
respect to the basis of exponential functions |Q 2 |- 1 / 2 e *< T + 27rm / Ar i a . x > i m g Z 2 5 i n L 2 (Q2) are given 
by the formula: 

_ |C n j, if m = jNi +p; 
[ 0, otherwise, 

p being defined from the relation t = r + 2irp/Nia, peP. Hence, matrices of the projections on 
ip n (t,x) and ip n (r,x) with respect to the above bases are simply related: 



Cm 



= I (- B n)mm, if m = jiVi + p, m = jN\ + p; 
13 I 0, otherwise, 



i? n and E n being projections in L2{Q2) and L2{Qi), respectively. 

Let us denote by Ey (r + 2irp/Nia) the spectral projection Ej\a, t) (see (|2.17[l ) with a = 1 
and i = r + 2irp/Nia, "extended" from L,2{Qi) to L2{Q 2 )- 



By analogy with (|2.12|) . ()2.13j) . we define functions g)- {k,r) and operator- valued functions 
as follows: 

gf\k,r) = t^Tr£ ((ffW(r) - z)~ V 2 ) r <fe, (3.13) 



Gr (k, t), r = 1, 2, • • • , as follows: 



IC 3 

-Jtt- V, rnh, . 



G< 2 >(fc,r) = f ((^(r) -z)-V 2 )' (^(r) -z)" 1 ^. (3.14) 



We consider the operators H a = + aW2 and the family H a (t), r G K2, acting in L,2(Q 2 )- 

Theorem 3.8. Suppose r belongs to the (eik~ 2l+1 ~ 5 )-neighborhood in K2 of the second non- 
resonance set X2(A, 5), < 75 < 21 — 11 — 16si ; t\ = e~± kn 1 . Then, for sufficiently large X, 
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(2) 

A > Ao(V, 5) and for all a, < a < 1, there exists a unique eigenvalue of the operator H a (r) in 
the interval S2(k) := (k 21 — ei/2, k 21 + e\/2). It is given by the series: 



\f(a,r) = Af (r + 2-jrp/N x a) + ' \k,r), j = j+p/N x , (3.15) 



r=l 

?2 



converging absolutely in the disk \a\ < 1, where p G P and j G Z are described as in Geometric 

(2) 

Lemma \S.h\ The spectral projection corresponding to (a, t) is given by the series: 



j 

oo 



E?\a,r) = Ef>(T + 27rp/N ia ) +Y,^ r G {2) (k,T), (3.16) 

r=l 

which converges in the trace class Si uniformly with respect to a in the disk \a\ < 1. 
The following estimates hold for coefficients g^\k,r), G^' \k,r) , r > 1; 



9r 



W(k,r)\ < ^i(4e?f, llG^^r)^ <6r(4e?r. (3.17) 



Corollary 3.9. T/ie following estimates hold for the perturbed eigenvalue and its spectral projec- 
tion: 



\r> (a, t) - \y> (r + 27rp/N 1 a) < 12aef, (3.18) 
E?\a,r) - Ef ) (T + 27rp/N 1 a) < 48ae?. (3.19) 



i 



(2) 

Remark 3.10. The theorem states that A^ (a,r) is a single eigenvalue in the interval E2(k,5). 
This means that \\?\a,T) - k 2l \ < ei/2. Formula (|3.18|) provides a stronger estimate on the 
location of \~ 2 \a, r). 

3 

Proof. The proof of the theorem is analogous to that of Theorem 12 . 31 and it is based on expanding 

the resolvent (H a 2 \r) — z)~ l in a perturbation series for z G C2. Integrating the resolvent yields 

(2) 

the formulae for an eigenvalue of H a and its spectral projection. In fact, it is obvious that 

{H a 2 \r) - z)- 1 = (# (1) (r) - z)-\l - aA 2 y\ A 2 := -W 2 {H^(t) - z)~ l . (3.20) 
Suppose z G C2. Using Corollary 13.61 and the estimate 1 1 W2 1 1 < e~ kVSl = ef, we obtain: 

||(jj(i)( T )_ z )-i||<±, p 2 ||<ffiM <4e ?<l. (3.21) 

The last inequality makes it possible to expand (/ — a^) -1 in the series in powers of dis- 
integrating the series for the resolvent, we obtain formulae (|3.15(l . (|3.16|) . Estimates (|3.17|) follow 
from the estimates (j3.21|) . | 

Next, we show that the series (|3.15j) . (|3.16j) can be extended as holomorphic functions of r in 
a complex neighborhood of xi\ they can be differentiated any number of times with respect to r 
and retain their asymptotic character. 



SPECTRAL PROPERTIES OF A POLYHARMONIC OPERATOR WITH LIMIT-PERIODIC POTENTIAL 21 

(2) (2) 

Lemma 3.11. The following estimates hold for the coefficients g). {k,T) and G r (k,r) in the 
complex (\e\k~ 2l+1 ~ 6 ) -neighborhood of the non-resonance set X2- 

\T(m)gi 2 \k,T)\ < ml -3 • 2 2j - 1+ H e f +i-M A .M(«-i+*) ) ( 3 .22) 
\\T(m)G^ (k,r)\\i < ml • 3r • 2 2r+1+ H e f "NlfcMCa-l+fl (3 . 23) 
Proof. Since 1)3.11(1 is valid in the complex (ei£;~ 2 ^ +1_5 )- neighborhood of the second non-resonance 

(2) f2l 

set, it is not hard to see that the coefficients g)- (k,r) and Gr (k,r) can be continued from the 
real (eiA;~ 2i+1 ~ 5 )-neighborhood of r to the complex (ei/c~ 2 ' +1_<5 )-neighborhood as holomorphic 
functions of two variables and inequalities 1)3 .17|) are hereby preserved. Estimating, by means of 
the Cauchy integral formula, the value of the derivative with respect to r in terms of the value 
of the function itself on the boundary of the (5 ghborhood of r (formulas 1)3.17(1 ). 

we obtain (l3~22l and (l3~23l . ■ 

From this lemma the following theorem easily follows. 

Theorem 3.12. Under the conditions of Theorem VJ. 8\ the series \3.15\) . \3.1b]) can be continued 
as holomorphic functions of two variables from the real {e\k~ 2l+1 ~ 5 ) -neighborhood of the non- 
resonance set X2 to its complex (e±k~ 2l+1 ~ s ) -neighborhood and the following estimates hold in the 
complex neighborhood: 

T(m) (\f\a,r) - \f\r + 2^p/N ia )) | < aC m e^ W fcH(W-l+«) f (3.24) 

|r(m) (Ef ] (a, t) -Ef\r + 27iy/JVio)) || < aC m el' H k^ 2l ~ l+s \ (3.25) 

here and below C m = 48m!2' m l . 
Corollary 3.13. 

V\f\a,T) = 2lk 2l - 2 k + o(k 21 - 1 ), k=p j (T + 2irp/N 1 a). (3.26) 

T(m)\f\a,T) = 0(k 21 - 2 ) if \m\ = 2. (3.27) 
The next lemma will be used in the third step of approximation where the operator (r) = 

(2) 

H\ (t) will play a role of the initial (unperturbed) operator. 

Lemma 3.14. For any z on the circle C2 and r in the complex (e\k~ 2l+1 ~ s )— neighborhood of 

X2, 

||(ff(!0( T )_ z )-l||<JL. (3.28) 

Proof. Considering the Hilbert relation 

(tf (2) (r) - z)- 1 = (F (1) (r) - z)- 1 + (H {1) (r) - z)-\-W 2 ){H^{t) - z)~\ 
and the estimate 1)3.11(1 . together with the estimate || W2 1| < (-1, we obtain: 

ll^ (2) (-) " ^ ; TiKuVVLn * 2|l( ^ (1)(r) " * f (3 - 29) 

1 — \\(HW (t ) — z ) L W2\\ £ i 
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3.3. Proof of the Geometric Lemma. 

3.3.1. Proof of Statement 1. At the end of Section 13.1.11 (page EJ) we defined the second non- 
resonance set, X2(A) C K2, by the formula X2(A) = %2 (x*) \ the set ^1 being given by 
(|3.6j) . Suppose r G X2(A). Since X2(A) C 3C 2 (Xi)> there is a p G P such that r + 2-np/Nxa <E x*- 
By Lemma 12.121 there is an eigenvalue \u\t + 2np/Nia) = A, which is given by perturbation 
series (|2.16|) . By the same lemma An (r + 2np/Nxa) is a unique eigenvalue of (r + 2irp/N\a) 
in e x = (k 21 - k 2l ~ 2 -^- s , k 21 + k 21 - 2 -^- 5 ). This means 

A n 1} (r + 2irp/Nia) - a! 1} (t + 2irp/Nxa) > k 2l ~ 2 -^- s , when n / n. (3.30) 

Using the definition of f2i, we obtain: 

|A^(r + 2-Kp/Nxa) - \§\t + 27rp/JVia)| > e x when p / p. 

Combining the last two inequalities, we obtain Statement 1(c). It remains to prove that p is 
defined uniquely by the relation r + 2irp/Nia G Xi- Suppose it is not so, and there is p ^ p: 
t + 2irp/N\a £ X*- Applying Lemma T2. 121 to both points in Xi, we obtain A n ^(r + 2irp/Nia) = 

A^ {t+2ttp / N\a) , the eigenvalues being given by perturbation series. The last equality contradicts 
to (l3~3ffl) . 

3.3.2. Proof of the Statement 2. If r is in the complex (ei/c~ 2i+1 ~ 5 )-neighborhood of X2> then 
there exists a tq in X2 such that |r — tq\ < e\k~ 2l+1 ~ & and tq satisfies Statement 1 of the lemma. 
Obviously r + 2irp/Nia is in the (/c _1_4si_2<s )-neighborhood of xi) where Theorem 12.31 holds. By 
Corollary ECU 



\f ) (T + 2irp/N 1 a)-\f ) (T + 2iTp/N 1 a) < sup f V T A} i; (f + 2irp/N ia ) \r - r | = C^iAr" 5 ). 
Considering that A^tq + 2irp/Nxa) = k 21 , we get (|33|) . 

3.3.3. Determinants. Intersections and Quasi-inter sections. Description of the set Sli in terms of 
determinants. We considered self-intersections of Si (A) belonging to OC2X1 ( see Definition 13 .21 and 
Lemma 13.41 page 117(1 . Here we describe self- intersections as zeros of determinants of operators of 
the type I + A, A G Si (see e.g. (21]). In fact, let us represent the operator (H^\t) - X)(H (t) + 
A) -1 in the form / + A\\ 

{H^(t) - X)(H (t) + A)" 1 = / + A x {t), Ax{t) = (Wx - 2X)(H (t) + A)" 1 . (3.31) 
Obviously, Ax(t) G Si. 

Remark 3.15. From properties of determinants and the definition of Si (A) it easily follows that 
the isoenergetic set Si (A) is the zero set of det^J + Ai (i)^ in K\. 

Remark 3.16. Matrices H (t), H^(t) and Ax(t) can be extended as holomorphic functions of 
two variables to a vicinity of R 2 . If y = Pj(0) + t, then matrices Ho(y), H^\y) and A x {y) differ 
from Ho(t), H^\t) and Ax(t) only by a shift of indices: Ho(y) mm = Ho(t) m+ j ym+ j, etc. Hence, 
IK^W^-z)" 1 !! = {{(H^it) - z^W and det(J + A x $)) =det(/ + Ai(tj). 
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Next, we recall that the set T>\(X) can be described in terms of vectors 5((<p) = >c{\,V)V, 
v = (cos ip, sin if), ip G 0i, see 1)2 .33(1 . (|2.42[) . By definition, Xi(A) = OCiDi(A). By Lemma 
12.121 page El Xi(A) does not have self- intersections, i.e., for every t G Xi(A), there is a single 
>c(if) G Di(A), such that t = %\?c{if). Next, if r € 3<C2X*(A), then there is p G P such that 
r + 2-Kp/N\a G X*(A). Note that p is not uniquely defined by r, since JC2X1W ma Y have self- 
intersections. Hence, every r G X2X1W can be represented as r = %2%(ip), here x{ip) is not 
necessary uniquely defined. The next lemma describes self-intersection of Si belonging to X2X1W 
as zeros of a group of determinants. 

Lemma 3.17. If r is a point of self-intersection of Si (A) belonging to 3C2X*(A), then r can be 
represented in the form r = %2^( ( p)> where if G @i and satisfies the equation 

det(7 + Ai(y(f))^j = 0, y(f) = 2(<p) + b, %=2irp/N x a, (3.32) 

for some p G P\ {0}. Conversely, if is satisfied for some p G P \ {0} ; then r = 3C2i?((/j) 

is a point of self-intersection. 

Proof. Suppose, r is a point of self-intersection of Si (A) belonging to 3C2X*(A). By Lemma l3~H 
page El there is p G P (here we use p instead of p in Lemma EH}, such that r + 2np/Nia G x*> 
and p G P, p 7^ p, such that r + 2irp/Nia G Si (A). This means that 

det(/ + ^i(r + 27rp/iVia)) = and det/ J + Ai(r + 2vrp/A^ia)) = 0. (3.33) 

By Lemma [2. 131 page !15l there is a single x(tf) G Di(A): %ix = r + 2wp/Nia (this means 3C2>? = 
r). According to Remark l3~T51 and formulae (|335|) . detfl-Mi(x)J = 0, det (/-Mi (>? + &)) = 0, 

where b = [2ir(p — p)/Nia\ mod K . Thus, we have obtained (|3.32|) . 

Suppose now that (|3.32|) holds for some p G P. Let us show that r = %2>c{if) is a point of 
self-intersection. Let p G P: r + 2np/Nia = %iH; such p, obviously, exists, since 3Ci>? G Ki, 
t G K 2 . Since G Di(A), we have r + 2-wp/Nia G x* C Si(A). Therefore, there is n G N: 

A^(r + 2np/Nia) = A. (3.34) 

Next, it follows from (|3.32jl and Remarks 13.15113. 1 51 that %iy(f) G Si (A). Hence, there is f G K2, 
p G P such that f + 2'npjNia = %iy(<f) G Si(A). This means that there is an n G N such that 

\ { }\t + 2-Kp/Nia) = A. (3.35) 

Note that t — f = ( X2V{ i p) — ^2^(f) = ^26 = 0, i.e., r = f. Considering this equality together 
with (|3.34jl . ()3.35j) . we see that r is a point of self-intersection, provided that p 7^ p. It remains 
to show that p 7^ p. In fact, r + 2wp/Nia = %ix, r + 2np/Nia = %iy, hence 2tt(p — p)/Nia = 
%i{jj — k) = 3Ci6 = b 7^ 0, the last equality holding, because b G K x . By Lemma l3~H r is a point 
of self-intersection of Si belonging to 3C2X*(A)- ■ 

Definition 3.18. Let <Pi be the complex {k~ 2 ~ isi ~ 25 ) -neighborhood of 61. 

By Lemma r2.ll! i?(v?) is an analytic function in <Pi, and, hence, det( I+A^y^ip)^ , y = x(f)-\-b, 
b G Ki, is analytic too. 
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Definition 3.19. We say that ^ € #1 is a quasi-intersection of OC2X1 with Si (A) if the relation 
(ET321) holds for some p G P \ {0}. 

Thus, real intersections correspond to real zeros of the determinant, while quasi-intersections 
may have a small imaginary part (quasi-intersections include intersections). 
Next we describe the non-resonance set in terms of determinants. 

Lemma 3.20. If r G then r = %2^( i p) where <p £ Ox and satisfies the equation 

det f + k2l J = 0, V{V) = %0 +6, b = ^p/N x a, (3.36) 

/or some p G P \ {0} and |e| < ei. Conversely, if \3. 36\) is satisfied for some p G P \ {0}, and 
|e| < ei, i/ien r = %2>c{(p) belongs to Q\. 

Proof. The proof is analogous to the proof of the previous lemma after we notice that the set of 
points r satisfying the inequality 

An (r + 2itp/N\a) — A? (r + 2irp/Nia) < e\ can be described as the set of points satisfying 

\&(t + 2irp/N 1 a) - \^\t + 2irp/N 1 a) = e for some |e| < e 1 . | 

Let us discuss separately the case W\ =0. If cp is a point of intersection, then ()3.32j) means 
that 

\2(<p) + 2nj/a + 2Trp/N ia \ 2 = k 2 (3.37) 
for some j G 1? and p G -P, and vise versa. Thus, (|3.32j) means that x((/j) is a point of intersection 
of the circle of the radius k centered at the origin (|x((/j)| 2 = k 2 when W\ = 0) and the circle 
of the same radius centered at a point 2irj/a + 2np/Nia of a dual lattice corresponding to the 
periods iViai, N\02- If there is a quasi-intersection, then (|3.37|) holds for some ip G <&\. We will 
refer to such ip as a point of quasi-intersection of two circles. Geometrically, such terminology is 
justified, since validity of (|3.37j) for some <p with a small imaginary part means, that two circles 
come close together, even if they do not intersect. 

We denote by u\ the set of (p G 61 corresponding to fii, i.e., u\ = {ip G 0i : X2>?i(y?) G Sli} C 
[0,2tt). 

3.3.4. Complex resonant set. Further we consider a complex resonance set wj[(A), which is the 
set of zeros of the determinants (|3.36|) in <Pi (p G P \ {0}, |e| < ei). By Lemma ?d.2Ul uj\ = 
fl 61. We prefer to consider quasi-intersections, instead of intersections, and the complex 
resonance set, instead of just the real one, for the following reason: the determinants (|3.32|) and 
()3.36|) . involved in the definitions of quasi-intersections and the complex resonance set lo\ are 
holomorphic functions of <p in (Pi. We can apply complex analysis theorem to these determinants. 
Rouche's theorem is particularly important here, since it states stability of zeros of a holomorphic 
function with respect to small perturbations of the function. We take the determinant (|3.32|) as 
a holomorphic function, its zeros being quasi-intersections: the initial determinant corresponds 
to the case W\ = 0, the perturbation is obtained by "switching on" a potential W\. Since there 
is no analog of Rouche's theorem (e.g. see |24| ) for real functions on the real axis, introducing 
the region (Pi and analytic extension of the determinants into this region is in the core of our 
considerations. We will use also a well-known inequality for the determinants (e.g. see |21|): 

det(I + A) - det(I + B) < \\A - fl||iexp(||4||i + + lV A, B G S x . (3.38) 
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Note that oj* = ^-> p ep\{o}^i p , where u)\ p corresponds to a fixed p in ((SHHU); and similarly, u\ = 
^ipeP\{o} u i,p- We fix p G P and study each u* >p separately. We start with the case W\ = 0. The 
corresponding determinant (|3,32|) is 

det (l + A) (yb (¥>))), I + Av{Uv)) = (^o(yoM) - a) (# (z7o(¥>)) + a)~\ (3-39) 

yo(<p) = k(cos <p, sin ip) + 6. This determinant can be investigated by elementary means. We easily 
check, that the number of zeros of this determinant in <I>i does not exceed cok 2+2si (Proposition 

l3~2Tl Corollary l3~22l . The resolvent (H (y (<p)) - k 2l ^j has poles at zeros of the determinant. 
The resolvent norm at <p £ <&\ can be estimated by the distance, which ip has to the nearest zero 
of the determinant ( Lemma 13 .23j) . Next, we introduce a set 0(b), which is the union of all disks 
of the radius r = &- 4 -6' s i-3<5 surrounding zeros of the determinant (|3.39f) in <P\. Obviously, any 
(p £ ^1 \0(6) is separated from zeros of the determinant (|3.39|) by the distance no less than r. This 

estimate on the distance yields an estimate for the norm of the resolvent ^Ho(yo(p)^ — k 2l%s j , 

when ip £ <P\ \0(b) (Remark l3.28j) . Further, we introduce a potential W\. Our goal in this section 
is to prove that the number of zeros of each determinant (|3.36f) is preserved in each connected 
component T(b) of 0(6), when we switch from the case W\ = 0, A\ = Aq to the case of nonzero 
W\ and from e = to |e| < ex- We also show that estimates for the resolvent are stable under such 
change when ip £ <P\ \0(6). We "switch on" potential W\ in two steps. First, we replace yb(¥>) by 

y(ip), i.e., we consider det(^I + Ao(y(p)^ and (^Ho(y(p)^ — k 2l ^j in <Pi. We take into account 

that y(<p) — yo(p) is small and holomorphic in <&\ (Lemma 12.11(1 . use (|3.38j) on the boundary of 
r and apply Rouche's theorem. This leads us to the conclusion, that the number of zeros of the 
determinant in r(6) is preserved, when we replace yo(<p) by y(p) (Lemma l3.29|) . Applying Hilbert 
relation for resolvents, we show that estimates for the resolvent in <Pi \ 0(b) are also stable under 
such change (Lemma l3.29l Corollary 13. 3U|) . In the second step we replace Ho(y(p)) by (y(<p)) 
and prove similar results (Lemmas 13.311 13.34[ Corollaries I3.32( I3.36JI . Further we consider the 
determinant (|3.36|) . Obviously it is equal to (|3.32j) when e = 0. Since, e is very small, its properties 
are similar to those of (|3.32j) . In particular, it has the same number of zeros (for a fixed e) as 
(|3.32|) at each connected component T(b) of 0(6) (Lemma 13.38(1 . From this lemma we see that 
u* C 0(6), 6 = 2-irp/Nia. Taking the real parts of the sets, we conclude: ui jP C 0(6) n ©1. We 
show that 0(6) is formed by no more than cok 2+2si disks of the radius r = ^ _4 -6si-3(5^ p rom 
this we easily obtain, that the total length of uj\ )P does not exceed k~ 2 ~ 4si ~ 3S . Since the set 
P contains no more than 4k 2s2 ~ 2si elements, S2 = 2s±, the total length of S7i does not exceed 

fc -2-2si-35_ Thig est i mate 

is proven in Section [3.3.51 
Let us introduce new notations: 

0* = UpgpuojQpTrp/JVia), $2 = $1 \ 0*, (3.40) 

^1 being given by Definitions 13.181 Note that the set 0* is formed by small disks centered at 
points tfmtf ^ ^1 j w hich are quasi-intersections of the circle of the radius k around the origin with 
the circles of the same radius k around points p m (0) + 2np/Nia, p 7^ 0, that is around points 
of the dual lattice corresponding to the periods Aqai, N\CL2- We will show (Corollary I3.39|) that 
uj\ C 0* and uj\ C 0* PI ©1. From now on we call 0* the first complex resonance set in <P\ and 
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Figure 7. The set 



0* fl 0i the first resonance set in 0i. Obviously, to obtain <I>2 we produce round holes in each 
connected component of <&\. The set has a structure of Swiss cheese (Fig. EJ; we will add 
more holes of a smaller size at each step of approximation. 

In the definition of the non-resonance set we take b = 2-Trp/Nia, p ^ 0. Let us consider now 
arbitrary b € K\ and investigate 0(b). 8 We denote by &o be the distance from b G ifi to the 
nearest vertex of i^i : 



&o = min 

m=(0,0),(0,l),(l,0),(l,l) 



6 — 2nm/a 



(3.41) 



Further we assume &o 7^ 0, i.e., b is not a vertex of K\. We denote by a narrow strip in the 
complex plane around the interval [0, 2ir): 

$0 = {tp : Slip £ [0,2vr),|^| < A;- 2 ~ 4si } . 

Obviously, <I>i C ^o- Let \yo(p) — xq\ 2 be an analytic function of p defined by the formula: 
\yo(<p) — xq\1 := (kcosp + bi — xqi) 2 + (ksinp + &2 — ^02) 2 > being a vector in M 2 . Note 

that |y (</?) - b\l = k 2 cos 2 p + /c 2 sin 2 99 = k 2 . Let |y (^) - »o|*' := (\yo(p) ~ »o|*) • Obviously, 
det^/ + ^(ybC^))! = if and only if |yo(y) +Pm(0)| 2i = k 21 for some m € I?. 

Proposition 3.21. Suppose b is not a vertex of K\. If p m (0) < 4k, then the number of solutions 
of the equation 



Wo(^)+Pm(0)\l l = k 21 (3.42) 
in ^0 does not exceed two. The solutions of ar &, in fact, the solutions <p^ of the equation 

Wo(<P)+Pm(0)\l = k 2 . (3.43) 

If the distance from tp £ <Pq to the nearest solution p^ of j3.43j ) is greater than en, en > 0, then 

(3.44) 



\m(v)+p m mf -k 21 



> ciboe^k 21 1 , 



c\ being a positive number depending only on I. If p m (0) > 4k, then the equation jS.4^j) does not 
have solutions in&o. Moreover, 

\M<P) +Pm(0)\l l ~k 2l \> ^P 2 1(0). (3.45) 



We need arbitrary 6 for next steps of approximation. 
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The proof is elementary and based on the fact that (|3.43j) is a quadratic equation with respect 
to j*. 

Corollary 3.22. For each b, the number of (p € satisfying \3.^ for at least one m 6 Z 2 does 
not exceed cok 2+2si , cq = 326162. 

Corollary VJ.2°A The equation (|3,42l) can have solutions in <Pq only if p m (0) < 4k. The number of 
such m's does not exceed 166i&2^ 2+2si • For a fixed m there exist no more than 2 solutions, tpf^. 
This gives the statement of the corollary. ■ 

According to the definition at the end of Section l3.3.31 points (p^ are points of quasi-intersections 
of the circle of the radius k centered at the origin and the circle of the same radius centered at 
-p m (0) - b. 

Lemma 3.23. If (p G satisfies the inequality min mgZ 2 \ip — ip^\ ^ e o for some < eo < 1, 
then the following estimates hold: 

(tfo(yoM) - k^Y'j < Cl k- 2l+ %% 2 , (3.46) 

(ilo(yoH) " k 2l Y\ < dk-^+^ef. (3.47) 

Proof. Estimate (|3.46|) immediately follows from (|3.44|) and l)3.45|) . Considering that -fTo(?/o( < ^)) 
is a diagonal operator and summarizing (|3.44j) . (|3.45j) over m, we arrive at (|3.47|) . | 

Now let us construct contours in C whose distance from all the solutions of the equations 
\i}o{}P) +Pm(0)\ 21 = i; 2 ', m G Z 2 , is at least /j- 4-68 ! -35 . "We recall that there are no more than two 
solutions in <Pq if p m (0) < 4k, and there are no solutions if p m (0) > 4k (Proposition 13.21]) . 

Definition 3.24. For each 99+ (or <^~) in <Po we make an open disk of the radius ^~ 4 ^6si-3<5 
centered at (or c/?~) and denote it by 0+ (or 0~). Let 0' the union of all such disks: 
0' = U m) ±0^. We denote by the part of 0' related to 3>\. Namely, is the union of connected 
components of 0' not being completely outside To stress that depends on b, we often will 
write 0(6). 

Remark 3.25. Obviously, 0(6) consists of no more than cok 2+2si disks (see Corollary I3.22|) . 
Thus set 0* given by 1)3.40(1 contains no more than 4cok 2Jr2s2 disks, since P contains less than 
4k S2 ~ Sl elements. 

Definition 3.26. Let T(6) be a connected component of 0(6) and 7 = dT. 

Remark 3.27. The set 0(6) consists of several components T(6). Since the number of points <p^ 
in ^0 does not exceed cok 2+2si , the size of any T is less than cok 2+2si .2/c _4_6,Sl_3<5 = 2cok~ 2 ~ 4:Sl ~ 3S . 
Clearly it is much smaller than the size of (see Definition I3.18|) . Some T(6) are strictly inside 
<Pi, while others can have their parts outside. However, any T(6) is in the (2cofc _2_4si_3<5 ) — 
neighborhood of <&\. This neighborhood has properties completely analogous to those of #1, since 
c/c _2 ~ 4si_3<5 << fc~ 2 ~ 4si ~ 2<5 , the latter value being the size of <Pi. Hence, even if a component 
r(6) is not strictly inside it has the same properties as those inside $\. Further we consider 
r(6) C <2>i, 0(6) C #1. 
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Remark 3.28. For any ip G <Pi \ 0(6) and any tp^ the estimate \tp — ip^l > k 4 6si 3,5 holds. 
Therefore, the estimates (jTHJ), and ({3~T7|) hold with e = k~ 4 ~ 6si ~ 35 for such <p. 



Obviously, zeros of det(^I+Ao(y(p))j are described by the equations \%{p) + b+p m (Q)\^ L = k 21 , 

m G Z 2 , G ^i, which can be rewritten as \yo((f) + h\(p>)v + Pm(0)| 2 ' = fc 2i , the functions 
ii(ip) = x(A, z?) and h\((p) = hi(\,is) being defined in Lemma 12.111 It is easy to see that the 
vectors >c(p>) satisfying this equation with a real p are just intersections of Di(A) with the circle 
of radius k centered at —b — p m (0). Since y(p) differ from yo(ip) by a small value h\(ip), the 
estimates IpOijl . (|3~l5j) and p^TH]) (|3~1?|) with e = A;~ 4 ~ 6si ~ 35 are preserved for ip G #i \ 0, when 
we replace yo(p) by y(p). We prove this result in the following lemma. 

Lemma 3.29. If b > fc7+8«i-2<yo+6<5 ; 7o = 2 / _ 2 - 4si - 25, then: 

(1) T/ie number of zeros of det(/ + Ao(y(p)) I in any T(b) is the same as that of detl / + 



^o(ybM))- 

(2) The following estimates hold for p G <fr\ \ 0(b): 

< cb^k- 21 ^ 12 ^ 66 , 

< cb^k- 2l+n+14s ^ 6S . 

Corollary 3.30. There are no zeros of det( I + A (y(p)) ) in <I>\ \ 0(b). 



H (y(p))-k 2L 
H (y(p))-k 21 



(3.48) 
(3.49) 



Lemma T*TM Note that det(/ + A (y (p)) j = if and only if \y (p) + Pm(0)\ 21 = k 21 for at least 

one m G 1? . The analogous statement holds for det^/ + j4o (y(p))^j • It is easy to show (Appendix 
|A~2|) that 



y6(<p) +Pm{0)\f - \y(<p) +p m (0)\ 



21 



1 

< 2 



\m(<p)+Pm(o)\ 2 J -k 21 



(3.50) 



for any ip in #x \ 0(6) and m G Z 2 . By Rouche's theorem [2D, the number of zeros of \yo(tp) + 
Pm(0)|f — k 21 in every T(6) is the same as that of \y(ip) + Pm(0)\f — k 21 . Therefore, the first 
statement of the lemma is proved. 

To prove the second statement, we use the Hilbert relation: 



(flo(ifaO) 



k 



21 



H Q (m(p>)) - k 



21 



+ (H (y(ip)) - k 21 ) (H Q {y Q (ip)) - H (y(p)) )(H {y (<p)) - k 



',21 



By Lemma with e = k~ 4 ~ 6si ~ 3S , 

' H (y (p)) - k 21 



Ho[yo(^>)) - k 



21 



< cb Q l k- 2l+9+l2s ^ 8 , 



< Cb lfc-2J+ll+14.i+W 
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We easily obtain from (j3.5()j) that 



H {y (<p))-Ho{y(<p)))(H (y {<p)) 



.,21 



< 1/2. 



(3.51) 



Now, using a standard argument as in Lemma 12.91 we arrive at (|3,48|) and Q3.49J1 . ■ 

Lemma 3.31. If bo > k~ 2l+9+12si+7S , then the following estimates hold for any tp £ #i \ 0(b): 

\(H^(y(<p)) - k 21 )-^ < cb^k- 21 ^ 12 ^ 68 , (3.52) 
\(H^(y( V )) - k 21 )' 1 ^ < cb Q ijfe-a+u+Wi+W (3.53) 
Corollary 3.32. The determinant of I + A\(y((p)) has no zeros in <P\ \ 0(b). 
Proof. We use Hilbert relation (^H^ (y(<p)) — k 2l ^j = 

(h (m) -k 2i y\ (H^(m) - k 2i y\-wi)(H,(m) 

Since b > k- 2l + 9+12si+7S > fc7+8 S i-2 70 +6<5 (70 = 21 - 2 - 4s x - 25), the previous lemma can be 
applied. It yields: 

Using a standard argument as in Lemma 12.91 we arrive at (|3.52j) and (|3.53jl . | 



Remark 3.33. Obviously (1332]) and (13331) hold for any fl^ 
in a. 

Further, we will use the notation ^ = TiTTT 1 for a pair of operators Ti, T2. 
Lemma 3.34. J/ bo > fc- 2i + 9 + 12s i+ 75 ; ^ e number of zeros of det( I + A\ 

T(b) is the same as that of det (^1 + Ao(y(<p) 

< H (y(<p)) + a x W x - k 21 



Hq + aWi , < a < 1 uniformly 



^n any 



Proof. Let us consider det 



det 



^o(^))+^i-fc 2 ' 
g (y(y))+a 1 ^i-fc 2 ^ / 



det I + 



, where < «i < 02 < 1. Obviously, 
(ai - 0:2)^1 \ 



^(^))+«2^1-i 2i 



Using (|3.53j) . which is uniform with respect to a, we obtain, by Rouche's theorem, that this 
determinant has the same number of poles and zeros if \a 2 — ol\\ is sufficiently small: |q2 — ol\\ < 
bok 2l ~ 11 ~ 14:Sl ~ 7S . We make a finite number of steps between a = and a = 1. Indeed, 



det 
det 



H {y(<p)) +a 2 W 1 -k 2 \ 
' H (y(<p)) + - k 21 ' 

H (y(<p)) +a 3 W 1 -k 21 



has the same number of zeros and poles and so does 



. Both determinants can be represented as fractions of holomor- 



phic functions: arm —j^t, respectively. The number of zeros being the same for the 

I2VP) 92m 
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/ H (y((f)) + axWt - k 2l \ f m 

pairs fi, fo and oi , q 2 - Note that det — — 4 — = — — . Obviously, it also 

\H (y(ip)) + asW! - ) f 2 g 2 

( H (y{^)) + W 1 -k 2l \ 

has the same number of zeros and poles. Hence, det . has the same 

\ H (y((p)) - k 21 J 

number of zeros and poles inside T(b). Considering that the last determinant is the quotient of 
det ^1+ Ai^y ((/?)) ^ and det^/+ J 4o(y('/ 3 ))) and these two determinants are holomorphic functions, 
we obtain, that they have the same number of zeros inside T(b). | 

Lemma 3.35. (1) If bo > k^ 2l+9+12si+7S , then the number of zeros of 

det + A\ (y ((/?)) ^ in any T(b) is the same as that of det (^1 + Aq (ybCv 9 ))) • 

(2) The distance from a zero of det (^1 + Ai(^y(ip))^j to the nearest zero of det (^1 + Ao(yo((p))^j 
does not exceed r/2. 

Corollary 3.36. The total number of zeros of det ^/+yli (y ((/?)) ^ m0(6) does not exceed cok 2+2si . 

The Corollary immediately follows from the first statement of the lemma, the definition of 0(6) 
as a collection of disks surrounding zeros of det( / + Ao(yo(p)) ) an d Remark 13.251 page EH 



Corollary 3.37. The number of self-intersections of X2X1W does not exceed cok 2+2si and 
L(X 2X t(X)) = L(xl(X)). 

The latter equality holds, since %2X\(X) has only a finite number of self-intersections. 

Lemm,a, [,V. The first statement immediately follows Lemmas 13. 291 fPart 1) and 13 .341 Next, note 
that Lemmas 13. 291 - l3~3"H hold not only for 0(6), but also for a set 0(6), constructed from disks of 
twice smaller radius r/2 with the same centers, since all estimates in the lemmas preserved under 

such change up to some insignificant constants. This proves that all zeros of det^I + ^1(2/(99))^ 
in 0(6) are, in fact, in the smaller set 0(6), i.e., the second statement of the lemma holds. | 

Let us recall that tu\ p is the set of zeros in ^1 of the determinant (|3.36|) when |e| < ex, here 
p € P \ {0} and fixed. The set U\ jP is the real part of u* : u± )P = u\ p fl 0i. 

Lemma 3.38. If b > k- 2l+9+12si+7S , \e\ < b k 2l - u ~ lAsi - 75 , then 



HW(y(<p))-k 



21 



€ 



(1) The determinant of — — — -pr, — has no zeros in <Pi \ 0(6). 

H (y{tp)) + k 21 

{ (y(<p)) - k 21 - e\ 

(2) The number of zeros of det ._. ,. — in T(b) is the same as that of 

\ H (y(ip)) + k 21 



det 



'H^(y(cp))-k 



21 



^ H (y(<p)) + k 21 J' 
Corollary 3.39. Ifb = 2Trp/N ia , u* lyP C 0(6), u 1>p C 0(6) n 0i. 
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Proof. To prove Lemmas l3.31l and l3.34[ we used the estimate || Wi [| = O(l). Obviously, ||Wi— e|| = 
O(l). Hence, the analogous lemmas hold when we replace W\ by W\ — e, i.e., when we consider 
H^(y((p)) — k 2l — e instead of just (y(p)) — k 21 . Modified Lemmas 13 . 3 1 1 and 13 . 341 immediately 
yield the statement of the present lemma. ■ 

3.3.5. Proof of statement 3 of Geometric Lemma. By Corollary 13.391 u\ jP C 0(b) n Br, b = 
2np/N\a. It is easy to see that bo > Cbk~ S2 , Cb = |min{&i,&2j-- Such b satisfies conditions of 
Lemmas 13.291 - 15. 381 for any p G P\ {0}. Taking into account that that 0(6) is formed by no more 
than cok 2+2si disks of the radius r = j^-^-^ 13 ^-^ ; we obtain that the total length of uj\ p does not 
exceed 2coA;~ 2 ~ 4si_3<5 . Since the set P contains no more than 4k 2s 2-' 2s i elements and S2 = 2si, 
the total length of u\ does not exceed 8co/c~ 2 ~ 2si ~ 35 . Therefore, the length of f2i does not exceed 
207rcoA;~ 1 ~ 2si ~ 3<5 . Considering that Xi(A) has a length 27rfe(l + o(l)V we obtain Statement 3. 

(2) 

3.4. Non-resonant part of the isoenergetic set of H a . Let 52(A) be an isoenergetic set 
of the operator H^: S 2 (X) = {r G K 2 : 3n G N : \£\oc,t) = A}, here {\£\a, r)}£° =1 is the 

(2) 

spectrum of H a (t). Now we construct a non-resonance subset X 2 (A) °f ^2 (A). It corresponds 

to non-resonance eigenvalues \~. (r) given by the perturbation series (|3,15|) , We start with a 

definition of 1>i(X) nonres . Recall that \2 C 9C 2 Xi(A) ( see the Geometric Lemma) and X*(A) = 
DCiDi(A), see (|2.48j) . Hence, %2 C 3C 2 T"i(A). Let Di(A) nonres be the preimage of \2 in Di(A): 

nonres — 

{xG Di(A) :0C 2 xGx 2 }. (3.54) 
Lemma 3.40. XTie formula 0C2T>i(X) nonres = \2 establishes one-to-one correspondence between 

nonres 

and X2- 

Proof. Suppose there is a pair xl,x 2 G 1 ) i(A)nonres such that K 2 xi = ^2^2 = t, t £ \2- We 
introduce also t\ = %\%\ and i 2 = 3Cix 2 . The definition (|2.48|) of x*(A) implies that £i,t 2 G Xi(A), 
since xl, x 2 G Di(A) nonres C Di(A). Clearly, % 2 t\ = %2t2 = t and, hence, t\ = t + 2-irpi/Nia, 
t 2 = t + 2np2/Nia for some pi,p 2 G P. Now, by Part la of Geometric Lemma 1331 p\ = P2, and, 
therefore, t\ = t 2 . Next, by Lemma 12.131 xi = x 2 . I 

We define 22(A) as the set of directions corresponding to 2?i(A) nonres : 

3 2 (A) = {v G Bi(A) : xi (A, v)v G T>i(A) nonrefl }, 

where xi(A, v) is defined by Lemma [2.111 xi(A, z?) = xi (</?). Note that 2 2 (A) is a unit circle 
with holes, centered at the origin, and 22(A) C iBi(A). We denote by 82(A) the set of angles cp, 
corresponding to 22(A): 

G 2 (A) = {ip G [0,2vr) : (cos^sin^) G 2 2 (A)}, G 2 C 6i. 

Let x G T>i(X) nonres . By (|3.54|) . r = 3C 2 x G X2(A). According to Theorem 13.81 for sufficiently 

large A, there exists an eigenvalue of the operator H a (r), given by (|3.15f) . It is convenient here 

to denote a! 2) (q,t) by \( 2 \a,x). We can do this, since, by Lemma I3.4fl( there is one-to-one 

correspondence between x G I>i (A) now . es and the pair (r, j), x = 2nj/a + r. We rewrite (|3. 15f) 
in the form: 

00 

A( 2 V,x) = A«(x) + / 2 (a,x), / 2 (a,x) = ^a^( 2 )(x), (3.55) 

■r=l 
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(9) 

here gr (x) is given by 1)3.13(1 . The function /2(a,x) satisfies the estimates: 

|/ 2 (a,x)| <12aei (3.56) 
|V/ 2 (a,x)| < 2Aae\k 2l - 1+s . (3.57) 

By Theorems 13.81 and 13.121 the formulas Q3.55JI - (|3.57|) hold even in (ei/c~ 2i+1_<5 )-neighborhood 
of Di(\) nonres , i.e., for any x = xv such that V G 232(A) and |x — xi(A,z?)[ < e\k~ 2l+1 ~ s . We 
define ©2(A) as a level set for A^ 2 ^(a,i?) in this neighborhood: 

D 2 (A) := jx = xz7: i7 g B 2) |x - *l(A, ^)| < eiAT 2 ^ 1 ^, A( 2 )(q, x) = a} . (3.58) 

Next two lemmas are to prove that D 2 (A) is a distorted circle with holes. 

Lemma 3.41. For every V E H2 and every a, < a < 1, i/iere is a unique x = x 2 (A, i?) m i/je 
interval I 2 '.= [xi(A, z?) - eifc - 2Z+1— <5 , xi(A, z?) + ei&;~ 2 ' +1 ~ <5 ] sitc/i that 

A (2) (a,x 2 z?) = A. (3.59) 

Furthermore, 

|x 2 (A,i?) - xi(A,z?)| < 2ae\k- 2l+l . (3.60) 

Proof. Let us show that there exists a x 2 (A,z?) satisfying (j3.59|) . Since xii? G 2)i(A) when 
z? G 232(A), for all x in the interval I 2 , there is an eigenvalue \( 2 \a, xv) given by (|3.55|) . In 
particular, A( 2 )(0,xii7) = AW(xiz7) = A. Let L^(a,u) := {A = A^(a,xi7) : x G I 2 }. It 
is an interval, since \( 2 \a, xv) is a continuous function of x. Considering that A^ 1 ) (xi?) is 
a continuous function and VX^ixv 1 + o(l)), we readily obtain: £/ 2 )(0, V) D 

[/c 2Z — cie\k~ s ,k 21 + ciei/c -5 ], < ci 7^ ci(/c). Since A( 2 )(a,xz7) is continuous in x and (|3.56j) 
holds, we have L^ 2 \a,v) D [A; 2 ' - ciei*;-" 5 + c 2 ef,k 21 + ciei*;" 5 - c 2 ef], < c 2 / c 2 (fc). The 
inequality ef < eik~ 5 yields: A( 2 )(a,i?) D [A; 2/ - cieik~ s /2, k 21 + cieiAr 5 / 2 ]- Thus for every 
V G 232(A), there exists a X2(A, V) G I 2 such that A^ 2 ^(a, x 2 i>) = A. 
Using (|3.26|) . we obtain 

A-J ^ = 2/fc 2 ^ 1 (l + o(l)), (3.61) 

ox 

and this implies that \^ 2 \a, xv) is monotone with respect to x in /2- Hence, there is only one 
x 2 (A, z?) satisfying l)3.59j) in I 2 . 

Now let us estimate |x 2 (A, v) — xi(A, v)\. By (|3.55() . we have A( 2 )(a,xiz?) = A^fxiz?) + 
/ 2 (a,xii7) = A+/ 2 (a, xii7) . Using (|3.56j) and (|3.59|) . we easily obtain: | A^ 2 ) (a, x 2 v) — A^ 2 ) (a, x\v) | < 
YLae\. Applying the Mean Value Theorem at the left-hand side and using (|3.61|) . we arrive at 
the inequality: lk 2l ~ 1 \x 2 — x\\ < Y2ae\. The estimate (|3.60f) immediately follows. ■ 

Let us recall that <F 2 is the second complex non-resonant set given by (|3.4()|) . see Fig. 
Further, we use the notations: x 2 (ip) = x 2 (A, z?), h 2 ((p) = ^(9?) — ^li^p), >? 2 (y) = ^li^p)^. 

Lemma 3.42. (1) The set T> 2 {\) is a distorted circle with holes: it can be described by the 
formula: 

V 2 (\) = {xeR 2 :x=x 2 (v), ^6 2 (A)}, (3.62) 
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where x 2 (<^>) = ^(v 9 ) + h 2 (ip), ^li^) is the "radius" o/©i(A) and h 2 ((p) satisfies the 
estimates 

M < 2ae\k~ 2l+ \ ^ < ±ae\k 1+s . (3.63) 

(2) The total length of '82(A) satisfies the estimate: 

L(Bi \S 2 ) < 4irk- 2 - 2si . (3.64) 

(3) T/ie function >c 2 ((p) can be extended as a holomorphic function of 92 into the complex 
non-resonance set <P 2 , estimates h'J. 6°J\) being preserved. 

(4) The curve 2) 2 (A) has a length which is asymptotically close to that o/T>i(A) in the following 

sense: 

l(d 2 (A)) ^^(^(A^l + O^- 2 -^)). (3.65) 

Proof. (1) By Lemma 13.411 \h 2 \ = |x 2 — < 2ae\k~ 2l+1 . Let us prove the second estimate 
in (|3.63j) . Considering that A = A^ 2 ) (a, x 2 ((/?)), A = A^ , differentiating both 

equations with respect to (p and subtracting the results, we obtain: 

-~ — (a, x 2 ) — = F(A, p), (3.66) 

OX Oif 

F(\,<p)= (3.67) 

(VA« (x x ) - VA( 2 ) (a, x 2 )) ^ + (xiVA« (xi) - x 2 VA^ (a, x 2 ) , /?), 

where jl = (— sin ip, cos 99). It is easy to check f Appendix IA.3|) that 

\F(X,<p)\ < 27ak 2l+s e\. (3.68) 

Then, the second estimate in (|3.63(l easily follows from formulae (|3.66|) and estimates 
(t3~6TTl and d3~68l) . 

(2) Considering that Q\ is the set of angles corresponding to 1*1 (A) (Lemma 12. 11(1 and 2 is 
the set of angles corresponding to r Di(X) nonres , we obtain 

L(T) 1 (X)\'D 1 (X) n(mres ) = [ J + xiitp) 2 dtp. (3.69) 

V ' J0x\0i V 

Using (|2.43j) . we easily show that the integrand is close to k. Hence, 

L(l>i(A)\Di(A) nonres ) = L(0 1 \0 2 )k(l + o(l)). (3.70) 
Let us show that 

L(© 1 (A)\D 1 (A) nonres ) <2vrA:- 1 - 2si (l + o(l)). (3.71) 



In fact, by Lemma \TM L(T>i(X) nonres ) = L(x 2 ). By Lemma l2~T51 L(T>i(A)) = L{ X \). 
By Corollary EHZ1 L(xl) = L(X 2 xl), and > therefore, L(T>i(A)) = L(3C 2 Xi)- It follows: 
L(D 1 (A)\Di(A) nonres nonres 

) = L{X 2X \)-L{ X2 ) = L(X 2X 1\X2), 
the last equality being valid since X2 is a subset of X 2 x*- Using the estimate ()3.10|) in 
the Geometric Lemma, we obtain 

l(Dx(X) \ T>i(X) nonres } < L(T>i(X))k~ 2 ~ 2si . (3.72) 
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Now considering that L(Di(A)) = 2irk(l +o(l)) (Lemma 12. llj) . we arrive at the estimate 
(l3~TTI) . Combining it with (l3~7771) . we get ifSlg]) . 

(3) By the third statement of Lemma 12.111 the function >c\{ip) is holomorphic in <Pi and 
\W(xi((p)) = A. If <p £ 6 2 , then A( 2 ) (xi(<^)) is defined by the series (|3.55j) . since 
xi(ty?) £ I ) i(A) ri onres for any (/9 £ 02- Now we show that A^ 2 ) (xi (</?)) can be analytically 
extended from O2 to ^1 \ 0. It is enough to check that the series (|3,55|) converges in 

#i\0. We need an estimate for (xi (</?)) , which is similar to (|3.17j) . Indeed, considering 
formula (j3.13|) for and taking into account that HW2II < e-K we easily conclude that 



it suffices to prove the estimate: 

ff«(*i(¥>)) 
From the definition of ff", we see: 



sup 

2GC2 



1 

< — . 



(3.73) 







r 


= max 







HW(Z 1 (<p)+2Trp/N 1 a) - 

To estimate the right-hand side, we start with the case p = 0. We obtain an upper bound 
on || (j?i(<^)) — z) || using results of Section^ Indeed, %\><i{ip) belongs to Xi(^); 
which is in the complex (/c _1_4si_2<5 )-neighborhood of the first non-resonance set Xi(^)- 
Hence, we can apply Lemma lA. II from Appendix IA.4I with t = %ixi{tp). By this lemma, 

(#W(j<l(p)) - Z y l has a single pole zq inside the contour C\. It is easy to see that 
zq = A, since A^ 1 ) fxi (<£>)) = A. Again, by Lemma [A. 11 



(z-k 21 ) (HV>{*y{<p))-Z 



for all z inside C\. Considering z £ C2, that is, \z 
(I3~7H) that 



< 16 (3.74) 
k 2l \ = ei/2, we easily obtain from 



< 



zeC 2 . 



ei 



(3.75) 



Next, we consider the case p £ P \ {0}. We use Lemma 13.311 with b = 2irp/Nia, p 7^ 0. It 
is possible to replace k 21 in (|3.52j) by any z £ C 2 , since the difference between z and k 21 
is small: |* - k 2l \ = ei/2, e x = o (6 ^" 9 " 12si " 6 ' 5 ) • Thus, 



J3'W(x 1 (^) + 27rp/iVia) 



< 



e 1 



(3.76) 



for any p £ P \ {0}, z £ C 2 . Combining (|3.75|) with the last estimate, we arrive at 
(|3.73j) . which means that the series for A^ 2 ) (a, i?i(<^)) converges. It is easy to see that 
estimates (|3.75|) . (|3.76j) are stable in a vicinity of f<i((p), namely, they hold at all points 
kv: \x—K\{ip)\ < eik~ 2l+1 ~ s . Hence, the perturbation series for the eigenvalue converges 
for such x too, the following estimates being valid: 



VA (2) (a,xz7 



2lk 2l ~ 2 k + o(k 2l - L ). 



< 12aef , 



(3.77) 
(3.78) 
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Next, we want to show that the equation \( 2 \a, xz?) = A has a solution X2(</j) 9 for every 
ip G <5i \ 0. It easily follows from (g^jj) - that 

^M =2i ^ 1(1 + 0(1)) , W!) =0(fc -, ,3.79) 
ax ax z 

for any </? G and x : \x—k\ < /c~ 1 ~ 4si ~ 2<5 . Using these estimates we readily obtain that 
jA^^xi?) — A| > k 2l ~ 1 ef for all x on the circle |x — >ci(p)\ = ef. Considering (|3,77|) and 
applying Rouche's theorem, we obtain that there is a single solution >C2{<p) of the equation 
A^(q,xz7) = A such that {^(p) — xi(f>)\ ^ e i- Applying (|3.78|) and Implicit function 
theorem, we obtain that ^(f) is a holomorphic function of ip in \ and estimates 
(l3~7l3l hold. 

Let us make here a remark for the future. Convergence of the series for the resolvent 
(i?( 2 ) (#i(¥>)) -*)~ , 2 G C*2, following from ((HZZHJ), means that the resolvent has a single 
pole z = (Ki(ip)i') inside C<i- Similar result holds when we replace £i(<p) by j<2(<p), 
since >ii((p) and %2{<p) are close: |x2((£) — xi((/?)| = o(ei). Considering that A( 2 )(x2(</?)) = 

A, we obtain that (z - A) (F (2) (x 2 (</?)) — z) is holomorphic inside C2 and the estimate 
similar to (|3.74[) holds: 



- A) (tf (2) (i? 2 (v?)) - zj < 32. (3.80) 
(4) First, we show that 

L(D 2 (A)) = L(lMA) nonres )(l + 0(fc- 2 - 2si )). (3-81) 

A— >oo V / 

Second, we check that 

L{V x {X) nmres ) = L(B 1 (A))(l + 0(fc- 2 - 2si )). (3.82) 

A— >oo V / 

From these two estimates (|3.65(l follows. In fact, considering that 

L(D 2 (A)) = / J x 2 + x 2 2 dip, L('Di(A) nonres ) = / a/x 2 + x^ 2 dp, 

and taking into account that li2(tp) = X2(<£>) — xi(</?) satisfies (|3.63|) . we easily get (|3.81|) . 
Formula (|3.82j) easily follows from (|3.72j) . 

■ 

Now define the non-resonance set X2W i n ^2 (A) by 

X^(A):=0C 2 2) 2 (A). (3.83) 

Lemma 3.43. The set X2W belongs to the (2aefk~ 2l+l ^ -neighborhood o/%2(A) in K2. If r G 

(2) 

X2*(A), then the operator (r) has a simple eigenvalue equal to A. This eigenvalue is given by 
the perturbation series \3.15\) . where p G P,j G 1? are described as in Geometric Lemma \3. 51 
part 2. 



In fact, it should be X2(a, 'fi)', we omit a for shortness. 
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Proof. By Lemma \'AA2\ T>2(\) is in the (2ae\k~ 2l+l ) neighborhood of T>i(X) nonres . Considering 
that X2W = ^^(A) and X2(A) = %2^>i(X)nonres (see ()3.54j) ). we immediately obtain that 
X<!>(A) is in the (2aef fc _2i+1 )-neighborhood of X2(A). The size of this neighborhood is less than 

e\k~ 2l+1 ~ 5 , hence Theorem 3.1 holds in it, i.e., for any r G Xs>(A) there is a single eigenvalue 

(2) 

of H£> (t) in the interval £2(^5 <5). Since xI(A) C £2 (A), this eigenvalue is equal to A. By the 
theorem, the eigenvalue is given by the series (|3.15j) . where p G P,j G Z 2 are described as in 
Geometric Lemma 13.51 part 2. ■ 

Lemma 3.44. Formula h'J. 8°J\) establishes one-to-one correspondence between x^iX) and T>2(\). 

Remark 3.45. From geometric point of view this means that X2W does not have self-intersections. 

Proof. Suppose there is a pair x. >r G ©2(A) such that %2>c = %2>c* = t, r € X2(A)- By the 
definition (|3.58j) of X>2(A), we have \^(a, x) = \^ 2 \a, H*) = A, i.e., the eigenvalue A of H^\t) 
is not simple. This contradicts to the previous lemma. ■ 

3.5. Preparation for the Next Approximation. 

3.5.1. Contracted set s (b). Let us recall that 0(6) is the union of disks of the radius r = 
£,-4-6si-3<5 cen tered at the zeros of the unperturbed determinant det(/ + ^4o(yb( ( / 9 ))) > see Defini- 
tion EHU We proved that the "perturbed" determinant det(/ + Ai(y(ip)) has the same number 
of zeros inside each r(6) as the unperturbed one, when 60 is big enough. To prepare the next 
approximation, we contract the set T(b) around zeros of the perturbed determinant. First, we 
consider the disks with the radius r' = rk~ 2 ~ As ' 2 ~ s centered at each zero <pj : o of the perturbed de- 
terminant det(/ + Ai(y (</?)) in the set 0(b). Obviously each of these disks is in 0(b), since r' « r 
and the distance between zeros of the perturbed and unperturbed determinants is smaller than 
r/2 (Xemma l3.35l Part 2). Next, we take the union of these smaller disks and denote it by s (b), 
here the index s stands for "small". Clearly, s (b) C 0(6). We denote by T s (b) a connected 
component of these new disks which intersects with $>i and by 7.5(6) its boundary, 7^ = dT s . Let 
us remind that in all our previous considerations 6 satisfies the conditons 60 > fc~ 2i + 9 + 12s i+ 7(5 . 10 

Lemma 3.46. (1) The set s (b) contains the same number of disks as 0(b). This number 
does not exceed J = cok 2+2si . 
(2) For any if G £>i \ s (b), 

-1 



.21 



< c6 k 



1t,-2Z+9+12si+6<5 



(2rJ/r) 



(3.84) 



H^(y(f))-k 21 
Corollary 3.47. If b > k- 2l+9 + 12s ^+ 7S , 

ff (1) (#))-fc 2r ' 



< cb 1 k- 2l+u+Us ^ e5 (2rJ/r') J . 



< 



(3.85) 



(3i 



10 



Indeed, we had two restrictions on bo: bo > k 



7+8si -27Q+65 



CLemma ET^ and b > k- 2l+9+12si+7S (Lemmas 



Recalling that 70 = 21 - 2 - 4si - 28, we get 7 + 8sj - 270 + 63 < -21 + 9 + 12si + 75, we can combine 
two restrictions into b > k- 2t+9+12si+7S . 
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21 



-1 



l ei 



k 2+2si 



(3.87) 



Corollary \3.47\ We use the condition on rj: 77 > 2 + 64/ (21 — 11) = 2 + 2/ s\ and the formulae for 
ei = e~4 1 , r/r' = /c 2+4s2+<5 . Simple computation shows that the right-hand side in p.84|) does 
not exceed 1/ei. Hence, the corollary holds. ■ 

Lemma \3.4b\ By construction, the set S (6) contains the same number of disks as 0(6). Since 

the number of disks in 0(6) does not exceed J = cok 2+2si (see the remark after Definition 13.241) . 

the same is true for s (b). 

If cp 6 #1 \ 0(b), then the estimates just follow from Lemma 13.311 Thus, we have to consider 

' H^(y(<p)) -k 2l \ . , , 

is holomorphic m every connected 



0\0 S . Note that (H^ (y((p)) - k 2l )~ l det 



component T(b) of 0(6). We denote zeros of det 



H (y(cp)) + k 21 

( H^(y(<p))-k 21 



\ H (jf(<p)) + k* 



in T by ip jj0 , j = !,-■■ , J*, 



J* < J. Obviously, (HW(ff(<p)) - k 21 )- 1 Ujliiv 

- ifjfl) is holomorphic in T. For any <p on 7, 

J* J* 
[HU(y(<p)) - k 21 )- 1 H(<p - ^ )| < \\(HW(m) ~ k 21 )^ I] \<P ~ ml 
j=l j=i 

Considering that \<p — <Pj.o\ < 2rJ* for any (p € 7 and the estimate ()3.52j) . we get 



||(if«(y») " k 21 )- 1 II(^ - ^o)l! < cb^k- 2l ^ 2 ^\2rJ 



(3i 



By the maximum principle the inequality (|3.88|) holds for any tp E T. Note that any ip € T \ T s 
satisfies the inequality \<p — ipj t o\ > r' and, hence, 



H^(y(<p))-k 



21 



< cb 1 k- 2l+9+12s ^ 6S (2rJ*/r 



holds for ip £ V \ T s . Considering that r/r' > 1 and J* < J, we obtain (|3.84j) . To prove ()3.85|) . 
we use the following Hilbert relation: 



(H^(m) 



k 



2i 



H (y(cp))+k 21 ) +{H (y(v))+k 



,21 



-1 



W 1 + 2k 2l )(HM(y(ip))-k 



21 



-1 



(3.89) 



Summarizing the terms of a diagonal operator, we easily get 



H (y(tp)) + k 



21 



< ck 



-2i+2+2si 



(3.90) 



Substituting ()3.84|) and the last estimate into (|3.89|) . we obtain H3.85|) . 
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3.5.2. The set s (b) for small bo- Everything we considered so far is valid for b obeying the 
inequality bo > fc— 2J+9+12si+7<5^ ^ere 5 ^ s foe distance from b to the nearest vertex of K\. We 
showed in Section f3. 3. 51 that b = 2irp/Nia, p € P\{0}, satisfies the requirement. However, in the 
next section and later, bo will be taken smaller, since the reciprocal lattice is getting finer with 
each step. To prepare for this, let us consider b being close to a vertex of K\\ 

< b < k- 2l+9+12s ' +75 . (3.91) 

We show that for such b the resolvent (H^{y[(p)) — k 21 ) 1 has no more than two poles ip^ in 
<Pi and even its small neighborhood. We surround these poles by two contours ^ and obtain 
estimates for ^H^\y(ip)) — k 21 ) 1 when ip is outside 7 1 * 1 . 

In fact, it easily follows from (j3.91|) that: b = o(fc~ 1 ~ 4s ' 1-2 ' 5 ) as k —> 00. Let i>\ be the 
(^A; _2_4si_2<5 ) neighborhood of <&\. Suppose |6| = bo, i.e., the closest vertex of K\ for & is (0,0). 
The perturbation series (|2.34j) converge for both A^ (v?)) and (y((p)) when tp € and 
both functions are holomorphic in #1, because and y(<p) are in the complex (2/c _1_4si_25 )- 

neighborhood of Do(A) nom . es for such ip (Lemma f2.11|) . Note that AW (i?i((^)) = k 21 for all cp € i>\. 
We base our further considerations on these perturbation series expansions. For b being close to 
a vertex e other than (0, 0), we take y(<p) = + b — e . 

We define <pf, € [0, 2ir) by the formula b = 6o(cos <pb, sin^b) when \b\ = bo, and by the analogous 
formula b — e = bo(cos <pb, sin ipb) when b is close to a vertex e other than (0,0). 

Lemma 3.48. Ifb satisfies \3. 91)) and |cq | < 6o^ 2 '~ 3 ~ 4si ~ 3<5 ; then the equation 

\W(y(ip)) = k 2l + e (3.92) 
has no more than two solutions, ipf^, in <P\. They satisfy the inequality 

\vf -(<p b ±n/2)\<h- 2 -^- 2S . (3.93) 

Proof. Suppose W\ = and \b\ = bo, i.e., the closest vertex of K\ for b is (0, 0). Then the equation 
(|3.92j) has the form \kv + b\ 21 = k 21 + eo- It is easy to show that it has two solutions <pf satisfying 
(|3.93|) . Applying perturbative arguments and Rouche's theorem, we prove the lemma for nonzero 
W\. A detailed proof is in Appendix IA.5I In the case when b is close to a vertex other than (0, 0), 
the considerations are the same up to a parallel shift. ■ 

Lemma 3.49. Suppose b satisfies h'J. 91\) and p £ <3>\ obeys the inequality analogous to h'J. 9J\) : 

\ip - {<p b ± vr/2)| < k- 2 - 4si ~ 25 . Then, 

^A« (y( v )) = fc ^oo ±2lbok 21 - 1 (1 + o(l)) . (3.94) 

Proof. Let W\ = and |6| = b . Then X^(y(cp)) = \kP+b\ 21 and 

\kv + b\ 21 = 2l\kv + b\ 2l ~ 2 (kv + b, kp), fx = (- sin tp, cos ip). 



dip 

For ip close to ipi, ± ir/2, we have (&,/?) = ±&o(l + o(l)). Considering also that (ft, u) = 0, we 
obtain -§^\kv + b\ 21 = ±2lbok 2l ~ 1 (1 + o(l)). Applying perturbative arguments, we get a similar 
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formula for nonzero W\. For a detailed proof see Appendix IA. 61 In the case when b is close to a 
vertex other than (0,0), the considerations are the same up to a parallel shift. | 

Definition 3.50. Let Tf(b) be open disks centered at tp^r € 3i 11 with the radius r' introduced 
in the previous subsection, r' = r A;" 2- 4s2 " 5 ; jf(b) be their boundary circles and s (b) = TfuTj. 



Lemma 3.51. For any tp in <P\ \ O s (b), 

\\W(y^))-k 2l \>b k^- s r\ 



(3.95) 



Proof. Suppose (|3.95|) does not hold for some tp in &i\Q s (b). This means that tp satisfies equation 
(EHI21) with some e : |e | < bok^r' < b k 2l - 3 - 4si ~ 35 . By Lemma I3~1H1 tp obeys (ETMj) . Thus tp 
could be either tpf Q or tp~ Q . Without loss of generocity, assume tp = tpf Q . Applying (|3.94|l in the 
segment between tp and tp$ , we obtain: eq = X^\y(tp)) — k 21 = ±2lbok 2l ~ 1 (l + o(l)) [tp — tp^. It 
immediately follows that \tp — tp^\ < r' , i.e., tp G Tf C O s (b), which contradicts the assumption 

¥>e$i\0 s (6). ■ 



Lemma 3.52. For any tp G $1 \ O s (b), 



H^(y{tp))-k 



21 



< 



16 



bQr l k 2l-l-8 ■ 



(3.96) 



H^(y(tp))-k 



2! 



< 



16cqA; 



2+2si 



1 bor'k 21 - 1 - 5 ' 

Corollary 3.53. If e x k- 2l+l - 28 < b < fc-2/+9+i2s 1+ 75 and (p€ 1 \ o a (&), then 



HW(&(<p))-k 



2i 



< 



:2 - 



(3.97) 



(3.98) 



H^(y(tp)) - k 



21 



-1 



< 



co 



k 2+2si 



(3.99) 



Corollary follows from the condition on 60 and the estimate r'k 38 > 16ei, which is obviously 
valid for sufficiently large k. 

Proof. Since y(tp) is in the (2/c- 1 - 4si " 2 ' 5 )-neighborhood of D (A) nonres-, we can apply Lemma 
IA. II from Appendix IA.4I By this lemma, 



\^{y{tp))-k 2l )(H^{y{tp))-k 



,2/ 



-1 



< 16. 



Using (|3~95|) . we easily get (|3~96|) . Using (|3~96|) and considering (|3~89|) . (|3~90|) . we obtain (|3~97|) . 



n,„± 



^0 = ff j when £ o = 0. 
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4. The Third Approximation 

4.1. The Operator . Choosing S3 = 2s2, we define the third operator by the formula: 

#(3) = H m + aW3j (0 < a < 1), W 3 = ^, 

r=M 2 +l 

where M 3 is chosen in such a way that 2 Ms « P 3 . Obviously, the periods of W 3 are 2 M3 " 1 (6i,0) 
and 2 M3 - 1 (0,6 2 )- We write them in the form: N 2 N 1 (a 1 ,0) and N 2 N 1 (0,a 2 ), here N 2 = 2 Ms ~ M2 , 
\k S3 ~ S2 < N 2 < Ak S3 ~ S2 . Note that 

M 3 M 3 

IIW3IU < 2 ||K|U< J] exp(-2^) <exp(-fc^). 

r=Af 2 +l r=M 2 +l 

4.2. Multiple Periods of W 2 (x). The operator, = Hi+W 2 (x), has the periods N\ax,Nia 2 . 
The corresponding family of operators, {H( 2 \t)} t <zk 2 > acts in L 2 (Q 2 ), where Q 2 = [0, Aqai] x 
[0, JVia2] and .K2 = [0, 2ir/N\ai) x [0, 2n /N\a 2 ). Since now on we denote quasimomentum t 
from the first step by quasimomentum r from the second step by t^ 2 '. Correspondingly, 
quasimomentum for Ha^ we denote by . Eigenvalues of (i*- 2 ^) are denoted by A^ {t^), 
n G N and its spectrum by A^ 2 ^ (t^). 

Next, let us consider W 2 (x) as a periodic function with the periods N 2 N±ai, N 2 N\a 2 . When 
changing the periods, the family of operators {i^ 2 - 1 (i^) } t (2) g ^ 2 is replaced by the family of 
operators {H® {t i3) )} ti3)&K3 , acting in L 2 (Q 3 ), where Q 3 = [0, JV2JV1O1] x [0,N 2 N ia2 ] and K 3 = 

[0,2Tr/N 2 N iai ) x [0,2ir/N 2 N ia2 ). We denote eigenvalues of (t®) by a| 2) (t^), n G N, and 
its spectrum by A^ 2 - 1 (i^)- We denote now by P^ the set P, introduced by (|3.4|) . its elements 
being pW. By Bloch theory (see e.g.|21|l. for any G X 3 , 

A< 2 > (V 3 >) = |J A( 2 ) (t® + 2tt P ( 2 ) /N 2 Niaj , (4.1) 

p(2) e p(2) 

where = <|V 2 ) = (p?\j4 2) ) G Z 2 : < pf ] < N 2 - 1, < p { 2 2) < N 2 - l}, 2-7rp( 2 '>/N 2 N 1 a = 

\2np^ /a\, 2irp^ /a^ /N 2 N\. An isoenergetic set £2(^0) C K 3 of the operator is defined 
by the formula: 

S 2 (\) = {t® e K 3 : 3n G N : \®(t®) = a} 

= {t^ eK3 :3neN, p^ G P^ : A^ (V 3 ) + 2irpW /N^a) = a} . 

Obviously, S 2 = X 3 S 2 , where % 3 is the parallel shift into K 3 , that is, % 3 : R 2 -> if 3 , 3C 3 (t( 3 ) + 2irm/N 2 Nia) = 

m G Z 2 , t^ 3 ) G if 3 . We denote index j, introduced in Part 1 of Geometric Lemma l2.11 by 
and j, introduced in (JSHSJ), by j( 2 \ = ;W +pW /N x . 
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4.3. Perturbation Formulae. The results of this section are analogous to those in the previous 
step, the index 2 just being replaced by 3. Let us start with establishing a lower bound for k. 
Since ns\ > 2 + 2s±, there is a number k* > e such that 

C(l + Sl )k 2+2si lnk < ^ 1 ,C = 400/(c + l) 2 , (4.2) 

for any k > k*. Assume also, that /c* is sufficiently large to ensure validity of all estimates in the 
first two steps for any k > k*. 

Lemma 4.1 (Geometric Lemma). For a sufficiently large A, A > k 21 , there exists a non-resonance 
set X3(A,(5) C K3X2 such that: 

(1) For any point t^ G X3> the following conditions hold: 

(a) There exists a unique p^ G P^ such that t^ + 2-Kp^ /N 2 Nia £ x\- 

(b) The following relation holds: 



A 



(2) 

;(2) 



(V 3) +2vrp (2) /N 2 N ia 



k 



21 



where A^ (i^ + 2np^ /N 2 Niaj is given by the perturbation series k&.15\) with 



a. 



and j( 2 ) = j + p/Ni, here j and p are defined by the point r = + 2irp^ /N 2 N\a 
as it is described in Part 2 of the Geometric Lemma for the previous step. 
(c) The eigenvalue \ { %{t^ + 2yrp( 2 ) /N 2 N 1 a) is a simple eigenvalue of (t^) and 

its distance to all other eigenvalues \ { n ] (t (3) + 2irp^ /N 2 Nia) of (i (3) ) is greater 



. l-h r l s 2 

than €2 = e ± 



A (V 3 ) + 2npW/N 2 N ia ) - Af (V 3 ) + 2^/N 2 N l 



> e 2 . 



(4.3) 



(2) For any t^ in the fak 21+1 & ) -neighborhood inC 2 ofx3, there exists a unique p^ € P^ 
such that + 2vrp( 2 ) /N 2 Ni a is in the (£2^ 21+1 5 ) -neighborhood in C 2 of x\ an( ^ 



(t^ +2KpW/N 2 N ia )-k 21 



< e 2 k~ 



(4.4) 



i( 2 ) 



j + p/N\, here j and p are defined by the point r = t^ + 2-Kp^ 1 /A^iVia as it 



is 



described in Part 2 of the Geometric Lemma for the previous step. 
(3) The third nonresonance set X3 has an asymptotically full measure on x\ * n the following 
sense: 

T i<V _ , ,* \ ,,_ "\\ 

si-2s 2 



L(X 3X *2\X3)) _ k . 4 



(4.5) 



The proof of the lemma is analogous to that for Geometric Lemma in the second step. It will 
be presented in Section l4~4l 

Corollary 4.2. If t^ belongs to the complex {e 2 k~ 2l+l ~ 5 ) — neighborhood of X3{^,°~), then for 
any z lying on the circle C3 = {z : \z — k 2l \ = €2/2}, the following inequalities hold: 

(i/< 2 >(i< 3 >) ' ' 



< 



\\(H( 2 Ht®)-z)-%< 



£2 

4c fc 2+2 *3 

C2 



(4.6) 
(4.7) 
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Remark 4.3. Every point 2irq/N2Nia (q G 1?) of the dual lattice for periods A^A^ai, N<2,N\a<2, 
can be uniquely represented in the form: 2 r Kq/N2N\a = 2imi / 'Nia+2irp^ /N2N1C1, where m G Z 2 , 
p( 2 ) £ p( 2 ). Note that 2imi/N\a is a point of a dual lattice for periods N\a\, N\a2 and 

p (2) G P (2) 

is responsible for refining the lattice. By Remark 13. 7\ 2Trq/N2N\a also can be uniquely represented 
as 2irq/N 2 N ia = 2itj/a + 2vrp( 1 VA r ia + 2itp^ /N 2 N ia , here j G Z 2 , pW G pW, p^ G P (2) . 

Let us consider a normalized eigenfunction ip n (i( 2 ),x) of .ff( 2 ) (i^ 2 ') in £2^2)- We extend it 
quasiperiodically to ^2(^3), renormalize and denote the new function by ^ n (t (3) , a?) , t (3) = X 3 t (2) . 
The Fourier representations of VVi i ^) m £2(^2) and ^(t^, x) in L2{Qs) are simply related. 
If we denote Fourier coefficients of V>n (t (2) ,x) with respect to the basis \Q 2 \-^/ 2 e i { 2 - Km l N ^ a+t(2) ' x ) ) 
m € Z 2 , in I2 (Q2) by C nm , then, obviously, the Fourier coefficients C nq of ?/> n (i( 3 ) , x) with respect 
to the basis \Q 3 \-i/2 e i(^q/N 2 N ia +t( 3 \x) ^ q £ ^ in £ 2 (q 3 ) are given by the formula 

C nm , if q = mN 2 +p (2) ; 
0, otherwise, 

p( 2 ) being defined from the relation = t^ + 2irp( 2 ^ /N1N2CL. Correspondingly, matrices of the 
projections on i/) n (r,x) and tp n (t^,x) with respect to the above bases are simply related: 

' (EnUm, if q = mN 2 + p (2) , q = mN 2 + p {2) ; 
0, otherwise, 



{E n )qq 

E n and J5„ being projections in ^2(^3) and £,2(0:2), respectively. 



We define functions gi 3 \k, t^) and operator-valued functions Gr\k,t^), r = 1,2, ••• , as 
follows: 

#(M (3) ) = ^Tr £ ((iT^ 3 )) -zY'wsJdz, 

Theorem 4.4. Suppose belongs to the (e2k~ 2l+1 ~ & )-neighborhood in K% of the third non- 
resonance set Xz{\$)> < 76 < 21 — 11 — 16si. Then, for sufficiently large A, A > kf and 
for all a, < a < 1, there exists a unique eigenvalue of the operator Ha^ (^ 3 ^) in the interval 
£3$) : = (k 21 — €2/2, k 21 + £2/2), 62 = e~^ kVS2 . It is given by the series: 

CO 

\%(a,t^) = X% (#) + 2^/N 2 N ia ) + £^(3)^(3)^ (4.8) 

r=l 

converging absolutely in the disk \a\ < 1, where := + p^ / N2N1 , p^ , being described 
in Geometric Lemma \4-l\ The spectral projection, corresponding to A^(a, i^), is given by the 
series: 

00 

E%(a,t^)=E% (t^ + 2^/N2N ia ) + ^ r G^(k,t^), (4.9) 

r=l 

which converges in the trace class Si uniformly with respect to a in the disk \a\ < 1. 
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The following estimates hold for coefficients g£ {k,t^), Gy (k,t^): 



fl£ 3) (M< 8) ) <^H) r , G( 3 )(M (3) ) <6r(4 e 3 )' 



(4.10) 



Corollary 4.5. The following estimates hold for the perturbed eigenvalue and its spectral projec- 
tion: 



X%(a,t^)-X%(t^+2^/N 2 N ll 



,tW)-E%(t®+27rpW/N 2 N ia 



(2). 



< Y2ae\, 



< 48aei. 



(4.11) 
(4.12) 



The series ()4.8fl , (|4.9(l can be extended as holomorphic functions of t^ in the complex ( — e 2 k 2l+1 ~ 

neighborhood of X3] they can be differentiated any number of times with respect to t^ and retain 
their asymptotic character. The results analogous to Lemma 13.111 - Lemma 13 . 1 41 hold . 

4.4. Proof of the Geometric Lemma. We define the third nonresonance set, X3(A) C X3X2W 
as follows: X3 (A) = ^(A) \ Q 2 (A), 

Q 2 (\) = |t( 3 ) e X 3X * 2 : 3n, h G N, p (2) ,p (2) G P (2) , p (2) ^ p (2) : 

A^ (i< 3 > + 27rp( 2 )/iV 2 iV 1 a) = A, t< 3 > + 2np^ /N 2 N ia G X *(A), 



A^ (t® + 2vrp( 2 )/^ 2 iV ia ) _ A f (V 3 ) + 2^ /N 2 N ia ) | < e 2 }- (4.13) 



Proofs of Statements 1 and 2 of Geometric Lemma are identical to those for the second step up 
to the shift of indices by 1. 

From now on we denote a vector b G K%, introduced in the second step, by b^\ Naturally, we 
denote the set 0(6) constructed in the previous step ( Section EOH) by (1) (6 (1) ). Correspondingly, 
T(b) is rW^ 1 )), 0* = 0l 1} (see (JSH3)), and s (b) is 0^\b^). Note that O^^feW) is defined for 
an arbitrarily small see Sections 13 . 5 . 1 1 and 13 . 5 . 21 Let us recall that the complex non-resonant 
set 02 is given by (|3.4()|) . 

Definition 4.6. By analogy with the definition of Wi(A) (Section I3.3.4J) . we consider a complex 
resonance set a;|(A), which is the set of ip G ( P 2 satisfying 



det 



'HW@ 2 \<p)) -k 21 - 

H Q {^){if))+m 



0, ^\ v >) = M^)+b {2) , 6 (2) 



2tt P ( 2 ) 
N 2 Nia 



(4.14) 



for some p^ G P*- 2 ) \ {0} and |e| < e 2 . Let uj 2 (X) be the set of cp G ©2 corresponding to £l 2 {\). 
Considering as in the proof of Lemma I3.2U1 it is easy to show that u 2 = u\ n ©2- 

By analogy with (|3.41j) , let us consider an arbitrary b^ G K 2 and its distance b^f 1 to the closest 
vertex of K 2 : 



l(2) 



mm 

m=(0,0),(0,l),(l,0),(l,l) 



b {2) - 2-Km/Nii 
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Now we construct 0^ (b^ 2 ^ C <P 2 , assuming > 0. Indeed, for each b^ + 2np^ /N\a, 
p(l) G p(l) 

, we take the set 0^ (6^ 2 ^ + 2np^ /N\a) described in Section \'6. 51 and consider a union: 

(2) (gW) = U p(1)6P(1) 0W (&( 2 > + 2^ 1 )/iV ia ) . (4.15) 

A connected component of 0^ 2 ^ f&^l we denote by T^ 2 ^ (M^J an d its boundary by 7^ / ' 

Remark 4.7. The set O^ 2 ) (^^ 2 ' ) ) consists of disks of the radius r' around all zeros of the deter- 
minant det^/ + A\ (y^ 1 )^))^ , here and below: 

& x \ip) = M<P) + b (2) , (4-16) 

This remark easily follows from the fact that each 0^ (b^ + 2np^ /N\a\ is built around zeros 
of det^I + A\ (^xi(ip) + fr( 2 ) + 2ttp^ /Nia^j, and A\ {if l \ifj) is a direct sum of the operators 
Al (x« + 6( 2 ) + 2vrp( 1 )/Af 1 a), pW G P« 

Remark 4.8. Let us estimate the number of disks in 0^ 2 ^ ^^ 2 ^ an d the srze OI a connected 

component T( 2 ) (&P0). Indeed, the set 0^ (b^ contains no more than Ac k 2+2s2 disks, since 

each 0^ consists of no more than cok 2+2si disks (Lemma l3.461 Definition 13 ,5U|) . and contains 
no more than 4fc 2 ( S2_ ' 51 ) elements. The set O^ 2 ) [b^ ) is formed by disks of the radius r', here 



r' = rk 2 4s2 s , r being the radius of disks constituting 0* , r = k 4 6si 3S . Hence, the size 



of each connected component of O^ 2 ) does not exceed 8cor'k 2+2s2 = 8cork~~ s ~ 2s2 . Thus, 

the size of I^ 2 ) (j)( 2 ^ is much smaller than the radius of disks constituting ol 1 ^. Here, even if a 
component I^ 2 ) (S^^j is not strictly inside <P 2 , it has the same properties as a component inside 
<5 2 - Further, we consider r( 2 ) (b^ C <P 2 , (2) (& (2) ) C <2> 2 - 
By analogy with ()3.4U[) . we introduce new notations: 

0l 2) = U p{2)€Pi2) \ {0} O {2) {27Tp^/N 2 N 1 a), $ 3 = $ 2 \ ■ (4.17) 
We will show ( Corollary that uj* C ol 2) and uj 2 C ol 2) fl 2 . From now on we call 0i 2) the 

(9} 

second (complex) resonance set in <& 2 and 0* H 2 the second resonance set in Q 2 . From (|4.17f) 
and (|4. 15|) we see: 

0l 2) = U p( x) 6 p(x) iP (2) 6P (2) UQ} 0« (2TT P V/N ia + 27rpW/N 2 N ia ), (4.18) 

each Og 1 '' (^irp^ /N\a + 2np^ /N 2 Nia^ being a union of disks with centers at zeros of the deter- 
minants det (i + Ai + 27rpW /Nia + 2irp^ /N 2 N 1 aj ) , the size of disks being r'. 



SPECTRAL PROPERTIES OF A POLYHARMONIC OPERATOR WITH LIMIT-PERIODIC POTENTIAL 45 




Figure 8. The set S . 



Remark 4.9. The set 0^ consists of no more than 16coA; 2+2s3 disks of the radius r', since each 
0^ is formed by no more than c k 2+2si disks and pW, contains no more than 4k ( S2 S1 ' 
and 4k 2 ( S3 ~ S2 ' ) elements, respectively. 

Obviously, 0* contains more disks than 0^ , 



(2) 

however, disks in 0; are of a smaller size than 



ones in 



Disks in are centered at zeros of unperturbed determinants, while smaller disks 



in are centered about zeros of the perturbed determinant (see Remark 14.7(1 . The shift of 
centers is of principle importance: we can reduce the size of disks, since they are more precisely 

(2) 

targeted. If W\ = 0, then O* is just a union of disks centered at quasi-intersections of the 
circle of radius k, centered at the origin, with circles of the same radius k, centered at points 
2irj/a + 2irpW /N x a + 2-npW /N 2 N ia , p^ 7^ 0, of the dual lattice for periods N 2 Nxax, N 2 N\a 2 . 

To obtain 0^ we delete 0* from 2 . Thus, the set ^3 has a structure of Swiss cheese (Fig. |HJ); 
we will add more holes of a smaller size at each step of approximation. 

By analogy with 1)3.31(1 . we define the operator A2 (y^ (</?)) by the formula: 

I + A 2 (y< 2 H^) = (H^(y< 2 \(p)) - A) (HoG^foO) +A)~\ 

(ip) being given in Definition [OJ H Q (y^ 2 ) (</>)) acting in L 2 (Q2)- The goal is to show that 
det(l + Ii(yW(y?))) and detfl + A 2 {^ 2 \ip))\ have the same number of zeros inside each 
r( 2 )(5(2)) anc j that the resolvents for HW (yW (if)) and -ff( 2 ) (y( 2 ) (</?)) satisfy similar estimates 
when tp G 0i \Q( 2 )(b( 2 )). 

Lemma 4.10. Ifb^ G K 2 and > eik- 2l+1 - 2S 
(1) For any (p G X \ 0^ (b^ 

'[^(y^i^-k 



then 



21 



2/ 



-1 



< 



4cnk 2+2s * 



(4.19) 
(4.20) 



(2) The number of zeros of det yl + A\ (y^ 1 ) (92) 



/// 



does noi exceed 4cok 2+2s2 . 



16 
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Corollary 4.11. There are no zeros o/det( I + Jli (^(ip)) j in <Pi \ (2) (^ (2) J • 
Proof. (1) Considering the relation between and H^ 1 ' , we easily see 



(i^^))-*; 2 ' 



max 

p (i) e p(i) 



■■21 



The condition 99 G £>i \ 0( 2 ) (&( 2 )) yields tp G #1 \ oi 1} + 27rp( 1 )/iV 1 a) for every 

p(i) G p(i). if 5(2) + 2 7 rpW/N 1 a satisfies 6 > A;- 2 '+ 9 + 12s i+ 7<5 , where b is the distance 
from + 2np^ /N\a to the nearest vertex of ifi, then we use Lemma l3.46l and Corollary 



Suppose 6^ 2 ) + 2irp^ /Nia satisfies b < k~ 



H {1) (y {1) (p) + 27rp^/N 1 c 



;2i 



1 1 

< — < -T7 



(4.21) 



2«+9+12si+7<5 



. It is easy to see that 60 > 



eik~ 2l+1 ~ 2S , since b^ > e\k~ 2l+l ~ 2S by the hypothesis of the lemma. Using Corollary 
13.531 to obtain an estimate similar to (|4.21() . Thus, ()4.19|) is proved. The estimate ()4.2U() 
can be checked in the analogous way from (|3.87|) and ()3.99j) . when we take into account 
that P^ contains less than 4fc 2s 2 -2s i e i emen t S- 
(2) The determinant of / + A\ (^^(c/?)) is equal to zero if and only if at least one of the 

determinants det(l + A 1 (y^(tp) + 2irp^/N 1 a)^j , G is zero. Since each of them 

has no more than cok 2+2si zeros (Corollary I3.36( Definition I3.50J) and P^- 1 ' contains no 

more than 4k 2 ( S2 ~ Sl ^ elements, the total number of zeros of det(/ + ii (yW (</>))) does 

not exceed 4cq/c 2+2s2 . 



The following five lemmas are the analogs of Lemmas 13.291 - 13.381 
Lemma 4.12. Ifb^ G K 2 and b^ > e^- 2 ^ 1 - 26 , then: 

(1) The number of zeros of det(^I + A\ {if 2 ^ ((/?)) ^ inside (b^^j * s the same as that of 
det(l + ii(tfWfa))). 

(2) For any <p G #2 \ 0^ (b( 2 ^ , the following inequalities hold: 

2 



< 



-2 ' 



H^\y< 2 \p))-k 



21 



-1 



< 



8cofc 2+2 



*2 



Corollary 4.13. There are no zeros of detyl + A\ {if 2 \<-p))j in <P 2 \ 0^ [b^ 
Proof. (1) Direct calculation gives: 



detjl + ii^ 



det(7 + Ii(yW)) = det(/ + T). 



(4.22) 
(4.23) 
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-1 



T=(W 1 -2X) (h {v< 2) ) + A (Ho&V) - H (yM) - A 



-i 



The estimate ()3.63f) yields: \if- 2 \ip) — y"(y>)| < 2aefk 21+1 . Using the last estimate we 
readily show that 



-i 



< 2laefk- 21 . 



(4.24) 



Considering (|4.2U|) . we obtain 

< 4A: 2/ (Ale\k~ 21 ^ (Ac k 2+2s2 /ej ) < 64lc ejk 2+2s2 < 1. (4.25) 

Hence, |det(J + T) — 1| < 1 for any (p E 7^ ■ Applying Rouche's Theorem, we 

complete the proof. 
(2) Let us consider the Hilbert relation: 

(P(/ 2 )(,)) - a)" 1 = - a) _1 + t,(P(#(,)) - a) -1 , 

To estimate ||T*||, we represent T* in the form T* = (I + ?i)T2, here 
Tx = (#« (y« fa)) - A) _1 (-Wi + 2A) , 

Using pn^|> . we get || < 4AeJ~ 2 . The norm of T 2 can be bounded as (|4.24|) . Multiplying 
the estimates for ||Ti|[, 1 1 1 1 , we obtain: ||T*|| < 8ale 2 < 1/2. Substituting the last 
estimate into the Hilbert relation and considering (|4.19|) . we arrive at (|4.22|) . Using 
(|4.2U[) . we prove ()4.23|) in the analogous way. 



Lemma 4.14. If € K 2 and b^ > e%k 21+1 25 , then the following estimates hold for any 

<pe$ 2 \ (2) (& (2)N 



H®(y< 2 Hv))-\ 
H( 2 \y< 2 \ V )) - X 



4 

< ~2 ' 



-1 



< 



16cq/c 



2+2s 2 



Corollary 4.15. There are no zeros of det( I + A 2 (y^ 2 \(p)) J in<P 2 \ ®^ (b^ ) . 



(4.26) 
(4.27) 



Proof. The proof is analogous to that of Lemma f.-i.^H we use Hilbert relation, the estimates (|4.22|) . 
Q4.23|) and the inequality \\W 2 \\ < ef, which immediately follows from (|3.2|) and the definition of 
ei- ■ 

Lemma 4.16. Ifb^ G K 2 and b^ > eik~ 2l+1 ~ 25 , then the number of zeros o/det (/ + A 2 (y^ 2 \(p))) 
inside r( 2 )(6 (2) ) is the same as that of det ( / + An 
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Proof. Using (|4.2H|) . we obtain that for any cp € 

H (2) (j^)) -fc»' 



det 



HW (y( 2 %)) — k 21 



det I 



det I + — 



W 2 



< 2IIW2I 



■21 



< 2ej 



Hi 1 ) (y(2)(yj)) -k 21 

A h 2+2s 2 



det I 



8c eik z 



< 1. 



By Rouche's Theorem, we complete the proof. | 
Lemma 4.17. (1) 7/6 (2) G K 2 andb^ > eik~ 2l+l ~ 2& , then the number of zeros of det (I + A 2 (y (2) (<p))) 



in any ri 2 \b^) is the same as that of det (l + ii^ 1 ) (</?)) 
(2) The distance from a zero of det (l + A 2 (y^ 2 \p))) to the nearest zero of 
det (j + does not exceed A 2 )/2. 

Corollary 4.18. The number of zeros of det (/ + A 2 (y( 2 ) (ip))) does not exceed 4cok 2+2s2 . 

The corollary is the combination of the first part of the lemma with the second part of Lemma 

mm 

Lemma \4-lb] The first statement immediately follows from Lemmas 14. 1 21 ( Part 1) and 14.161 Next, 
note that Lemmas l4.1UI - H.16l hold not only for 0( 2 )(&( 2 )), but also for a set 6 (2) (6 (2) ), constructed 
from disks of twice smaller radius A 2 ) /2 with the same centers, since all estimates in the lem- 
mas preserved under such change up to some insignificant constants. This proves that all zeros 
of det (1 + A 2 {y {2) {if))) in 0( 2 )(&( 2 )) are, in fact, in the smaller set 6 (2) (6 (2) ), i.e., the second 
statement of the lemma holds. ■ 

Lemma 4.19. Ifb^ G K 2 and bf ] > e^ 21 ^ 26 , \e\ < ef, then 

H&) (y( 2 ) Up)) - k 21 - e /- \ 

(1) The determinant of . . — — — — has no zeros in # 2 \Q(2) 6(2) . 

Ho(ifl 2 > (tp)) -\- k V / 

f HW(y< 2 Xp))-k [ 



,21 



(2) The number of zeros of det 
/HW(y< 2 Xp))-k 2lS 



V H (yi 2 )(p)) + k 21 



in 



r( 2 ) 6( 2 ) is the same as that 



of det 



H (yi 2 )(p))+k 2 ^ 



(2) 



Corollary 4.20. The following relation hold: uj^ C 0* 

Lemma \4-l^\ and Corollary \4-2b\ Proof of the lemma is completely analogous to that of the previ- 
ous two lemmas up to replacement of the estimate HW2II < e i by ||W2+e|| < 2ef. Corollary follows 
from the lemma and Definition 14.61 In fact, uj 2 (X) = U p (2) e p(2)\{ }W 2 p {2) , where u 2p (2) is the set 

of zeros of the determinant (|4.14|) for a fixed p^ 2 \ By Lemma 14.191 (Part 1), uj 2tP (2) C 

-, since e 2 < ef. Note that such &( 2 ) satisfies the condition of the lemma, since 



N 2 Nia\ :2 < 2k S3 . Using the definition (|4.17|) of 0* and the formula ^(A) = U p (2) g p(2)\j- } w 2 p( 2 ) 
yields the corollary. ■ 
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Prom now on we denote the radius r of the disks in ol 1 ^ by and the radius r' of the disks 
in 0(2) byr( 2 ). 

Lemma 4.21. (proof of Statement 3) Under conditions of Geometric lemma \%~l\ , estimate \4-5\ ) 
holds. 

Proof By Corollary l4~2Ul C oi 2) . Hence, uj 2 C ol 2) n 9 2 , the set cu 2 being given in Definition 
14,61 Considering that that ol 2 ^ is formed by no more than 16cofc 2+2s3 disks of the radius r^ 2 ) (see 
Remark 14. we obtain that the total length of lo 2 does not exceed 32cofc 2+2s3 r^ 2 \ Taking into 
account that L(0 2 ) « kL{u 2 ) and r^ = r (l)fc-2-4 S2 -.5 = fc -6-6 Sl -4s 2 -45 ? S3 = 2g2 = 4s ^ we 

obtain that the the length of f2 2 does not exceed 64co/c~ 3_2si ~ 2,S2 ~ 4<5 . Considering that x 2 W nas 
a length 27r£;(l + o(l)) we arrive at Statement 3. | 

(3) 

4.5. Nonresonant part of the isoenergetic set of H a • This section is analogous to Section 
13.41 for the second step. Indeed, let 53(A) be an isoenergetic set of the operator H a 3 ^: 53(A) = 
{#) G K 3 : 3n G N : \ { n ] (a,t^) = A}, here {A^ 3) (a, t (3) )}^° =1 is the spectrum of H a 3 \t^). Now 
we construct a non-resonance subset X3 (A) of S3 (A). It corresponds to non-resonance eigenvalues 
A .(3j(t( 3 )) given by the perturbation series l|4.8j) . We start with a definition of D 2 (\) nonres . Recall 
that X3 C 3^3X2(A) (see Geometric Lemma l4.1j) and x 2 (A) = 3C 2 I> 2 (A), see (|3.83|) . Hence, 
X3 C 3C 3 D 2 (A). Let T> 2 (\) nonres be the preimage of X3 in £2 (A): 

D 2 (A) nonres — 

{x G ©2(A) : JC^x E X3}- (4.28) 
The following lemma is analogous to Lemma 13.401 in the second step. 

Lemma 4.22. The formula 3C 3 D 2 (A) nonres = X3 establishes one-to-one correspondence between 
©2 (A) nonres 

and X3- 

Proof. The proof is analogous to that of Lemma r3.4Ul up to the shift of indices by 1, i.e., X2 — ► X3> 
X*(A) -> x|(A), r = t (2) -> t (3) ; we use formula (|3.83|) instead of (|2.48jl . Part la of the Geometric 
Lemma for the third step instead of Part la of the Geometric Lemma for the second step, and 
Lemma E3U instead of ■ 

We define 23(A) as the set of directions corresponding to T> 2 (X) nonres : 

B 3 (A) = {v G 3 2 (A) : x 2 (A, v)v G V 2 (\) nonres } 

where x 2 (A, v) is defined by Lemma 13.421 x 2 (A,z?) = x 2 {ip). Note that '33(A) is a unit circle 
with holes centered at the origin and '33(A) C 32(A). We denote by 63 (A) the set of angles ip 
corresponding to S3 (A): 

6 3 (A) = {<p G [0,2tt) : (cos 93, sin v?) G B 3 (A)}, G 3 C 9 2 . (4.29) 

We define 1>3(A) as a level set for A^ 3 ^(a, x) in the (e 2 A; _2 ' +1-<5 )-neighborhood of D 2 (\) nonres : 

D 3 (A) := [x = xv: G 3 3 (A), |x - x 2 (A, z7) | < e 2 A;" 2m ^, A (3) (a, i?) = a} . 

Next two lemmas are to prove that 2) 3 (A) is a distorted circle with holes. Their formulations and 
proofs are analogous to those of Lemmas 13.411 and 13.421 
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Lemma 4.23. For every V E ^(A) and every a, < a < 1, i/iere is a unique x = x 3 (A,z?) in 
i/ie interval I 3 := [x 2 (A, z?) - e 2 fe~ a+1 ~' 5 , x 2 (A, z?) + e 2 A;~ 2Z+1 ~ 5 ] suc/t i/iai 

A (3) (a,x 3 z7) = A. (4.30) 

Furi/iermore, 

|x 3 (A,z7) - x 2 (A,z?)| < 2aetAT 2m (4.31) 

Let us recall that ^ 3 is the third complex non-resonant set given by Q4.17JI . see Fig. |HJ Further 
we use the notations x 3 (</?) = x 3 (A,z?), h$(<p) = x 3 (y>) — x 2 (<£>), x 3 (</2) = *t${ip)v. 

Lemma 4.24. (1) The set ©3(A) is a distorted circle with holes: it can be described by the 
formula: 

£ 3 (A) = {xeR 2 :x=x 3 (^), c/?ee 3 (A)}, (4.32) 

where >c 3 {ip) = x 2 (<^) + h 3 ((p), x 2 (y>) is the "radius" of D 2 (A) and h 3 ((p) satisfies the 
estimates 



< 2aeik' 



-21+1 



dh* 



< 4aelk 1+S . (4.33) 



dip 

(2) The total length of !B 3 (A) satisfies the estimate: 

L (B 2 \ S 3 ) < 4vr£r 4 ~ 2si - 2s2 . (4.34) 

(3) The function x 3 {(p) can be extended as a holomorphic function of <p to the complex non- 
resonce set ^3, estimates {4-3$ being preserved. 

(4) The curve 2) 3 (A) has a length which is asymptotically close to that o/D 2 (A) in the following 

sense: 

l(d 3 (A)) a = qo L(© 2 (A))(i + 0(A;- 4 - 2si - 2s2 )). (4.35) 

Proof. The proof is analogous to that of Lemma 13.421 Note only that, in Part 2, when proving 
convergence of the series for the resolvent (i?^ 3 ^(x 2 (y?)) — z) , we use the estimate 

sup (H&ifaW-zY 1 <—, (4.36) 
zeC 3 v ' £ 2 

analogous to (|3.73|) . the operator acting in L 2 (Q 3 ). The estimate (|4.36|) follows from (|4.26|) 
for 6 = 2np^/N 2 N l a, p^ ^ 0, and estimate (|3.8U|) . which takes care about the case p' ' = 0. 
As a side result of these considerations, we obtain an estimate analogous to (|3.8U|) for the new 
resolvent (ff( 3 )(x 3 (ip)) - z) . ■ 

We define the non-resonance set, xI(A) in S3 (A) by the formula analogous to (|3.83|) : 

X5(A):=0C 3 D 3 (A). (4.37) 
The following lemmas are analogous to Lemmas 13.431 and 13.441 

Lemma 4.25. The set X 3 (A) belongs to the (2ae2k~ 2l+1 ) -neighborhood of xsW in K 3 . Ifv® G 

X 3 (A), then the operator Ha\t^) has a simple eigenvalue equal to A. This eigenvalue is given 
by the perturbation series j4.&\) . 

Lemma 4.26. Formula establishes one-to-one correspondence between X 3 (A) and D 3 (A). 

Remark 4.27. From geometric point of view this means that X 3 (A) does not have self-intersections. 
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4.6.1. Contracted set 0^ \ S^J. For any such that > e\k 21+1 2S , we have constructed 

the set O^ 2 )^ 2 )) C <p2, which surrounds zeros of the determinant det^/+Ai (^^(v?))! ; V^if) = 
+ b^ 2 \ see Remark 14.71 The set consists of disks of the radius A 2 ^ = r ( 1 )^-2- 4s 2-<5 ) r (i) _ 
^.-4-6si-3<5 anc j con tains no more than 4cok 2+2s2 disks, see Remark 14.81 We proved that the 
new determinant det^I + ^(y^Hv 7 )))) j7^(v) = ^(v 3 ) + °*( 2 \ has the same number of zeros 
inside each T^ 2 \b^) as det^I + .^(y^ 1 ) (</?)) ^ (Lemma 15.111 Part 1) and both have no zeros 
in G <P 2 \ 0( 2 )(6( 2 )) (Corollaries l4"TTl l4~m To prepare the next approximation, we contract 
the set di 2) (6( 2 )) around the zeros of the new determinant det ^/ + A2 (y^ 2 ^ ((/?)) ^ . We consider 
the disks with the radius = r ( 2 )^-2-4.s3-<5 centered at each zero (fjQ of the determinant 
det^I + A2 (y^ (<p)) J ■ Obviously each of these disks is in 0( 2 )(&( 2 )), since r^ « A 2 ^ and the 
distance between zeros of the old and new determinants is smaller than A 2 " 1 /2 (Lemma l5.111 Part 
2). Next, we take the union of these smaller disks and denote it by oi 2 ^(6^ 2 -*). We denote by 
(&) a connected component of these new disks, which intersects with and by 7! (&^) its 
boundary, js = dT^ ■ 

Lemma 4.28. (1) The set 0^ (b^) contains the same number of disks as 0( 2 )(&( 2 )). This 
number does not exceed 4cofc 2+2s2 . 
(2) For any <pe<P 2 \ (2) (b^), 



(y< 2 H^) - k 21 



< 



V 2 )j( 2 )' 

r (3) 



J(2) 



j( 2 ) = 4cn k 2+2s \ 



(4.38) 



-1 



< 



4J (2) / 2r (2)j(2) 
r(3) 



J(2) 



(4.39) 



The proof is analogous to that of Lemma 13.461 Lemma 14.141 being used instead of Lemma 
13.311 Considering the expressions 

for ,.(2)^(3)^(2) 

, €2 , we easily obtain that the right-hand part 
of (|4.38|) does not exceed l/e2- Hence, the following corollary holds: 

(2) 



Corollary 4.29. I/bjf > eik~ 2l+1 - 25 , then 



2/ 



< 



e-2 



(4.40) 



< 



4cok 2+2s2 



£2 



(4.41) 
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4.6.2. The set 0^ (b^A for small b^. Everything we considered so far is valid if b^ > t\ k 21+1 2S . 



(2) 

However, in the next chapter and later, 6q is taken smaller, since the reciprocal lattice is get- 
ting finer with each step. To prepare for this, let us consider b^ being close to a vertex of K2: 
< b { Q ] < eik- 2l+1 - 2S . We show that for such b^, the resolvent (h® (y^((p)) - k 21 ^ has 
no more than two poles in and even its small neighborhood. We surround these poles by two 
contours 7^ ^ and prove an estimate for the norm — k 2l ^j when ip is outside 

these contours. 

Let <£>2 be the ^^--neighborhood of <I>2- It is easy to see that §2 C&x- Suppose |&( 2 )| = b^\ i.e., 
the closest vertex of K2 for b^ is (0,0). We consider the functions A^f^ 1 ) (</>)) and X^ (if- 2 ^ (p)) 
defined by perturbation series (|2.34|) and (|3.55|) for p G <?>2- The convergence of these series 
can be easily justified. In fact, by Lemmas 12.111 and 13.421 H\{ip) and oki^P) are holomorphic 
functions of (p in ^2- It is easy to show that the lemmas hold even in ^2, since all estimates 
involved are preserved in ^2- The perturbation series ()2.34|) and (|3.55|) converge for X^\x x {p>)) 
and A ( 2 )(x2 (</?)), respectively, because >c\{p) and ^(p) are m small neighborhoods of T>o(X) nonres 
and Hi (A) 

nonres (see (2.75), (3.75)). Since the estimates involved are stable with respect to a 
change of xi^ not exceeding eik~ 2l+1 ~ 5 , the perturbation series for A^ (y^ ((,?)) and A^ 2 ) [if 2 ^ (<p)) 
also converge, both functions being holomorphic in ^>2- We base our further considerations on 
these perturbation series expansions. For 6( 2 ) being close to a vertex e other than (0,0), we take 
y( 2 \tp) = x 2 M+6 (2) -e. 

From now on, we denote the solutions pf of the equation A^ (y^ (y>)) = A; 2 ' + £o> introduced 
in Lemma 13.481 by pt ^ ■ As in the previous subsection, we choose = r ( 2 ) k~ 2 ~ iS3 ~ s here. 

Lemma 4.30. If < b^ < ei£r 2m ~ 2<5 and |e | < 6q 2) fc 2i-1 "-M 3 ), tfien i/ie equation 

X®(f 2 \p))=k 2l + e (4.42) 

/ias no more £/ian two solutions cpt ^ in For any tpf ^ there is p^ 1 ^ € such that 

<r( 3 )/4, (4.43) 



±(2) ±(1) 



here and below p^ 1 ^ is pt} 1 ^ for eo = 0. 

Proof. Let \b^\ = b^ ] . First, we expand X^ (y^ 2 \p)) and AW(yW(<p)) near the point b^ = 
and consider that X^{5<2 (</?)) = A^ (i?i (</?)) = A. Then, using l|3.57j) . the estimates of the type 
C22J, HOB and (pTHOj) . we check that 

A (2) (y< 2) (¥>)) - AW(^(^))| < 4 2) ei (4-44) 

in <^2 and even in its /-^-neighborhood, the neighborhood being a subset of S±. Suppose (|4.42|l 
holds for some p € & 2 . By &m . X^(y^{ip)) = k 21 + e' , e' = O (^fc^-M 3 )) . Hence, p> 

satisfies conditions of Lemmas 13.481 and 13.491 Surrounding p> by a circle C of the radius r( 3 )/4, 
we see that {X^ (^(p)) - k 21 - e' \ = \lk % ~ x bf\^ (1 + o(l)) > \e' Q \ on this circle. Applying 
Rouche's theorem, we obtain that there is a solution of X^ l \if l \ip)) = k 21 inside this circle. Thus, 
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any solution of Q4.42JI is in the circle of the radius /4 surrounding <p^^\ the point ip^^ being in 
the (r( 3 )/4)-neighborhood of <p2- It remains to check that Q4.42|) has no more than two solutions 
and no more than one in a vicinity of each (Pq^ ■ We construct the disk of the radius 

centered at (p^^, described above, and note that (<p)) - k 2l \ = i//c 2 ' _1 &Q 2 M 3 ) (1 + o(l)) 

on the circle. Using Q4.44|) and Rouche's Theorem, we obtain that there is only one solution of 
Q4.42|l in the disk. If b is close to a vertex other than (0,0), the considerations are the same up 
to a parallel shift. | 

Lemma 4.31. Suppose < < eik~ 2l+1 ~ 2S and <p £ <p2 obeys the inequality analogous to 



±(1) 
<p - vv 



< r^ . Then 




dtp' 



A (2) (y< 2 ^)) =^oo ±2lb®0-\l + o(l)). 



(4.45) 



The proof is completely analogous to that of Lemma 13. 491 

Definition 4.32. Let T^ 2 -* (b^j be the open disks centered at (p^^ with radius r®; ^yf^ (b^) 
be their boundary circles and oi 2) (b^) = T^ (2) U rj {2) . 

Lemma 4.33. For any ip in 

|A( 2 )(y( 2 )( 95 )) - k 2l \ > b^k 21 - 1 - 5 ^. (4.46) 
The proof is analogous to that of Lemma 13.511 
Lemma 4.34. For any ip <E ^2 \ 0^ (b 1 - 2 ^ , 



< 



16 



■■21 



b (2) r (3) fc 2J-l-i' 

64c fc 2+2s2 
1 " bPrWkF- 1 - 6 ' 



< 



Corollary 4.35. If e 2 k 



-21+1-28 ^ t(2) 



< C < erk- 21 ^- 25 , and <p£<P 2 \ 0? ] (& (2) ), then 



,.2/ 



-1 



(4.47) 
(4.48) 



(4.49) 



(tf (2) (y^M)-^)" 1 



< 



4c £; 2+2s2 



(4.50) 



The corollary follows from the condition on b^\ for the formula A 3 ^ = r^k 2 4f53 5 = 
£.-8-6si-4s 2 -4s3-5<5 anc [ estimate (|4.2|) . The proof of the lemma is analogous to that of Lemma 
13321 
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5. The u-th Step of Approximation. Swiss Cheese Method. 

5.1. Perturbation Formulae. On the n-th step, n > 4, we choose s n = 2s n _i define the 
operator by the formula: 

M n 

H (n) = H (n-l) + aWnt ( <a<l), W n = V ^ 

r=M n - i+l 

where M n is choosen in such a way that 2 M ™ w k Sn . Obviously, the periods of W n are 2 Mn ~ 1 (bx, 0) 
and 2 M ™~ 1 (0,6 2 )- We write them in the form: iV n _i • ... • Ni(ai,0) and iV n _i • ... • iVi(0,a 2 ), here 
-/V"n_i = 2 M -~ M -\ J***— < A T n _ 1 < 4fc^-'»-i. Note that H^JU < Er=A/ n _ 1+ i II^IU < 
exp(-/c^- 1 ). 

The n-th step is analogous to the third step up to replacement the indices 3 by n, 2 by n — 1, 
the product A^iVi by iV n _i • ... • iVi, etc. 

We note that k, satisfying 1)4.2(1 . obeys the analogous condition for the n-th step: 

(7(1 + s n - 2 )k 2+2Sn - 2 lnk< k^ Sn - 2 (5.1) 

with the same constant C. The inequality (|5.1|) can be obtained from (|4.2|) by induction. This is 
an important fact: it means that the lower bound for k does not grow with n, i.e., all steps hold 
uniformly in k for k > k*, /c* being introduced by (|4.2|) . Further we assume k > k*. 

The formulation of the geometric lemma for n-th step is the same as that for Step 3 up to a 
shift of indices, we skip it here for shortness. Note only that in the lemma we use the set Xn-iW 
to define Xn(A). In fact, we started with the definition of Xl(ty an d then use it to define xf(A) 
(Step 1). Considering Xi(A), we constructed X2(A) fLemma l3.5() and later used it to define X2W 
( Section I3.4|) . Using X2W1 we introduced X3(A) (Lemma 14. 1|) and, then X3W (Section 14. 5|) . 
Thus, the process goes like xi ~ * Xi ~~ * Xi ~ > X2 ~^ X3 ~ > X%- Here we consider the set Xn-iW 
defined by (|4.37|) for n = 4 and by (|5.18|) for n > 4. 12 The estimate (|4.5|) for n-th step takes the 
form: 

L(XnXU\Xn)) <fc ^ Sn = 2j2(l + s t ). (5.2) 
L Un-lJ i=1 
It is easy to see that S n = 2(n — 1) + (2™ — 2) si and S n 2 n s\ i=a s n for large n. The formulation 
of the main results (perturbation formulae) for n-th step is the same as for the third step up 
to the shift of indices. The formula for the resonance set £l n -i and non-resonance set \ n are 
analogous to those for f2 2 , Xi ( see <|4.13|) ). 

5.2. Proof of Geometric Lemma. Proofs of the first and second statements of Geometric 
Lemma are simple and completely analogous to those in the second step and third step. We skip 
them here for shortness. Proof of the third statement is completely analogous to that for the 
third step. For the purpose of rigor we present a short version of it. 



12 Strictly speaking we assume that there is a subset Xn-iW of the isoenergetic surface S„-i(X) of iJ^™ -1 ) 
such that perturbation series of the type (14.811 . 14. 91 converges for i*-™ -1 ' £ Xn-i(-M an d Xn-i(^) has properties 
described in Section f4.5l up to replacement of 3 by n — 1. In particular, we assume that Xn-i(^) — % n -i'Dn-i(X), 
where D n _i(A) satisfies the analog of Lemma 14.241 and that the analogs of Lemmas 14. 251 and 14.261 hold too. Here, 
%n-i is the parallel shift into K n -\. Further in this section we describe the next set XnW which has analogous 
properties. 
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In the second and third steps we defined the sets Qto (P-A and O^ 2 ) (b^ ( see Definition 

EHjand formula Considering (33B> , we define O^" 1 ) f^"" 1 )), fc^ -1 ) G n > 4, 

by the analogous formula 

(»-i) (gC*- 1 )) = U p( n_ a ) 6 p(„- a) 0( n - 2 ) (fc^ 1 ) + 2vrp("- 2 )/iV n _ 2 a) , (5.3) 

here and below iV„_ 2 = N n _ 2 • ■■■■ ' N±, 

P (m) = {p (m) = ( p M p M) , o < pS m) < N m - 1,0 < < 7V m - 1}; 

the set oi"^(6( m )), m > 1, is a collection of disks of the radius r( m+1 ) = r M&- 2 - 4s m+i-< 5 around 
zeros of the determinant det^I + A m (jf™) (<p))j in <£ m , here y( m \ip) = % m ((p) + b^ m \ x m {(p) 
being defined by Lemmas l2~TTl Cm = l),E3U(m = 2), gH (m = 3) and EHTl Cm > 4). Let 

= U p( „_ 1)eP( „_ 1)U0} O("- 1 )(27r^"- 1 )/iV ri _ 1 a), <P n = ^ n _! \ Ol^ 1 '. (5.4) 

We denote by r^™ -1 ) ^b^ 1 ^ 1 ^ a connected component of O^ 1 ) ^6( n_1 )^ . Note that the complex 
non-resonance set is defined by the recurrent formula analogous to (|4.17|) . 

Lemma 5.1. The set oi"^ (b( m ^ , 6^ m ) G K m contains no more than A m ~ 1 c k 2+2Sm disks. 
Corollary 5.2. TTie sei 0( n_1 ) (b^ 1 ^ contains no more than 4 n ~ 2 c k 2+2Sn - 1 disks. 

Corollary 5.3. The set 0* ^ contains no more than 4 n ~ 1 c /c 2+2Sn disks. 

The lemma is proved by induction procedure. For m = 2 it holds by Part 1 of Lemma 14.281 
and by Definition 14.321 For m > 2 it follows from Corollary 15.211 and Definition 15.261 Rigorously 
speaking, we assume for now that the lemma holds for m < n — 2. At the end of this section, we 
obtain the analog of the lemma for the next step (n — > n+ 1), this will complete the induction 
procedure. Corollaries 15.21 and 15.31 are based on the fact that contains no more than 

4^.2(s„-s„_i) e l emen ts and a similar estimate holds for p( n-2 ). 

Obviously, <I> n has a structure of Swiss cheese, more and more holes of smaller and smaller radii 
appear at each step of approximation (Fig. [SJ. Note that the disks are more and more precisely 

"targeted" at each step of approximation. At the n-th step the disks of 0* are centered around 
the zeros of the determinants 

det(/ + A„_ 2 (x„_ 2 (¥>) + 27rp("- 2 ViV n _ 2 a + 27rp(^ 1 )/iV„_ 1 a)) , 

where p( n-2 ) g pi n ~ 2 ) ; j/™- 1 ) g p( n_1 ) f the corresponding operator H^ n ~ 2 ^ being closer and closer 
to the operator H at each step of approximation. If W\ = Wi = ... = W„_2 = 0, then ol™ ^ is 
just a union of circles centered at quasi-intersections of the circle of radius k, centered at the origin, 
with circles of the same radius centered at points 2nj/a + 2np^ / 'N±a + .... + 27rp( n ~ 1 )/iV n _ia, 
which are points of the dual lattice corresponding to the periods iV n _iai, iV n _ia 2 . 

The following is a proof of third statement of Geometric Lemma for n-th step, namely, estimate 
(|5.2[) . It is completely analogous to that in the third step, we only have to provide that condition 
(j5.1j) is sufficient for all necessary estimates. We start with an analog of Lemma 14.101 we shift 
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indices I —> n — 2, 2 — > n — 1 and replace the coefficient 4co in the estimate (|4.2()[) by 4 n ~ 2 co- The 
operator fl"( n ~ 2 ) is i/( n ~ 2 ) extended from Q n ~2 to Q n -\- By analogy with ()4.16j) . if n ~ 2 \if) = 
x n _ 2 (c£>) + b^ n ~ l \ the function being defined in ^ n _ 2 . Thus, we have the following lemma and a 
corollary: 

Lemma 5.4. IfP^V G K n ^ and b { ™~ 1] > e n -2^ 2l ^ 2i \ then 
(1) For any if G n _ 2 \ ( " _1) (& (n ~ 1} ) : 

^- (5.5) 



f("- 2 )(^- 2 )(^))-^) 1 

jy(n-2) (j/(»- 2 )( ¥ ,))-jfea)" 1 



< 



c n-2 



4 n - 2 c /c 2+2s "- 1 



c n-2 



(5.6) 



(2) The number of zeros of detyI+A n ^2{y( n 2 \ip))j in <P n ^ 2 does not exceed 4 n 2 cok 2+2Sn - 1 . 

Corollary 5.5. There are no zeros o/det^J + A n _ 2 (y( n ~ 2 )(92))^ in 
&n-2 \ (n_1) (b {n ~ l) ^} . 

Proof. Proof is analogous to that of Lemma I4.1UI We use Corollary 14.291 for n = 4 and Corollary 
15.231 for n > 4, as well as Corollary 14.351 (n = 4) and Corollary 15.291 (n > 4). 13 Part 2 of the 

lemma follows from the definition of o( n_1 ) ( frl™" 1 ) ) and Corollary 15.21 | 



Next, we repeat Lemmas 14.121 - 14.211 up to replacement of indices. Proofs are completely 
analogous to those for Step 3. In Lemma 14.211 we show that the total length of 0„_i does not 
exceed k~ Sn , here J7 n _i is the analog of fii, Q 2 , see (j3~5|) . (JUTSJ). Let ]f n ~ x \<f>) = i? n _i(^)+6 (n " 1) , 
the function being defined in ^ n -\. 

Lemma 5.6. IfU n ~^ G K n -i and b^ 1] > e n ^k' 2l+1 ~ 2S , then: 

(1) The number of zeros o/det^J + A n -.2(y^ n ~ 1 \ ( p))^j in <iny 

r (n-l) pn-i)^ ja ^ e 

same as 

tfiai /det(/ + i n _ 2 (y< n - 2 )(^))). 

(2) For any (p G <£ n -i \ O^ -1 ) (6^ n_1 M ; i/ie following inequalities hold: 



Hin-2),^n-l)^_ k 2l 
H(n-2)^n-l)^_ k 2l 



< 



< 



fc n-2 

2 • 4 n - 2 c /c 2+2s ' 1 - 1 



e 2 

e n-2 



(5.7) 
(5.8) 



Corollary 5.7. There are no zeros of det + A n _ 2 (if n ^((p) 



13 



Rigorously speaking, we have to assume that the lemma holds for n-th step. For (n + l)-th step it follows 
from Corollaries 15 . 231 and 15 . 29l 
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Proof. The proof is analogous to that of Lemma 14.121 In Part 1 we use the estimate jy^™ -1 -* (<p) — 
y( n ~ 2 \<p)\ < 2ae^_ 2 ^ ~ 2l+1 , which follows from assumed properties of Xn-l> see l|4.33|l with 3 — > 
n — 1. We apply (|5,6|) instead of (|4.2U|) and produce an analog of Q4.25JI . the estimate 16 • 
4 n ~ 2 coe 2 l _ 2 k 2+2Sn - 1 < 1 following from (|5.1|) . In Part 2 we use the obvious estimate e^„ 2 < 
1/4. " " | 

Lemma 5.8. Ifb^ 1 ^ G Kn-i and &q > e n _2k~ 2l+1 ~ 2S , then the following estimates hold for 
any ip G <2>„_i \ O^ 1 ) (b^A : 

Hl^Ujfi^vj) - &Y 1 (5.9) 



ra-2 

(*<-•> < irvfpn. (5 . 10 , 

Corollary 5.9. There are no zeros o/det (^1 + A n _\ (yt™- -1 ) (v))J ^ n 

*fi-l \ (n_1) (6 {n_1) ) • 

Lemma 5.10. Ifb^ 1 ^ G i^n-i and b^" 1 ^ > e n -2k~ 2l+1 ~ 2S , then the number of zeros ofdet(l + 
4^1 (j/*"- 1 ^))) in fl ™2/ r(n_1) (& (n-1) ) is fte same as i/iai o/ det(/ + i n _ 2 (y (n " 1) (^) 



Lemma 5.11. (1) If b^ n ^ G K n -\ and ^ > e n _2^ 2 ' +1 25 > iaen the number of ze- 
ros of det(l + ^ n _i(y< n-1 )(^))J in any r^ 1 )^" -1 )) is t/ie same as f/ia£ o/detfl + 

i n - 2 (^- 2) M)). 

(2) The distance from a zero of det^/ + A n -\i^ n ~^ {<p)yj to f/ie nearest zero o/det^/ + 
^4n-2 (y^ n_2 ' ) (<£>)) ) does no£ exceed r( n-1 )/2. 

Corollary 5.12. T/ie number of zeros of det (J + Ai-i (^""^(v 9 ))) does not exceed ^ l ~ 2 CQk 2+2Sn - x . 

The corollary is the combination of the first part of the lemma with the second part of Lemma 
15.41 The following is an analog of Definition 14.61 

Definition 5.13. A complex resonance set w*_ 1 (A) is the set of ip G &n-i satisfying 

for some b {n ~^ = 2tt P ^-^ /N n -ia, p (n " 1} G P {n - l) \ {0}, and |e| < e n - X . We denote b y Un- i(X) 
the set of cp G n -i corresponding to f2 n _i(A). Considering as in the proof of Lemma T3. 201 it is 
easy to show that uJ n -\ = n Q n -i- 

Lemma 5.14. Ifb^^ G if n _i and b ( Q~ 1] > e n ^ 2l ^ 2 \ |e| < e*_ 2 , i/ien 
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-k 2l -e 

(1) The determinant of y—, — -r r — has no zeros in 

H (y("-i)(<p))+k 21 



(H^-V {#)) -k 2l -e\ 

(2) The number of zeros of det j— ■ — . ' — in any 

\ H (yy n L )(ip)) + k M J 

r (n-l) 5(n-i) is the same as that of det .J 1W \^ J „, . 

Corollary 5.15. The following relation holds: C oi" . 

Lemma 5.16. The estimate h5.ty) holds. 

Proof. By Corollary EM C oi n_1) . Hence, w n _i C ol n_1) nG n _i. By Corollary Ol ol n ~ 1} 
is formed by no more than 4 n_1 coA; 2+2Sn disks of the radius r*™ -1 ) = r ( n - 2 )k~ 2 ~ 4:Sn - 1 ~ 5 . Hence, 
the total length of u) n -\ does not exceed 2 • A n ~ 1 cor ( - n ~ 2 ^k~ s . Considering the recurrent formula 
f or r (m+X)^ r (m+l) = r (m)^-2-4s m+ i-<5^ we eas Qy get ^hat 4 n c r (n ~ 2 ) < /c" 5 ", n > 1. Therefore, 

the length of fi n _i does not exceed l/c -5 ' 1 . Considering that Xn-i(A) has a length 27r/c(l + o(l)) , 
we obtain (|5.2|) . ■ 

(n) 

5.3. Nonresonant part of the isoenergetic set of H a • Now we construct a nonresonance 
subset x^(A) of the isoenergetic surface S n (X) of H a n ^ in K n , S n (X) C K n . It corresponds to 
nonresonance eigenvalues given by perturbation series. The sets X*(A), xjK^); X;s(A) are defined 
in the previous steps as well as the non-resonance sets Xi(A), X2(A), X3(A)- Let us recall that 
we started with the definition of Xi(A) and then use it to define x*(A)- Considering Xi(A), we 
constructed X2(A) (Step 2). Next, we defined X5>(A), using X2(A). Using xl(^)> we introduced 
X3(A) and, then xI( A )- Thus, the process looks like xi X* X2 -> X2 X3 -» Xl- The 
geometric lemma in this section gives us X4 an d every next Xn if Xn-l i s defined. To ensure the 
reccurent procedure we show now how to define Xn(A) using Xn(A). 

As in previous sections, we start with the definition of D n -i(X) nonres . First note that Xn-l i s 
defined by the formula x* n -i = ^n-i^ n-ii where D n _i is a distorted circle with holes, and the 
shift gives one-to-one correspondence between Xn-l an d ^n-i (it is proven for n = 2,3,4 and 
we we assume this here for n > 4 to provide an induction procedure). By analogy with previous 
steps, we define D n _i(A) nonres as the preimage of Xn in Dn-i(A): 

Dn-l(A) nonres — {xGD n _x(A) :X n xGXn}, (5.12) 

Dn-i(A) nonres being defined since Xn C Xn— i- We define !B n as the set of directions corresponding 
to D n —\ (A) nojir . eiS . 

S n (A) = {z? € Si : x n _i(A, z?)z? G ©n-i(A)„ „ res }, 

where x n _i(A, z7) is the "radius" of the distorted circle D n _i(A) defined by the analog of Lemma 
14.241 for step n — 1. Note that S n is a unit circle with holes centered at the origin and 23 n (A) C 
1>n-i(A). We denote by G n the set of angles (f in polar coordinates, corresponding to r B n : 

9 n = {(/? G [0, 2vr) : (cos <p, sin tp) G @ n C 9 n _i. (5.13) 
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We define D n (A) as a level set for X^ n '(a, x) in this neighborhood: 

D n (A) = {x = xz7: z7e t B n ,\x->c n _ l (\,i?)\ < e n -ik~ 2l+1 - S , \V)(a,2) = a} . 

Considering as in the previous step, we prove the analogs of Lemmas 14. 231 and !4.24l For shortness, 
we provide here only the second lemma. By analogy with previous sections, we shorten notations 
here: x n ^i{p) = x n _i(A,z?), x n _i(v?) = x n _i(A,z7)z7. 

Lemma 5.17. (1) The set D n (A) is a distorted circle with holes: it can be described by the 
formula: 

B n (A) = {xeK 2 :i = ^), ^e n (A)}, (5.14) 
where >c n {ip) = x n ^\{ip) + h n (cp), and h n (p) satisfies the estimates 

< Aael^k 1 ^. (5.15) 

(2) The total length ofB n (X) satisfies the estimate: 

\B n ) <47rAT s ". (5.16) 

(3) Function >e n ((p) can be extended as a holomorphic function of <p to the complex non- 
resonance set <P n , estimates h5.15\) being preserved. 

(4) The curve T) n (X) has a length which is asymptotically close to that of D n _i(A) in the 
following sense: 

L(D n (A)) A = oo L(D„_ 1 (A))(l + 0(A : - 5 ")). (5.17) 

Now define the nonresonance set, x«(A) in S n (X) by the formula analogous to (|4.37l) . Indeed, 

X * n (X):=X n V n (X). (5.18) 
The following lemmas are analogous to Lemmas 14.251 and 14.261 

Lemma 5.18. The set Xn(A) belongs to the (2ae^_ 1 /c _2 ' +1 ) -neighborhood of Xn(A) in K n . If 
f( n ) £ Xn(A), then the operator Ha '(i^) has a simple eigenvalue equal to X. This eigenvalue is 
given by the perturbation series analogous to \4-8j) - 

Lemma 5.19. Formula i5.18\) establishes one-to-one correspondence between xJj(A) and D n (X). 
5.4. Preparation for the Next Approximation. 

5.4.1. Contracted set di™ ^ (b^ 1 ^) . We constructed the set Q( n_1 ) (6(" -1 )) , which surrounds 
zeros of the determinants 

det(/ + A n -2 (x n _ 2 (^) + & (n " 1} + 27r^™- 2 ViV„_ 2 a)) , p^ G P {n ~ 2 \ 

or, shorter, zeros of the determinant det^/+>l n _2 (^"^(v 3 )) Ji here y^ n ~ 2 \ip) = x n _2( ( /?)+& ( ' n_1 \ 
A n _2 corresponds to if( n ~ 2 ). The set consists of disks of the radius r'^ 1 ) and contains no more 
than 4 71-2 co/c 2+2Sn - 1 disks. We proved that the new determinant detfl + A„_i(y( n_1 )((^))^, 
y^ n ~ l \ip) = i? n _i(<^)+fr( n_1 ), has the same number of zeros as the old one inside each r( n_1 ) (ft( n-1 )) 



< 2aef l _ 1 k 



-21+1 



dh r , 



dip 



(iO 
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(Lemma 15.11(1 : and both have no zeros in € ^ n -i \ (& ^) (Corollaries 15 .5( 15.9(1 , To pre- 

pare the next approximation, we contract the set oi™ 1 - ) (^ 71-1 )) around the zeros of the new 
determinant. We consider the disks with the radius = r i n ~ l )k^ 2 ^ iSn ^ & centered at zeros of 
the new determinant. Obviously each of these disks is in 

(n-l) (p(n-l)^ 

since « r^ n ^ 

and the distance between zeros of the old and new determinants is smaller than A n ~^ /2 (Lemma 
15.11(1 . Next, we take the union of these smaller disks and denote it by 0^™ ^ (6^ n_1 ^). Obviously, 
the following lemma holds. 

Lemma 5.20. The set Qs ^' (b^ n ~ 1 ^ belongs to CK" -1 ) (fr^™ -1 )) and consists of the same number 
of disks as Q( n - 1 )(6( n - 1 )). 

Corollary 5.21. The set oi n ~^ (&( n_1 )) contains no more than 4 n ~ 2 c$k 2+2Sn ~ 1 disks. 

The corollary follows from the lemma and Corollary 15.21 
Lemma 5.22. For any ip G <2> n _i \ oi n-:L) (fr^" 1 )) , 

J(n-l) 



< 



-n-2 



r (n) 



:2i 



-i 



< 



4j(n-l) / 2 r (n-l) j(n-l) 



j(n-l) 



~n-2 



j(n-l) = 4 n-2 CQ/ 



-O) 



e * 



(5.19) 
(5.20) 



The proof is analogous to that of Lemma 13.461 Using the recurrent formula for A n \ condition 
(|5.1|) and a simple induction procedure, we easily obtain that the right-hand part of (|5.19|) does 
not exceed l/e n _i. Hence, the following is true: 

Corollary 5.23. 



< 



1 



tfCn-D(^-l) (¥0) 



.2/ 



j(n-l) 



< 

1 £n-l 



5.4.2. The set { P 1] (b^- 1 ^) for small b^' 1 ^. Everything we considered so far is valid if 6q > 



(5.21) 
(5.22) 

(n-l) 



e «-2^ 2i+1 25 - However, at each step 6g™ 1 ^ is taken smaller, since the reciprocal lattice is 
getting finer. For this reason we now consider ft( n_1 ) being close to a vertex of if n _i: < 



,(n-l) 



< e n _ 2 A:- 2 '+ 1 - 25 . As in previous steps, we show that for such ft( n ^ the resolvent 



(y("- 1 )( V p)) - fc 2Z 



has no more than two poles in <? n _i and even its small neigh- 
borhood. We surround these poles by two contours 7 S and obtain an estimate for the norm 



of (^H^ n 1 ) (if n ^((f)) — k 2l ^j when cp is outside these contours. 



14 



14 If |g(n-i)| = fo ("- 1 > ; th e n ■jf-"-^( lf ) = H„-i(ip) + b^-V. For b 1 ™" 1 ' being close to a vertex e of K n -i other 
than (0,0), we take y {n ~ x \<p) = J? n _! (<p) + b^" 1 ' - e. 
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Jn-2) 



By analogy with Step 3, ^ ra _i is the — ^ — neighborhood of ^ n -l- Obviously, <P n -i C @ n -2- We 
assume that the analog of Lemma 14,301 holds when we replace the index (2) by (n — 2). Namely, 
we assume that for any |eo| < b^ ^pZ-i-^O-i) the equation \( n ~ 2 \if n ~ 2 \(p)) = k 21 + eo has 
no more than two solutions </?i^ n 2 ^ in $> n -i and the inequality 1^^™ ^ — ^o^™ < A n ~^/A 

holds, here <p^ n ^ is <pt for eo = 0. We also assume that the estimate analogous to (|4,45j) 
holds when we replace (2) by (n — 2). Let us formulate the corresponding results for the next step 
of approximation (n — 2 — > n — 1). The formulations and proofs of the following four lemmas and 
are analogous to those for Lemmas 14,301 - 14.341 Definition 15 ,261 is analogous to Definition I4.32| 
and Corollary 15.291 is an analog of Corollary 14.351 

< e n ^2k~ 2l+1 ~ 2S and |e | < 6 "~ 1) k 21 ' 1 ' 5 ^, then the equation 

A (n-1) (y^- 1 )^)) = k 21 + e (5.23) 

has no more than two solutions, (pt ^ n *\ in <P n -i- For any ipf ^ n there is ^ n 2 ^ G <& n ~2 
such that 



Lemma 5.24. If < b Q n 1] 



u±(«-i) 



^±(n-2) |<r (n) /4 



(5.24) 

Lemma 5.25. Suppose < b^ ^ < e n -2k ~ 2l+1 ~ 2S and ip € $ n -i obeys the inequality analogous 
to flTER ): \<p - {n ~ 2) \ < rH Then, 
d 



(5.25) 



Definition 5.26. Let vf {n l) {b^' 1 ^) be the open disks centered at f with radius A n ^; 



±(n-2) 



7 ±(n-i)p„_i)^ be their boundary circles and oi" -1 ^^- 1 )) = r+ (n_1) ur; M) . 

Lemma 5.27. For any <p in <2> n _i \ O^ 1 '^ (b^ 1 ^), 



A (n-D (y(^ {ip ))-k 

Lemma 5.28. For any ip G <2> n _i \ oi n_1) (U n ~^) , 



> to 



(n-l) k 2l-l-6 r (n)_ 



H (n~l) (y(«- 1 )( (/3 ) 



.2/ 



< 



16 



:2i 



-i 



< 



16 • 4 n - 2 c k 2+2s ^ 



(5.26) 

(5.27) 
(5.28) 



Corollary 5.29. J/e n _ifc- 2 '+ 1 - 2 ' 5 < &J n-1) < e n _ 2 A;- 2 ' +1 - 25 , andp belongs to # n _ x \oi n_1) (ft^ 1 )), 



then 



H (n-l) Un-l)^\_ k 2l\ 1 



< 



-2 ' 



.2/ 



< 



4 n - 2 c A; 2+2s ™- 1 



-n-l 



(5.29) 
(5.30) 



The corollary follows from the condition on b^ 1 ^, estimate (|5.1|) and the formula r*™) 
r (n-l)^.-2-4s n -<5_ 
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6. LlMIT-ISOENERGETIC SET AND ElGENFUNCTIONS 

6.1. Limit-Isoenergetic Set and Proof of the Bethe-Sommerfeld Conjecture. At every 
step n we constructed a set B n (A), 23 ra (A) C ® n _i(A) C Si (A), and a function x n (A, v), v G 25 n (A), 
with the following properties. The set D n (A) of vectors x = x n (X,V)V, V G 25 n (A), is a slightly 
distorted circle with holes, see FigEJ Fig0 formula 1)1.13)1 and Lemmas 12.111 13.42| 14.241 15.171 

For any x n (X, v) G D n (A) there is a single eigenvalue of H^ n \x n ) equal to A and given by a 
perturbation series. 15 Let 53 M (A) = (Xi=l ^n{X). Since 23 n _|_i C 23 n for every n, Boo (A) is a unit 
circle with infinite number of holes, more and more holes of smaller and smaller size appearing 
at each step. 

Lemma 6.1. The length of r B 00 (X) satisfies estimate U.9\) with 73 = 5/2. 

Proof. Using (|2.44jl . 1)3. 64|) . 1)4.34)1 and (|5.16|) and consideirng that S n rs 2 n si, we easily conclude 
that L(B„) = (1 + 0(A;- 5 / 2 )), k = X 1 ' 21 uniformly in n. Since B n is a decreasing sequence of 
sets, (HSJ) holds. ■ 

Let us consider x^A, z?) = hnin^oo x n (X, u), v G ^^(A). 

Lemma 6.2. T/ie ^mii Xoo(A, z?) exists for any v G 200(A) and i/ie following estimates hold when 
n > 1: 

|x oc (A, l 7)-x n (A,^)| <4e^-2m 5 en = exp(-^»), s n = 2 n - 1 a 1 . (6.1) 

Corollary 6.3. For every v G Boo (A) estimate ll.lUj) holds, where 
74 = (41 - 3 - 4*1 - 3<5)/2Z > 0. 

The lemma easily follows from the estimates (|3.63(l . ()4.33() and ()5.15|) . To obtain corollary we 
use 1)2.43)1 and take into account that 70 = 21 — 2 — 4si — 26. 

jm=l dip ' 



Estimates ()3.63)) . ()4.33)) and ()5.15)1 justify convergence of the series Ylm=i "75!? j ano ^ nence ) of 



the sequence We denote the limit of this sequence by 

Lemma 6.4. The estimate with 75 = (4/ - 5 — 8s\ — 4d)/2l > 0, holds for any V G "Boo (A). 

We define Doo(A) by 1)1. 8|) . Clearly, 2) 00 (A) is a slightly distorted circle of radius k with infinite 
number of holes. We can assign a tangent vector Q^v + xp,, fl = (— sin ip, cos (p) to the curve 
Doo(A), this tangent vector being the limit of corresponding tangent vectors for curves D n (X) at 
points x n (X, v) as n — > 00. 

Remark 6.5. We easily see from ()6.1)l . that any x G Dqo(A) belongs to the (4e n k~ 2l+1 )- 
neighborhood of D n (A). Applying perturbation formulae for n-th step, we easily obtain that 
there is an eigenvalue X^ n \x) of H^ n \x) satisfying the estimate X^ n \x) = X + 6 n , S n = O (e*)> 
the eigenvalue X^ n \x) being given by a perturbation series of the type 1)3.55)1 . Hence, for every 
x G T> 00 (\) there is a limit: 

lim A (n) (i?) = A. (6.2) 
n— +00 

Theorem 6.6 (Bethe-Sommerfeld Conjecture). The spectrum of operator H contains a semi-axis. 

15 The operator H {n \x) is defined for every 3? £ K 2 as expiained in Remark 13.161 page 122) The perturbation 
series is given by a formuia analogous to 113.551 . which coincides with 13,151 up to a shift of indices corresponding 
to the parallel shift of x into K n . 
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Proof. By Remark I6.5| there is a point of the spectrum of H n in the ^-neighborhood of A for 
every A > kf, k* being introduced by (|4.2|) . Since \\H n — H\\ < e^, there is a point of the spectrum 
of H in the ^-neighborhood of A, <5* = S n + e^. Since it is true for every n and the spectrum of 
H is closed, A is in the spectrum of H. | 

6.2. Generalized Eigenfunctions of H. A plane wave is usually denoted by e l( ~ k ' x \ k E M. 2 . 
Here we use ok instead of k to comply with our previous notations. We show that for every ok in 
a set 

Soo = Ua>a„£ ) 00(A), K=kf 1 
there is a solution ^^(ok, x) of the equation for eigenfunctions: 

(-A) 2 ^oo(x, x) + y(x)^ oc (x, x) = A 00 (x)^ 00 (x, x), (6.3) 

which can be represented in the form 

a:) = e'^' x) (l + u 00 (x,x)), IhooO^a^L^) < c|x| -71 , (6.4) 

where Uoc(ok,x) is a limit-periodic function, 71 = 2/ — 4 — 7si — 2(5 > 0; the eigenvalue Xoo(oi) 
satisfies the asymptotic formula: 

Aoo(i<) = |x| 2Z + 0(|i?i~ 72 ), 72 = 2Z-2-4si-35>0. (6.5) 

We also show that the set Soo satisfies (jl.7|) . 

In fact, by (jHHJ, any ok G D 00(A) belongs to the (e n /c~ 2 ' +1_<5 ) -neighborhood of D n ( A). Applying 
the perturbation formulae proved in the previous sections, we obtain the following inequalities: 

\\E^\ok)-E { - Q \ok)\\ l <2k-^\ \\E^ n+1 \>t)-& n \ok)\\ x <A%el, n > 1, (6.6) 

|AW(i?) _ |x| 2/ | < 2k-~<\ |A (n+1) (x) - A (n) (i<)| < 12et n > 1, (6.7) 

where £( n+1 ), #( n ) are one-dimensional spectral projectors in L2(Q n +i) corresponding to poten- 
tials W n+ i and W n , respectively; \^ n+l \ok) is the eigenvalue corresponding to _E( n+1 )(i<), E^(ok) 
corresponds to V = and the periods a\, 02- This means that for properly chosen eigenfunctions 
ty n+1 (ok,x): 

11*1 " *o|U 2 ( Ql ) < 4£^|Q 1 | 1 / 2 , * (x) = (6.8) 

||* n+1 - *„|U aW „ +l) < lOOe^ [Qn+il 1 / 2 , (6.9) 

where is \I/ n extended quasi-periodically from Q n to Q n +i- Eigenfunctions ^> n , n > 1, are 
chosen to obey two conditions. First, || , I'n||L 2 (Q n ) = \Qn\ l ^ 2 \ 16 second Im(^ n , ^ n _i) = 0, here 
(•,•) is an inner product in L2(Q n )- These two conditions, obviously, provide a unique choice 
of each fy n . Considering that ^ n +i and *$> n satisfy equations for eigenfunctions and taking into 
account (|6.7|) . (|6.9|) we obtain: H^n+i — ^! n \\ W 2i^Q n+i ^ < cfc 2/ e^|Q n+ i| 1//2 , n > 1 and, hence, 

||* n +l - *n||i 00 (Q n+1 ) < c/c 2/ e^|(5n+i| 1/2 - Since * n+1 and § n obey the same quasiperiodic con- 
ditions, the same inequality holds in the whole space M 2 : 

||^n+i - *n\\ Loo (R2) < ck 2l el\Q n+1 \ 1 ' 2 , n > 1, (6.10) 

16 The condition ||*„||l 2 (q„) = \Qn\ 1/2 implies ||*„||l 2 (q„ +i) = \Q n +i\ 1/2 ■ 
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where \l/ n +i, \]/ n are quasiperiodically extended to M. 2 . Obviously, we have a Cauchy sequence in 
L OQ (M 2 ). Let \&oo(x, x) = lim n ^oo ^ n (x, x). This limit is defined pointwise uniformly in x and in 

Theorem 6.7. For every sufficiently large A and ok € D 00 (\) the sequence of functions ^ n (ok, x) 

converges in L OQ (M?) and W?, (JS?). The limit function (ok, x), ^ r 00 (x, x) = limn^oo ^ n (ok, x), 
satisfies the equation 

(-A) 2/ ^ 00 (x, x) + V(x)^ 00 (ok, x) = A^ oc (x, x). (6.11) 
It can be represented in the form 

x) = S*^ (1 + Uoo (ok, x)) , (6.12) 
where Uoo(ok,x) is a limit-periodic function: 



x,x) = ^2<u n (ok,x), (6.13) 



n=l 



u n (ok,x) being a periodic function with the periods 2 Mn 1 b±,2 Mn l b2, 2 Mn ~ /c 2 ™ 81 , 

||«i||r (m <ck~^, 7l = 2Z-4-7 Sl -25>0, (6.14) 



\\un\\ Loo (Ri)<ck 2l e 3 n _ 1 \Q n \ 1 /\ n>2. (6.15) 
Corollary 6.8. Function Uoo(ok,x) obeys the estimate J 6'. 

Remark 6.9. If V is sufficiently smooth, say V € C l (R), then estimate (|6.14|) and, hence, (|6.4() 
can be improved by replacing 71 by 70. 

Proof. Let us show that is a limit-periodic function. Obviously, = ^ o+Y^=o(^ n+i— ^n), 
the series converging in L QO (R 2 ) by ()6.1(J|) . Introducing the notation u n+ i = e -1 ^'^ (*$t n+ i — \J/ n ), 
we arrive at (|6.12|) . q6.13j) . Note that u n is periodic with the periods 2 Mn_1 6i, 2 Mn_1 6 2 - Estimate 
(|6.15j) follows from (|6.1U|) . We check l|6.14j) . Indeed, by (|6.8|1. Fourier coefficients (iti)j, j G Z 2 , 
satisfy the estimate |(iti)j| < cA; -70 ^!! 1 / 2 < c/c~ 7o+Sl . This estimate can be easily improved for 
j : Pj(0) > 2k: \(ui)j\ < c\j\~ 21 . Summarizing these inequalities and taking into account that the 
number of j : Pj(0) < 2k does not exceed cok 2+2si , we arrive at ()6.14|) . It remains to prove ()6.11|) . 
Indeed, ^ n (ok,x), n > 1, satisfy equations for eigenfunctions: H n '^ n = (ok)^ n . Considering 
that ^ n (ok,x) converges to x) in W 2 Uoc and relation (|6.2j) . we arrive at (|6.11|) . ■ 

Theorem 6.10. Formulae \6.3\) . \6.J\) and 16. fy) hold for every ok E Soo- The set Soo is Lebesgue 
measurable and satisfies \1. 7| ) with 73 = 5/2. 

Proof. By Theorem HT71 flEfl, (jH31) hold, where X^k) = A for x G D 00(A). Using (|TTT7|l . 
which is proven in Corollary E31 with Xqo = \ok\, we easily obtain (|6.5j) . It remains to prove 
p.7|) . Let us consider a small region f/„(Ao) around an isoenergetic surface D„(Ao), Ao > kf . 
Namely, U n (X ) = U| A _ Ao | <rj ^n(A), r n — (-n—ik 2 ^ , k — Aq^ . Considering an estimate of 
the type (|3,61|) for X^ n \ok), which holds in the (e„_ifc~ 2i+1 ~ 2<5 )-neighborhood of !D„(Ao), we see 
that U n (\o) is an open set (a distorted ring with holes) and the width of the ring is of order 
e n -ik- 2l+1 - 25 . Hence, |J7 n (A )| = 2nkr n (l + o(£r 5 / 2 )). It easily follows from Lemma EH] that 
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U n+ i C U n . Definition of Doo(Ao) yield: Doo(A ) = n^L 1 [/„(A ). Hence, 9oo = n£° =1 3 n , where 
Sn = UA>A*I'n(A). Considering that U n +i C U n for every Ao > A*, we obtain 9n+i C 9n- 
Hence, |9oc>nBR,| = limn-^ |Sn H BrJ. Summarizing volumes of the regions U n , we easily 
conclude |9„nB R | = |B R | (l + 0(i?-' 5 / 2 )) uniformly in n. Thus, we have obtained 1)1. 7JI with 
73 = 5/2- I 

Appendix A. Appendices. 
A.l. Proof of Corollarv l3.6l Letro€x2- Taking into account the relation A^(ro+27rj»/A r ia) = 
k 21 and the definition of C 2 , we see that \\f\ro + 2tt P / N x a) - z\ = e x /2. Using and 
the last equality, we easily obtain: \\n\ro + 2irp/Nia) — z\ > ei/2 for A^ (to + 2irp/N\a) 7^ 
\j (tq + 2irp/Nia). Therefore, for any z € C 2 , 

||(#( 1 )(r )-z)- 1 ||<2/ ei , (A.l) 

i.e., (f.-t. llj> is proved for To € X2- Now we consider r in the complex 
(e\k~ 2l+1 ~ s )— neighborhood of X2- By Hilbert relation, 

(#W(r) - z)- 1 = (# (1) (r ) - z)- 1 + T^H^ir) - z)~\ 
T x = (H (t ) - zrHHoiro) + k 21 ), T 2 = (H (t ) + k 21 )' 1 (H (r ) - H (r)). 

Suppose we have checked that HT1T2H < k~ s . Then, using (jA.lj) and a standard argument as in 
Lemma 12.91 we easily arrive at p.llj) . The estimate < 4k 21 /e\, easily follows from (|A.1|) . 
The estimate ||T 2 || < 2Zei£r 2 ^ 5 easily follows from |r - r | < eik- 2l+1 - & . Thus, ||TiT 2 || < 8lk~ s 
and, hence, (|3.11|) is proved. Estimate ()3.12j) follows from the Hilbert relation 

(i?«(r) - z)- 1 = (Ho(r) + k 21 )' 1 + (H (r) + k 2l )-\-Wi + k 21 + z){H^(r) - z)'\ 

an elementary estimate \\{Hq{t) + k 21 )' 1 ^ < hb 2 k~ 2l+2+2s2 and (|5~TT|) . 

A. 2. Proof of estimate (|3.50|) . Note that 

\V6(<P)+Pm@)\l ~ \V(<P)+Pm(0)\? 



\yo(<p) +Pm(0)\l - Wo(<p) + hx(tp) +p m (0)\l =0[\hi\- \ya(<p) + p m (0)\ + \hi 



(A.2) 



If p m (0) < 4k, then, obviously, \yo(<p) + Pm(®)\1 = 0(k 2 ). Considering (|2.43|) . we arrive at 



\vo(<p) +p m (o)\? - W&) +Pm(o)\l l = o(\h\ ■ k 21 - 1 ) = 0{k 

Combining this relation with (|3.44|) where eo = k~ 4 ~® Sl ~ 3S , we obtain 
+ftn(0)|? " \y(<p) +ftn(0)|? 



,2/-2-4si-27o-5n 



|y ^)+Pm(o)|?-^ 

c\hx\k 2l ~ l 



< 



£2Z-9-12si-6<5 



O(6 - 1 |/ ll |fc 8+12si+6<5 ) = O(6 - 1 A; 7 + 8si - 2 ^+ 5 ' 5 ) = o(l), 



(ill 
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the last estimate following from the restriction on bo in the statement of the lemma. Thus, Q3.5UJ) 
holds for p m (0) < 4k. Suppose p m (0) > 4k. Using (|3.45j) and (|A.2|) . we easily obtain: 



M<P) +Pm(0)|f - \y(<p) +p m (0)\ 



21 



\M<P)+Pm(0)\ 2 J-k 21 

This means (|3.5Uj) holds in this case too. 



< 



Pm(0) 



0(k 



-2-4si-2 7o -<5^ 



A.3. Proof of estimate (f3"U5|). We use (j3~B7|) . First, we represent VAW(xi) - X7\W(a, x 2 ) as 
the sum of and ^4 2 : 

Ai = VA (1) (xiz7) - VA (2) (a,xiz7), A 2 = VA (2) (a, xiz7) - VA (2) (a, x 2 i7). 

Considering ()3.57|) . we obtain \A\\ < 24aef k 2l ~ 1+s . Using the Mean Value Theorem and (|3.27j) . 
we get: |-A 2 | < c/c 2 '~ 2 |x 2 — jt\\. Applying H3.6Uj) yields: |A 2 | < cae^A; -1 . Adding the estimates for 
A\ and A2 we arrive at: 



|VA (1) (xii7) - VA (2) (a,x 2 i?)| < 25ae?fc 



3r,2i-l+<5 



Similarly, 



|xiVA (1) (xii7) - x 2 VA (2) (a,x 2 z7)| < 2Qae\k 



3 tJ2/+<5 



(A.3) 



(A.4) 



■270+1- 



. Substituting (|A.3|) . (|A.4|) and the last estimate into 



By Lemma EH1 % = '% = 0(Af 
the formula for F(X, u), we obtain 

A.4. Lemma lA.lt 

Lemma A.l. Ift is in the complex {2k~ l ~ Asl ~ 26 )-neighborhood of the nonresonance set Xi(A, 6), 
then the operator-function (H^ (t) — z) 



,-l-4si-2<5-j 

has a single pole z$ inside the contour C\ and 
-1 



(z-zq) (H^(t)- Z y 



< 16, 



(A.5) 



for every z inside C\. 
Proof. By 



det 



< 2l|Wi| 



(A.6) 



H (t) - z 

for every z on the contour G\ . Using the estimate (|2.1Uj) and considering that Hq (t) is a diagonal 
operator, we easily obtain ||(i?o(i) — z)~ || < k~ 2l+4:+6si+s . It immediately follows that the right- 
hand part of ()A.6|) is less then 1. Applying Rouche's theorem, we conclude that the determinant 
has the same number of zeros and poles inside C\. Considering that (-ffo(^) ~ z )~ has a single 
pole, we obtain that (t) - z) has a single pole too. We denote it by Zq. Obviously, 

(z — Zq) (H^(t) — z) 1 is holomorphic inside C\. The definition of C±, given in Corollary 12.21 
yields: | z — zq \ < 2k 2l ~ 2 ~ 4si ~ s . Using the last inequality and estimate (|2.30j) . we obtain (jA.5|) for 
any z € C\. By the Maximum principle it can be extended inside C\. | 
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A. 5. Proof of Lemma 13.481 By Lemma 12.111 Part 2, the function >t\{ip) is holomorphic in the 
(/c~ 2 ~ 4si ~ 2<5 ) — neighborhood of <P\ and A^(i<i ((f)) = k 21 . Hence, the equation (|<3.92[) is equivalent 
to A (1) (?/((/?)) = A (1) (y(</?) - b) + e . We use perturbation formula (|2.34|) : \y{^)\f + fi(y(<p)) = 
\y(ip) — b\ 21 + fi{y{f) — b) + eo- This equation can be rewritten as 



(2< y », b% - 16| 2 ) (\m\f~ 2 + ■■■ + \m - b\ 2 j- 2 



+ fM<p))-fMv)-b)-e = 0, (A.7) 

where (y, 6)* = y\bi +2/2^2- Using the notation b = bo (cos ipi,, sin ip^), dividing both sides of the 

equation (|A.7|) by 2bok(\y(ip)\ 2l ~ 2 H h |y(y>) - b\f~ 2 ^J , and considering that y((p) = j?i(^) + 6 = 

(A; + hi)v + b, we obtain: 

cos(<£ - <^ 6 ) - e c/i(v?) + 32 (v) = °' ( A - 8 ) 

-1 



where gi(<p) = 2b k( \y(f)\ 



\2l-2 



+ 



+ \y(<p) - b\i 



21-2 



and 



92 ^ = ^bdk^* + ^ + (MvW) ~ /i(^) ~ k))#l(^)> ^i(v) = hi((p)P. 

Obviously gi = O (b^ 1 k~ 2l+1 ) . Let us estimate 32 ((f). Using the inequality (|2.43|) for hi, and 
considering that bo = o(k~ 1 ~ 4si ~ 2S ), we easily obtain: 



(hi(cp),b), 



bok 



<^ = 0(k- 2 - is ^- 5 ), ^ = o(k~ 2 -^- 2S ) 



By oa . 

/i(ifoO) - AGfa) - 6) I < sup 1 v/i|6o = oibok 21 - 1 - 2 ^ 5 ) = O(b k- 2l+3 + 8s ^ 5S ), 

21 - 11 



and therefore, 



{fM<p))-fM<p)-b))9i(<p) = 0(k-^+ Ss ^ s ). Since Sl 



32 



wc 



have g 2 (<p) = o(k- 2 - 4si ~ 2S ). Using gi (<p) = 0{b Q l k- 2l+1 ) and e < & ^~ 3 ~ 4si ~ 3<5 , we obtain 



-2-4S1-25 



!-^). Thus, 



7T 



32 M - eo3l(y) = °( k 



-2-4si-2i5 



(A.9) 



Suppose that db — is in the 2 4si 2<5 ) -neighborhood of <^i. We draw two circles C± 

centered at (fj, ± ^ with the radius |fc~ 2 ~ 4si ~ 2<5 . They are both inside the complex l^- 2 - 45 !- 2 " 5 . 

neighborhood of <Pi, the perturbation series converging and the estimate (|A.9j) holds. For any 92 on 

C±, \(p-(tp b ±7r/2) 



l^-2-4si-2<5 



1 25 and, therefore, | cos(p-<#,)| > ±k 2 4si 2,5 > |32(</>)- £ o3iMI 



for any </? G C±. By Rouche's Theorem, there is only one solution of the equation 



inside each 



G±. If^+vr/2isnotinthe(i/c- 2 - 4si - 25 )-neighborhoodof^i,then|cos(v9-99 fe )| > ±k~ 2 ~ 4si ~ 2S 



in <P\ and, hence, equation ljA.81) has no solution. Thus, there are at most two solutions in <P\ 
and - (<p b ± tt/2)| < ifc- 2 " 4 ^- 25 . 
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A. 6. Proof of Lemma I3.49L Using the perturbation formula (|2.34j) . we obtain: 







a 



d_ 

dip 




AW(^))-A«(^)-6) 



V\V&)\? - ^W) - b\l l , {k + h l )fi + h' l u) + 
where V = (cos p, sin (p) and jl = V = (— sin p, cos p). Note that 



v) 



(A.10) 



V|jfoO|? " V|yM - b\l l = 2l\y(p)\ 2 ! *M ~ 2l\y(p) - b\ 2 J \y{p) - b) 



2l\y(p)\l l 2 b + 2l \y(p)\l 



\2l-2 



\y(<p) 



\2l-2 



{y(p)-b). (A.11) 



Substituting (|A~Tll into (|A~10|) . we get ^aJ 1} (y(^)) = Ti + T 2 + T 3 , 



21 



\v(p)\* 



21-2 



|y(v) 



■|2Z-2 



y(p) -b,(k + h 1 )fl+h , 1 i7) , 



T 3 = (Vh{y(<p)) - Vh{y(<p) -b),(k + h l )fl+h' l v) . 



Considering that p is close to pb±n/2, we readily obtain: (b, V)* = o(bo), {b,fi)* = ±&o(l + o(l)). 
Using also estimates l|2.43|) for hi, we get 7\ = ±2Z& & 2 ' _1 (1 + o(l)). Substituting y((p)-b= (k + 
h\{p))v into the formula for T 2 and taking into account that (p,, z?)* = 0, we arrive at the estimate 



T 2 = o^ok 21 - 1 ). By (l2~3fill with 



m 



V/i(y») - V/i(y(^) - b) = O(& fc- 2 '+ 4 + 12s i+ 75 ) 



Hence, T 3 = o (&o^ 2i X )- Adding the estimates for Ti,T2,T 3 , we get (|H.94|) . 
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